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Abstract 

Vehicles  entering  planetary  atmospheres  at  high  speed  reqnire  an  ablative  heat 
shield  in  order  to  withstand  the  high  thermal  energy  flnx  to  the  body.  The  interaction 
between  the  ablative  prodncts  and  the  flow  held  is  not  well  characterized.  Nnmerical 
simnlations  were  condncted  to  investigate  the  inhnence  of  carbon  ablation  on  shock 
layer  radiation.  Data  collected  from  experiments  performed  in  the  X-2  expansion 
tnnnel  at  the  University  of  Qneensland  was  nsed  to  compare  to  the  simnlations.  The 
model  was  a  short  half-cylinder  made  of  isomolded  graphite  and  was  tested  in  8.6 
km/s  Earth  entry  how.  The  model  snrface  was  heated  within  a  temperatnre  range 
of  1770-3280  K.  The  radiation  emitted  from  the  CN  violet  bands  was  measnred  by 
nltraviolet  spectrometry  in  a  spectral  range  from  353-391  nm.  This  research  devel¬ 
ops  a  novel  hnite-rate  snrface  kinetic  model  for  determining  the  chemical  state  of 
an  ablating  bonndary  layer.  The  proposed  ablation  model  acconnts  for  competing 
snrface  reaction  processes  snch  as  adsorption/desorption,  Eley-Rideal  mechanisms, 
oxidation,  nitridation,  and  snblimation.  The  inclnded  oxidation  mechanisms  predict 
CO  as  the  primary  oxidized  prodnct  at  the  considered  snrface  temperatnres,  which  is 
in  agreement  with  experiment  and  theory.  A  previons  model  had  incorrectly  predicted 
CO2  as  the  primary  oxidized  prodnct  for  a  majority  of  the  tested  snrface  tempera¬ 
tnres.  The  ablative  gas  species  predicted  by  this  new  snrface  model  resnlts  in  better 
agreement  with  experimental  spectral  measnrements  than  predictions  provided  by 
legacy  ablation  models,  and  represents  a  signihcant  improvement  in  cnrrent  modeling 
capabilities  for  hypersonic  noneqnilibrinm  ablating  re-entry  flows. 
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A  Nonequilibrium  Finite-Rate  Carbon  Ablation  Model  For  Radiating 

Earth  Re-entry  Flows 

I.  Introduction 

f  I  '^HE  Air  Force  confronts  an  nncertain,  dynamic  fntnre  as  it  faces  strategic  pol¬ 
icy  changes,  force  strnctnre  overhanls,  technological  advancements,  and  severe 
resonrce  constraints  that  in  2010  prompted  the  United  States  Air  Force  Chief  Scien¬ 
tist  to  release  his  vision  for  fntnre  science  and  technology  (S&T)  investments.  The 
“Technology  Horizons”  vision  docnment  lays  the  fonndation  for  how  the  Air  Force 
was  going  to  maintain  its  technological  advantage  thronghont  the  world.  Within  this 
docnment,  the  Chief  Scientist  identihes  30  Potential  Capability  Areas  (PCAs)  that 
can  connter  different  potential  threat  environments  and  keep  the  United  States  at  the 
technological  forefront  [2].  Among  the  potential  threat  environments  the  Air  Force 
will  likely  enconnter  in  the  fntnre  is  that  of  an  anti-access/area-denial  environment. 
One  of  the  key  technology  areas  that  wonld  enable  the  Air  Force  to  operate  in  snch 
an  environment  is  that  of  hypersonic  systems.  In  fact,  fonr  of  the  30  PCAs  identify 
hypersonic  systems  or  technologies  as  being  critical  capabilities  to  ensnre  operational 
snccess,  which  are  prompt  theater-range  ISR/strike  systems,  penetrating  persistent 
long-range  strike,  high-speed  penetrating  crnise  missile,  and  rensable  airbreathing 
access-to-space  lannch  [2]. 

The  need  to  penetrate  highly-defended  adversary  airspace  to  achieve  desired  ef¬ 
fects  will  likely  remain  an  endnring  reqnirement  past  the  intended  vision  time  frame 
of  2030.  The  Air  Force  S&T  leadership  reinforced  this  vision,  led  by  the  Air  Force 
Research  Laboratory  (AFRL)  commander,  by  making  long  range  precision  strike  an 
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increased  emphasis  technology  area  [3],  Investment  areas  inclnde  high-speed  propnl- 
sion  systems,  high  temperatnre  materials  and  strnctnres,  and  high-speed  system  in¬ 
tegration  and  demonstration.  The  goal  of  long  range  precision  strike  is  to  develop 
a  high-speed  weapon  or  aircraft  to  engage  high-valne,  time  sensitive  targets  in  anti¬ 
access/ area-denial  environments. 

Development  of  hypersonic  systems  encompasses  a  mnltitnde  of  scientihc  disci¬ 
plines,  bnt  the  design  of  a  thermal  protection  system  (TPS)  is  one  of  the  most  crit¬ 
ical  investment  areas  becanse  its  failnre  conld  resnlt  in  a  total  loss  of  the  vehicle. 
Becanse  a  TPS  is  a  critical  investment  area,  the  Air  Force  had  a  joint  program  with 
the  Defense  Advanced  Research  Projects  Agency  (DARPA)  called  Falcon  aimed  at 
advancing  capabilities  for  access  to  space  and  hypersonic  flight  [121,  122].  As  part 
of  the  Falcon  program,  a  series  of  hypersonic  demonstration  vehicles  were  developed, 
and  specihcally  the  vehicle  designated  as  Hypersonic  Technology  Vehicle  2  (HTV- 
2)  was  snpposed  to  demonstrate  enabling  technologies  for  fntnre  hypersonic  re-entry 
operational  systems  [30].  The  technologies  to  be  evalnated  were  thermal  protection 
systems,  aerodynamic  shapes,  manenverability,  and  long-range  commnnication  for 
hypersonic  crnise  and  re-entry  applications.  The  design  of  the  external  TPS  proved 
to  be  one  of  the  most  challenging  technology  areas  for  the  HTV-2  program.  Unfor- 
tnnately,  the  HTV-2  had  two  nnsnccessfnl  test  flights  that  resnlted  in  an  incomplete 
evalnation  of  pertinent  hypersonic  technologies.  Other  high-speed  programs  snch  as 
the  Space  Shnttle,  NASP,  X-33,  Genesis,  X-37,  and  Hy-Fly  also  spent  a  considerable 
amonnt  of  effort  designing  effective  TPS  systems  [86].  All  of  these  vehicles  ntilized 
carbon-based  materials  for  their  critical  high  temperatnre  aerodynamic,  strnctnral, 
and  propnlsion  applications. 

One  of  the  reasons  why  a  TPS  is  a  challenging  technology  area  is  becanse  a  TPS 
is  snbjected  to  severe  thermal  and  mechanical  loads  when  exposed  to  hypersonic  re- 
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entry  environments  and  mnst  be  designed  to  prevent  excessive  heat  from  damaging 
the  vehicle.  The  materials  nsed  for  a  TPS  interact  with  the  flow  throngh  varions 
thermochemical  and  thermophysical  processes  snch  as  ablation,  spallation,  thermal 
condnction,  and  radiative  transport.  In  most  flight  regimes,  radiation  is  ignored  be- 
canse  the  flow  does  not  reach  very  high  temperatnres.  However  for  hypersonic  flows, 
strong  shocks  develop  that  lead  to  high  temperatnres  aronnd  the  vehicle.  To  ac- 
cnrately  analyze  the  overall  heat  transfer  in  these  flows,  modeling  of  air  chemistry 
effects  and  radiation  from  hot  shock  layers  mnst  be  considered.  Additionally,  re-entry 
conditions  can  have  thermal  and  chemical  noneqnilibrinm  within  different  regions  of 
the  flow  and  these  relaxation  processes  mnst  be  acconnted  for  as  well.  In  recent  years, 
compntational  flnid  dynamics  (CFD)  capabilities  have  made  great  progress  in  simn- 
lating  the  effects  of  ablation  and  radiation.  However,  the  degree  of  complexity  and 
hdelity  varies  within  different  research  codes  in  how  to  model  ablation  and  radiation 
since  the  compntational  cost  of  implementing  the  most  general  theories  is  prohibitive 
and  not  practicable.  In  reality,  all  of  these  phenomena  are  highly  conpled  [14]. 

One  consideration  that  is  vital  to  developing  a  TPS  is  the  accnrate  modeling  of  gas- 
snrface  interactions  nnder  hypersonic  conditions,  which  remains  to  be  a  challenging 
and  complex  problem.  Snrface  reactions  both  degrade  TPS  materials  and  contribnte 
to  aerothermal  heating.  It  is  crncial  to  have  reliable  predictive  capabilities  that  can 
accnrately  compnte  the  thermophysical  environment  snrronnding  the  vehicle  as  well 
as  the  chemical  kinetics  occnrring  on  the  snrface.  However,  detailed  information 
regarding  gas-snrface  interactions  is  lacking  from  experimental  data  and  the  ablation 
process  is  mainly  characterized  by  recession  rates  or  reaction  efficiencies  [28,  33,  60, 
95,  114,  116,  135].  Therefore,  nnmerical  models  nsed  for  the  snrface  kinetic  processes 
are  of  low  hdelity  and  large  safety  factors  are  applied  to  TPS  designs.  For  example, 
the  prediction  of  the  convective  heat  load  on  the  Mars  Science  Laboratory  developed 
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Figure  1.  Surface  mass  flow  rates  from  a  collection  of  graphite  ablation  experiments  as 
function  of  surface  temperature  at  1  atm.  Adapted  from  Gosse  et  al.  [36] 

by  NASA  had  a  60%  uncertainty,  which  required  a  40%  increase  in  TPS  thickness 
[18],  Additionally,  for  strongly  radiating  species  like  CN  [118],  the  accurate  modeling 
of  species  formation  rates  is  critical  for  calculations  of  radiative  heating.  Greendyke 
and  Hartung  [38]  noted  that  the  radiative  contribution  to  the  overall  heat  transfer 
for  blunt  re-entry  vehicles  could  be  as  high  as  50%  or  more  in  high  nonequilibrium 
flow  environments. 

To  improve  modeling  efforts,  the  characterization  of  TPS  materials  must  be  fully 
understood.  The  fundamental  material  found  in  most  carbon-based  TPS  is  graphite, 
which,  in  hypersonic  re-entry  flows,  can  react  due  to  oxidation,  nitridation,  and  sub¬ 
limation.  The  dominant  mechanisms  of  graphite  ablation  over  a  range  of  surface 
temperatures  at  1  atm  are  shown  in  Figure  1.  Data  from  different  graphite  ablation 
experiments  is  compared  to  equilibrium  ablation  results  predicted  using  the  NASA 
Chemical  Equilibrium  with  Applications  (CEA)  thermo-chemical  database  [70].  Eor 
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temperatures  less  than  approximately  2000  K,  the  primary  surface  reaction  is  carbon 
oxidation.  Between  temperatures  of  approximately  2000  and  3000  K,  the  ablation 
rate  plateaus  indicating  a  diffusion-limited  ablation  regime.  In  this  regime,  all  of  the 
available  oxygen  at  the  surface  is  entirely  consumed  by  the  oxidation  reactions.  Above 
3000  K,  graphite  sublimation  is  the  dominant  mechanism  and  ablation  rates  expo¬ 
nentially  increase.  Oxidation  and  nitridation  are  exothermic  reactions  that  consume 
solid  carbon  and  inject  gas  into  the  boundary  layer.  The  coupling  between  the  surface 
and  the  flow  held  is  a  non-linear  process  under  hypersonic  hight  conditions,  and  the 
balance  between  catalytic  and  surface-participating  reactions  are  often  unclear  [67]. 
A  correct  understanding  and  accurate  modeling  of  all  of  these  ablative  phenomena 
play  an  integral  part  in  the  design  of  TPS  for  re-entry  vehicles. 

Ablation  is  typically  modeled  under  the  assumption  that  the  surface  and  surround¬ 
ing  gas  are  in  chemical  equilibrium.  Previous  comparisons  to  experiments  and  hight 
data  have  shown  that  assumptions  of  equilibrium  provide  a  conservative  estimate  of 
the  heat  huxes  on  the  surface  of  vehicles  such  as  the  Apollo  capsule.  Pioneer- Venus 
probes,  Galileo  entry  probe,  and  Stardust  return  capsule  [98,  49,  112].  For  the  Star¬ 
dust  return  capsule  that  had  an  Earth  entry  velocity  of  12  km/s,  equilibrium  surface 
recession  was  over  predicted  by  50%  at  some  locations  [112].  When  the  analysis  was 
performed  assuming  nonequilibrium  at  the  surface,  peak  ablation  rates  decreased  by 
39%  and  the  total  heat  load  decreased  by  32%.  The  nonequilibrium  surface  reces¬ 
sion  was  closer  to  the  measured  recession  and  was  only  over  predicted  by  11%  [13]. 
Nonequilibrium  had  a  signihcant  effect  on  the  prediction  of  the  thermal  and  chemical 
environment  surrounding  the  Stardust  return  capsule.  For  very  high  energy  flows,  a 
nonequilibrium  treatment  using  hnite-rate  surface  chemistry  models  is  required. 

Another  challenge  to  TPS  design  is  that  there  is  no  universally  supported  surface 
kinetic  model  for  the  air-carbon  system.  There  is  a  vast  amount  of  reference  material 
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on  the  chemistry  of  a  carbon  snrface  [17,  16,  33,  106,  15,  89,  8,  48],  and  in  fact, 
the  hrst  theory  on  snrface  adsorption  of  gases  is  attribnted  to  Langmnir  [51]  back  in 
1918.  There  have  been  several  attempts  to  model  the  snrface  chemical  kinetics,  bnt 
investigators  have  not  reached  a  common  opinion  abont  what  the  main  processes  are 
on  the  carbon  snrface  [136,  44,  14],  As  a  resnlt,  mnch  nncertainty  remains  on  this 
topic.  Generally,  experiments  investigating  these  types  of  snrface  reactions  will  report 
loss  coefficients  describing  the  fraction  of  gas-snrface  collisions  that  remove  species 
from  the  gas  phase  [37,  9,  66].  Usnally  these  loss  coefficients  are  cnrve-ht  as  fnnctions 
of  temperatnre,  which  can  then  be  incorporated  as  a  bonndary  condition  into  CFD 
simnlations  [113].  However,  this  approach  is  pnrely  empirical  and  does  not  acconnt 
for  any  physics-based  mechanisms  cansing  the  snrface  reaction  [67]. 

It  is  important  to  note  that  the  radiation  emitted  away  from  a  re-entry  vehicle  is 
of  interest  to  not  only  design  engineers  bnt  also  the  intelligence  commnnity  (IC) .  It  is 
reasonable  to  assnme  that  near-peer  conntries  to  the  United  States  (US)  are  pnrsning 
stealthy  hypersonic  weapons  and/or  conntermeasnres  to  hypersonic  technologies  in 
response  to  onr  own  system  developments.  Therefore,  it  is  imperative  that  the  IC  have 
the  capability  to  detect  and  identify  enemy  hypersonic  weapons  as  a  matter  of  national 
secnrity.  As  Martin  [69]  points  ont,  while  it  may  be  possible  to  rednce  radar  cross- 
sections  and  mask  propnlsion  signatnres,  it  is  impossible  to  conceal  radiation  emitted 
from  vehicles  moving  at  high  Mach  nnmbers.  Hence,  any  improved  thermochemical 
models  will  enhance  the  capability  of  the  IC  to  correctly  identify  hostile  vehicles  that 
pose  a  threat  to  the  US  or  its  interests.  It  is  dne  to  all  of  these  considerations  that 
the  model  which  is  developed  in  the  following  chapters  is  proposed. 

The  goal  of  the  research  described  in  this  dissertation  is  to  advance  onr  nnder- 
standing  of  the  gas-snrface  interactions  for  the  air-carbon  chemical  system  in  a  high 
temperatnre,  reacting  flow  environment  by  improving  cnrrent  ablation  models  to  bet- 
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ter  match  experimental  measurements.  To  evaluate  the  ablation  models,  a  series  of 
high-speed  flow  experiments  were  conducted  in  the  X-2  facility  at  the  University  of 
Queensland,  Australia.  The  X-2  facility  is  a  shock  expansion  tunnel  that  can  be 
used  to  test  subscale  models  at  realistic  flight  temperatures  and  enthalpies.  The 
model  used  in  the  experiments  was  a  short  half-cylinder  made  of  isomolded  graphite 
and  was  tested  at  8.6  km/s  Earth  entry  flow  monitored  by  ultraviolet  (UV)  spec¬ 
trometry.  The  experiments  pre-heated  the  model  to  high  temperatures  to  stimulate 
surface  reactions  and  increase  ablation  during  the  microseconds  of  available  test  time 
[133].  Further  development  of  pre-heating  techniques  and  capabilities  have  enabled 
the  testing  of  carbon  models  approaching  sublimation  surface  temperatures  [56] .  The 
nonequilibrium  chemistry  occurring  in  the  shock  layer  was  investigated  by  making 
calibrated  measurements  of  the  radiance  emitted  by  the  CN  molecule.  The  choice  of 
analyzing  CN  was  ideal  since  it  is  a  strong  emitter  and  radiates  in  the  UV  spectrum 
for  which  the  X-2  facility  is  equipped  to  measure.  Also,  an  attempt  was  made  to 
capture  radiation  spectra  from  C2  and  C3  at  near  sublimation  conditions.  The  ex¬ 
periments  provided  benchmark  data  to  test  the  validity  of  the  physical  models  used 
in  the  numerical  simulations. 

The  work  outlined  in  this  dissertation  aligns  itself  with  the  following  four  research 
objectives: 

Research  Objective  1:  Perform  an  in-depth  computational  study  by 
comparing  current  state-of-the-art  ablation  models  and  identify  a 
preferred  model  at  conditions  of  interest. 

Research  Objective  2:  Propose  improved  ablation  model  by  developing 
original  finite-rate  surface  reactions  based  on  experiments  and  theory. 

Research  Objective  3:  Perform  sensitivity  study  to  determine  which 
surface  reactions  have  most  impact  on  radiative  heating  and  ablation  rates. 

Research  Objective  4:  Validate  developed  ablation  model  against  the 
benchmark  X-2  experimental  results. 
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Chapter  II  provides  an  overview  of  the  experiments  performed  in  the  X-2  facility.  In 
Chapter  III,  the  theory  and  compntational  methodology  is  presented  describing  the 
implementation  of  the  ablation  bonndary  condition  into  a  CFD  code.  Before  applying 
the  ablation  models  to  the  X-2  experiments,  Chapter  IV  characterizes  the  models 
in  a  representative  binnt  body  re-entry  flow.  Chapter  V  then  presents  the  resnits 
and  analysis  of  the  radiation  comparisons  to  the  X-2  measnrements.  Additionally, 
proposed  modifications  to  a  cnrrent  snrface  model  are  implemented,  providing  better 
agreement  with  the  experimental  data.  Finally,  Chapter  VI  provides  concinsions  and 
proposes  fntnre  work. 
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II.  Experimental  Testing 


Understanding  hypersonic  re-entry  remains  one  of  the  greatest  challenges  to  aerospace 
engineers  dne  to  the  many  different  types  of  phenomena  present  in  the  flow  field,  which 
inclnde  viscons-inviscid  interactions,  snrface  ablation,  noneqnilibrinm  thermodynam¬ 
ics,  finite-rate  chemical  reactions,  and  radiative  heat  transfer  [108],  The  only  way  to 
investigate  these  phenomena  in  a  conpled  manner  is  throngh  flight  testing  becanse 
gronnd  test  facilities  are  nnable  to  replicate  all  aspects  of  hypersonic  flight  conditions. 
However,  the  expense  of  flight  test  programs  has  become  so  great  that  more  emphasis 
is  being  placed  on  gronnd  testing  for  the  design  of  fntnre  hypersonic  systems. 

The  premise  of  this  dissertation  is  that  integrating  compntational  flnid  dynamic 
(CFD)  simnlations  with  gronnd  experiments  will  greatly  improve  onr  nnderstanding 
of  the  phenomena  seen  in  hypersonic  flight  and  mitigate  some  of  the  limitations  of 
gronnd  testing.  For  example,  shock  expansion  tnnnels  can  reprodnce  hypersonic  flight 
speeds  and  enthalpy  conditions,  bnt  the  test  times  are  very  short,  100  ps  to  2  ms, 
which  does  not  allow  for  any  model  thermal  response  [132].  Arc  heated  facilities 
enable  steady-state  testing  of  the  thermal  environment  bnt  do  not  correctly  simnlate 
hypersonic  conditions  [59].  However,  by  attempting  to  piece  together  the  different 
regimes  of  hypersonic  flight  in  gronnd  test  facilities,  and  conpling  those  resnlts  with 
CFD  simnlations,  heavy  reliance  on  flight  test  programs  may  not  be  needed. 

To  stndy  the  effects  of  ablation  in  expansion  tnnnels,  a  recent  pre-heating  concept 
was  developed  at  the  University  of  Qneensland  [132,  133,  56].  The  concept  involves 
electrically  pre-heating  carbon-based  materials,  snch  as  those  commonly  nsed  for  a 
TPS,  np  to  temperatnres  approaching  the  snblimation  ablation  regime.  An  environ¬ 
ment  is  created  where  the  snrface  bonndary  condition  is  matched  to  that  seen  in 
re-entry  flight,  and  the  gas-snrface  interaction  can  be  simnlated.  It  is  then  possible  to 
nnmerically  simnlate  these  flows  and  compare  radiation  intensities  to  those  measnred 
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in  the  experiments.  The  following  sections  describe  the  experimental  setnp  nsed  for 
the  expansion  tnnnel  tests. 

2.1  X-2  Facility 

The  hypersonic  testing  was  condncted  in  the  X-2  shock  expansion  tnnnel  at  the 
University  of  Qneensland,  Anstralia.  The  tnnnel  is  a  free  piston  driven  facility  capable 
of  simnlating  re-entry  conditions  for  different  atmospheres  at  speeds  of  np  to  13  km/s. 
The  operation  of  the  X-2  facility  is  illnstrated  in  Fignre  2.  The  tnnnel  is  a  single 
tnbe  separated  into  three  sections  by  two  diaphragms.  A  2.0  mm  thick  pre-scored, 
steel  primary  diaphragm  separates  the  driver  section  and  shock  tnbe,  while  a  thin 
alnminnm  sheet  secondary  diaphragm  separates  the  shock  tnbe  from  the  acceleration 
tnbe  [29]. 

The  piston  is  accelerated  by  a  large  gas  reservoir  that  compresses  the  driver  gas 
nntil  the  pressnre  is  snfficiently  high  to  rnptnre  the  primary  diaphragm.  The  shock 
wave  that  is  prodnced  propagates  throngh  the  shock  tnbe  accelerating  the  test  gas 
nntil  the  secondary  diaphragm  is  rnptnred.  The  test  gas  is  then  free  to  flow  down¬ 
stream  into  the  lower  pressnre  of  the  acceleration  tnbe.  The  test  gas  expands  and 
accelerates  as  it  travels  throngh  the  nozzle,  greatly  increasing  its  energy  in  an  isen- 
tropic  nnsteady  expansion  [29].  This  whole  process  creates  high  enthalpy  test  flows 
snitable  for  aerodynamic  testing. 

A  steady  test  flow  is  established  throngh  the  nozzle  exit  and  is  passed  over  a 
snbscale  model  or  instrnmentation  package  in  the  test  section.  A  typical  available 
test  time  is  on  the  order  of  tens  to  several  hnndred  microseconds.  Optical  windows 
at  the  sides  of  the  test  section  allow  non-invasive  optical  diagnostic  techniqnes  snch 
as  flow  visnalization  throngh  high-speed  imaging  or  emission  spectroscopy.  Fnrther 
information  on  the  X-2  facility  and  its  capabilities  can  be  fonnd  in  References  [29,  71, 
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s  =  shock  cs  =  contact  surface  usx  =  unsteady  expansion 

rs  =  reflected  shock  sx  =  steady  expansion  «  +  a  =  reflected  u  +  a  cliaracteristic 


Figure  2.  Schematic  diagram  of  the  X-2  expansion  tunnel  and  ideal  x-t  diagram  [57]. 
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Table  1.  Calculated  X-2  freestream  conditions. 


Parameter 

Value 

Total  Enthalpy  (MJ/kg) 

38.2 

Pressure  (Pa) 

847 

Temperature  (K) 

2040 

Velocity  (m/s) 

8500 

Density  (kg/m^) 

1.45  X  10-3 

2/N2 

0.751 

2/O2 

0.225 

Vno 

8.53  X  10-3 

yco 

1.26  X  10-3 

yco2 

4.65  X  10-^ 

yAr 

0.013 

yo 

2.37  X  10-3 

Table  2.  X-2  fill  conditions  for  each  section  of  the  expansion  tunnel. 


Reservoir 

Driver 

Shock  Tube 

Acceleration  &  Test  Section 

6.85  MPa 

742  mBar  He 

3  kPa  Air 

10  Pa  Air 

186  mBar  Ar 

75], 

The  tunnel  condition  used  for  these  tests  was  a  8.6  km/s  flight  speed  equivalent, 
which  is  representative  of  Earth  re-entry  conditions.  The  freestream  conditions  of 
the  test  flow  are  shown  in  Table  1.  An  estimate  of  the  freestream  conditions  was 
detailed  by  Zander  [131],  which  was  a  one-dimensional  nozzle  simulation  that  used  an 
equilibrium  assumption  to  calculate  the  chemical  composition  of  the  flow.  The  X-2 
All  pressures  for  each  section  of  the  expansion  tunnel  are  specified  in  Table  2. 


2.2  Graphite  Test  Model 


The  conducted  tests  used  a  half-cylinder  model  that  had  a  50  mm  outer  diameter, 
a  46  mm  inner  diameter,  and  was  10  mm  thick.  The  material  used  for  all  test  models 
was  GM-10  grade  isotropic  graphite  manufactured  by  Graphel,  LEG.  The  graphite  is 
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Figure  3.  SEM  image  of  the  ISO-63  grade  graphite  surface  at  a  magnification  of  100  /im. 


isostatically  pressed  to  produce  a  very  fine  grain,  high  density  sample.  The  graphite 
has  an  amorphous  structure  with  no  defined  orientation  as  shown  in  Figure  3.  The 
surface  was  imaged  using  a  scanning  electron  microscope  (SEM)  at  a  magnification  of 
100  pm.  The  resistivity  and  density  of  the  graphite  is  1500  pO/cm  and  1.78  g/cm^, 
respectively.  The  amount  of  resistive  heating  is  driven  by  the  resistance  of  the  model, 
which  is  a  function  of  the  density.  A  higher  resistance  translates  to  larger  achievable 
wall  temperatures. 

The  model  is  positioned  in  the  center  of  the  core  how  exiting  the  nozzle,  which 
generated  a  three-dimensional  how  field.  Figure  4  displays  a  schematic  and  actual 
view  of  the  model  located  in  the  tunnel  during  testing.  A  new  model  was  used  for 
each  test  because  the  model  gets  destroyed  after  arrival  of  the  driver  gas. 
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(a)  Schematic 


(b)  Actual 


Figure  4.  Schematic  and  actual  view  of  carbon  model  in  the  X-2  test  section  [133]. 
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2.3  Techniques 


2.3.1  Model  Pre-Heating. 

Pre-heating  of  models  can  be  nsed  to  eliminate  the  time  reqnirements  necessary  to 
correctly  simnlate  the  effects  of  ablation  and  snrface  chemistry  in  expansion  tnnnels. 
The  time  reqnired  to  observe  ablation  effects  in  expansion  tnnnels  is  on  the  order  of 
seconds  for  models  initially  at  room  temperatnre  [74].  However,  expansion  tnnnels 
only  have  test  times  on  the  order  of  microseconds.  Zander  et  ah  [133]  has  shown 
that  pre-heating  of  the  model  can  remove  this  time  reqnirement  by  achieving  re-entry 
snrface  temperatnres  that  allows  for  the  stndy  of  ablation  and  gas-snrface  interactions 
in  hypersonic  flow  fields.  The  pre-heating  is  done  by  sending  an  electric  cnrrent 
throngh  the  model,  which  gets  condncted  throngh  the  entire  cross-section  resnlting 
in  a  nniform  snrface  temperatnre.  Prior  to  each  test,  the  model  was  raised  to  the 
desired  temperatnre  and  held  at  that  temperatnre  for  approximately  5  seconds  before 
firing.  The  model  temperatnre  was  controlled  by  selecting  an  appropriate  amperage 
with  the  power  snpply. 

2.3.2  Two  Color  Ratio  Pyrometry. 

Determination  of  the  model  snrface  temperatnre  was  achieved  by  the  two  color 
ratio  pyrometry  (TCRP)  techniqne.  It  is  necessary  to  know  the  temperatnre  of  the 
model  snrface  before  the  flow  arrives  in  order  to  replicate  these  experiments  with 
nnmerical  simnlations.  The  TCRP  techniqne  is  described  by  Zander  et  ah  [133],  and 
it  can  be  implemented  with  a  commerically  available  digital  single-lens  reflex  (DSLR) 
camera. 

Two  color  ratio  pyrometry  relies  on  the  ratio  of  light  intensities  at  two  different 
known  wavelengths.  Digital  cameras  provide  this  capability  by  captnring  intensities 
for  red,  green,  and  bine  wavelengths.  The  two  color  ratios  can  be  nsed  to  calcnlate  a 
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Figure  5.  Representative  two-color  ratio  pyrometry  analysis  displaying  temperature 
contours  derived  from  a  DSLR  image  of  the  graphite  model.  Image  taken  immediately 
before  test  time. 


temperature  for  each  pixel  using  the  image  data.  The  graphite  model  was  assumed  to 
be  radiating  as  a  grey  body  (i.e.  constant  emissivity),  which  is  consistent  with  other 
research  involving  carbon-carbon  materials  in  near- vacuum  conditions  [133,  10],  The 
camera  used  for  the  tests  was  a  Canon  400D  DSLR  with  a  Canon  EF  75-300  mm 
lens.  The  images  taken  targeted  the  front  surface  of  the  model.  Typical  settings  used 
during  the  tests  were  ISO  100,  1/4000  exposure,  f/10  aperture  and  a  focal  length  of 
300  mm.  An  example  TCRP  analysis  is  shown  in  Figure  5  illustrating  a  fairly  uniform 
temperature  across  the  entire  surface. 

2.3.3  High  Speed  Imaging. 

Recordings  of  each  experimental  test  were  done  with  a  Shimadzu  HPV-1  Hyper¬ 
vision  high  speed  CCD  camera.  The  high  speed  camera  is  capable  of  frame  rates  up 
to  1  MHz  and  records  the  luminosity  of  the  flow  over  the  model.  The  high  speed 
video  from  each  test  was  used  to  check  the  timing  of  the  spectrometer’s  exposure 
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with  respect  to  shock  arrival,  flow  stabilization,  and  the  end  of  steady  test  time.  The 
videos  for  all  the  tests  were  recorded  at  a  frame  rate  of  500  kHz,  which  allowed  for 
all  stages  of  flow  development  to  be  observed  with  a  total  recording  time  of  202  /rs. 

An  example  analysis  is  shown  in  Fignre  6  with  a  model  at  a  snrface  temperatnre  of 
1920  K.  Fignre  6b  exhibits  a  time  where  steady  conditions  were  achieved  and  6c  shows 
when  the  spectrometer  was  first  exposed.  Any  test  that  had  an  excessive  amonnt  of 
Inminosity  dne  to  contaminants,  snch  as  iron  and  alnminnm,  dnring  the  spectrometer 
exposnre  time  were  rejected.  Fignre  6e  indicates  the  end  of  steady  test  time  and 
shows  driver  gas  arrival. 

2.3.4  Ultraviolet  Spectrometry. 

Measnrements  of  shock  layer  emissions  were  performed  nsing  an  Acton  Research 
SpectraPro  2300i  spectrometer  conpled  to  a  Princeton  Instrnments  PI-MAX  ICCD 
UV-sensitive  camera.  The  spectrometer  was  conflgnred  nsing  a  grating  with  a  groove 
density  of  1800  g/mm,  centered  at  372.5  nm  that  covered  a  wavelength  range  of 
353-391  nm.  The  wavelength  range  was  chosen  to  provide  the  maximnm  spectral 
resolntion  of  the  primary  radiative  transitions  for  the  CN  violet  bands.  The  camera 
nsed  a  20  /rs  delay  and  an  exposnre  time  of  15  /is  dnring  steady  test  time. 

Dnring  initial  testing  at  snrface  temperatnres  in  excess  of  3000  K,  attempts  were 
made  to  observe  radiation  from  the  C2  Swan  and  C3  Swings  bands.  For  these  tests, 
a  600  g/mm  grating  centered  at  415  nm  was  nsed,  which  allowed  observation  from 
350-480  nm.  No  discernible  signal  from  these  emission  bands  conld  be  measnred. 
Therefore,  it  was  decided  to  focns  only  on  the  previons  spectrometer  settings  that 
targeted  CN. 

The  high  temperatnre  gas  emits  light  throngh  a  pair  of  UV  grade  fnsed  silica 
windows  monnted  on  both  sides  of  the  X-2  test  section.  The  UV  optical  path  is 
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Figure  6.  Example  analysis  of  high  speed  footage  (times  shown  are  with  respect  to 
start  of  video)  [57]. 
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diagrammed  in  Figure  7.  A  concave  focusing  mirror,  flat  turning  mirror,  and  periscope 
were  used  to  redirect  the  light  onto  the  50  /rm  wide  entrance  slit  of  the  spectrometer. 
The  periscope  changes  the  height  of  the  light  path  from  the  test  section  window 
to  the  height  of  the  spectrometer  entrance  slit,  and  reorients  the  image  by  90°  as 
geometrically  required  to  capture  the  model  surface,  shock  layer,  and  freestream  flow. 
All  optical  components  were  UV  enhanced  magnesium  fluoride/aluminium  coated 
mirrors  to  ensure  good  reflectivity  at  the  wavelengths  of  interest  [29]. 

2.4  Acquisition  of  Spectra 

The  radiation  measured  by  the  UV  spectrometer  during  testing  targeted  the  stag¬ 
nation  streamline  region  of  the  model.  Figure  8  shows  the  capture  area  of  the  UV 
spectrometer  that  was  imaged.  Inside  this  capture  area,  the  spectrometer  camera 
accumulates  an  image  of  the  model  surface,  shock  layer,  and  some  of  the  freestream. 
Note  that  the  vertical  dimension  of  the  capture  area  identified  in  Figure  8  has  been 
magnified  for  clarification. 

The  resultant  image  is  a  two-dimensional  plot  of  pixel  intensity  with  wavelength 
and  spatial  distance  as  the  two  axes.  An  example  raw  data  image  is  shown  in  Figure  9 
with  the  shock  front,  boundary  layer  edge,  model  edge,  and  flow  direction  labeled. 
The  calibration  process  is  described  in  Lewis  et  al.  [57],  which  converts  the  raw 
spectral  data,  measured  in  the  arbitrary  units  of  “counts”,  into  spectral  radiance  with 
units  of  W/(cm^-/im-sr).  Figure  10  shows  the  calibrated  spectra,  which  have  been 
spatially  averaged  within  the  boundary  layer  region,  for  the  surface  temperatures 
ranging  from  1770-2410  K.  The  primary  CN  radiative  transitions  are  clearly  seen  in 
this  figure. 
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Figure  7.  Schematic  of  UV  spectrometer  layout  used  during  X-2  testing  (adapted  from 
Eichmann  [29]). 


20 


1000mm 


Figure  8.  Side  view  of  UV  spectrometer  capture  area  on  the  model  (slit  height  not  to 
scale,  adapted  from  Zander  et  al.  [133]). 
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Figure  9.  Raw  spectral  data  for  2170  K  test  (scale  adjusted  for  visibility,  adapted  from 
Lewis  et  al.  [57]). 
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Figure  10.  Spectra  averaged  within  the  boundary  layer  [57]. 


2.5  Cases 

There  were  two  separate  experimental  test  efforts  or  “campaigns”  that  focused 
on  different  ablation  regimes.  The  hrst  campaign  focused  on  the  oxidation  and  ni- 
tridation  reactions  with  a  surface  temperature  range  of  1770-2410  K.  The  second 
campaign  attempted  to  achieve  sublimation  surface  temperatures  with  a  range  of 
2610-3280  K.  Between  the  two  campaigns,  there  were  a  total  of  eight  experimental 
tests  that  produced  “clean”  spectral  results.  The  measured  wall  temperature  along 
with  its  associated  uncertainty  are  provided  in  Table  3.  It  should  be  noted  that  the 
DSLR  images  taken  for  case  3  were  overexposed  and  therefore  unusable.  Therefore, 
the  temperature  supplied  for  case  3  was  taken  from  another  test  with  identical  power 
settings  which  was  rejected  due  to  flow  contamination  [57].  Cases  1-4  were  from  the 
hrst  experimental  campaign  and  cases  5-8  were  from  the  second. 

The  maximum  temperature  for  the  hrst  campaign  was  limited  by  the  available 
power  supply.  The  power  supply  used  for  the  hrst  campaign  was  a  10.75  V  DC 
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Table  3.  Measured  wall  temperature  for  each  test  case. 


Case  Number 

Wall  Temperature  (K) 

Uncertainty  (K) 

1 

2410 

±280 

2 

2170 

±180 

3 

1920 

±180 

4 

1770 

±180 

5 

3280 

±20 

6 

3190 

±40 

7 

2760 

±180 

8 

2610 

±140 

rectifier,  which  could  provide  up  to  250  A  with  this  test  model.  To  achieve  higher 
temperatures,  a  low-ripple  DC  rectifier  capable  of  supplying  up  to  1500  A  was  pur¬ 
chased.  Analysis  by  Lewis  et  al.  [56]  found  that  supplying  approximately  450-500  A 
could  achieve  surface  temperatures  of  3000-3300  K.  A  comparison  of  SEM  images  of 
before  and  after  heating  to  3300  K  is  provided  in  Figure  11.  Heating  the  graphite 
model  to  3300  K  resulted  in  a  significant  change  to  the  surface  structure,  particularly 
an  increase  in  porosity,  which  had  unforeseen  consequences  during  testing  that  will 
be  discussed  later. 

The  presence  of  iron  contamination  in  the  spectrum,  as  seen  in  Figure  10,  was 
problematic  because  it  caused  integrated  radiance  levels  to  be  higher  than  that  of  a 
clean  how.  However,  Lewis  et  al.  [57]  showed  that  the  error  due  to  iron  contamination 
was  not  significant.  Emission  from  CN  was  estimated  to  contribute  75-85%  of  the  total 
measured  radiance  within  the  boundary  layer  region.  The  cases  identified  in  Table  3 
were  selected  for  analysis  because  they  had  the  lowest  levels  of  contamination.  The 
next  chapter  will  discuss  the  numerical  theory  used  to  simulate  these  cases. 
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(a)  Before  heating. 


(b)  After  heating  to  3300  K. 


Figure  11.  SEM  images  taken  at  a  magnification  of  20  fim  comparing  the  model  surface 
before  and  after  heating  to  3300  K. 


24 


III.  Theory  and  Methodology 


Hypersonic  re-entry  flows  inclnde  many  physical  phenomena  that  cannot  be  mod¬ 
eled  by  the  perfect  gas  form  of  the  Navier-Stokes  eqnations.  The  types  of  physical 
phenomena  expected  in  these  flows  inclnde  chemical  and  thermal  noneqnilibrinm, 
vibrational  and  electronic  excitation,  and  weak  ionization.  The  nse  of  a  real  gas 
model,  where  molecnlar  interactions  are  considered,  is  necessary  to  acconnt  for  these 
phenomena. 

The  following  sections  will  describe  the  CFD  and  radiation  solvers  nsed  in  this 
dissertation.  The  CFD  solver  is  called  US3D,  which  is  a  parallel  implicit  nnstrnctnred 
solver  that  was  developed  at  the  University  of  Minnesota  [80,  81].  The  radiation 
solver  is  the  Noneqnilibrinm  Air  Radiation  (NEQAIR)  code  version  13.2  developed  at 
NASA  Ames  Research  Center  [125].  The  CFD  and  radiation  solntions  were  compnted 
in  an  nnconpled  manner  becanse  the  incoming  Mach  nnmber  was  aronnd  9.4  and 
radiative  heating  was  not  expected  to  contribnte  mnch  to  the  overall  wall  heating  rate. 
Additionally,  a  description  of  the  snrface  reaction  model  and  how  it  is  implemented 
as  a  CFD  bonndary  condition  is  provided. 

3.1  Assumptions 

The  Knndsen  nnmber  is  the  ratio  of  the  mean  free  path  to  a  characteristic  flow 
length,  snch  as  the  diameter  of  the  graphite  test  model,  and  can  be  nsed  to  distingnish 
between  continnnm  or  non-continnnm  flow  conditions.  The  continnnm  formnlation 
reqnires  that  the  Knndsen  nnmber  be  mnch  less  than  one  (nsnally  less  than  0.1)  so 
there  are  a  snfflcient  nnmber  of  molecnles  within  a  compntational  volnme.  Hence, 
there  is  little  statistical  variation  at  any  point,  and  as  a  resnlt,  the  continnnm  formnla¬ 
tion  of  the  viscons  flnxes  is  consistent  [22].  The  Knndsen  nnmber  near  the  stagnation 
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region  for  the  X-2  flow  conditions  ranged  from  10“'^  to  10“^.  Therefore,  the  flow  fields 
are  assnmed  to  be  consistent  with  a  continnnm  formnlation.  It  is  also  assnmed  that 
there  are  snfflcient  nnmber  of  collisions  of  the  gas  molecnles  with  the  wall,  so  there 
is  no  velocity  or  temperatnre  slip.  For  continnnm  calcnlations,  the  Navier-Stokes 
eqnations  may  be  solved  to  determine  flow  field  solntions. 

The  thermal  state  of  the  gas  was  assnmed  to  be  in  noneqnilibrinm  and  follows 
Park’s  two-temperatnre  model  [92,  91],  The  energy  in  the  translational  and  rotational 
modes  of  all  gas  species  are  characterized  by  a  single  temperatnre,  T.  A  common 
translational-rotational  temperatnre  is  a  reasonable  assnmption  becanse  typically  ro¬ 
tational  eqnilibration  with  translation  is  very  fast  and  occnrs  within  a  few  collisions. 
The  energy  of  the  vibrational  and  electronic  modes  of  all  species  and  the  electron 
translational  energy  mode  is  described  by  a  single  temperatnre,  T^.  The  jnstiflcation 
that  these  three  energy  modes  can  be  described  by  a  single  vibrational  temperatnre 
is  based  on  the  rapid  energy  transfer  between  the  translational  motion  of  the  free 
electrons  and  the  vibrational  motion  of  the  polyatomic  species,  and  on  popnlations 
of  the  low-lying  electronic  states  of  heavy  particles  tending  to  be  in  eqnilibrinm  with 
the  gronnd  electronic  state  at  the  electronic  temperatnre  [91].  These  assnmptions 
regarding  energy  transfer  simplify  the  conservation  eqnations  considerably  by  elimi¬ 
nating  the  need  for  separate  translational  and  vibrational  energy  eqnations  for  each 
polyatomic  species  and  an  energy  eqnation  for  electrons.  While  these  simplifications 
may  not  be  adeqnate  for  all  high  speed  problems,  they  provide  accnrate  resnlts  for 
aerodynamic  coefficients  and  convective  heating  rates  of  re-entry  vehicles  [31]. 
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3.2  CFD  Solver 


3.2.1  Conservation  Equations. 

The  gas  dynamic  conservation  eqnations  for  a  individnal  gas  species,  s,  in  a 
noneqnilibrinm  flow  are  presented  in  this  section.  The  eqnations  that  describe  the 
conservation  of  mass,  momentnm,  and  energy  for  a  gas  species  follow  the  work  in 
References  [7,  120,  54],  The  species  mass  conservation  eqnation  is  given  by: 


(1) 


where  ps  are  the  species  densities,  uj  is  the  velocity  in  the  j  direction,  Vsj  is  the 
species  diffnsion  velocity,  and  Ws  are  the  species  mass  prodnction  rates  dne  to  chemical 
reactions.  The  momentnm  conservation  eqnations  are  written  as: 


(2) 


where  p  is  the  mixtnre  density,  p  is  the  pressnre,  and  Tij  are  the  viscons  shear  stress 
components.  There  is  a  separate  momentnm  eqnation  for  each  spatial  dimension  of 
the  problem.  The  total  energy  conservation  eqnation  is  represented  by: 


dE  d  .  .  d  .  .  d  . 


where  E  is  the  total  energy,  qj  is  the  translational-rotational  heat  flnx  in  the  j  di¬ 
rection,  is  the  vibrational-electronic-electron  heat  flnx  in  the  j  direction,  Hg  is 
the  species  enthalpy,  and  ns  is  the  total  nnmber  of  species.  Finally,  the  vibrational- 
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electronic-electron  conservation  eqnation  is  given  by: 


dE, 

dt 


dx  {E,u,)  + 
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(4) 


where  E^j  is  the  vibrational-electronic  energy,  is  the  vibrational-electronic  energy 
of  a  polyatomic  species,  Wy  is  the  vibrational  energy  sonrce  term,  and  nd  is  the  total 
nnmber  of  polyatomic  species.  The  following  sections  provide  more  details  on  how 
these  terms  are  modeled. 


3.2.2  Viscous  Terms. 


The  viscons  shear  stresses  are  modeled  assnming  a  Newtonian  flnid  and  are  given 


by: 


/  dui  du. 

Tij  =  -/i  7^  + 


OXk 


(5) 


\dxj  '  dxi  j  uxk 
where  fi  is  the  mixtnre  coefficient  of  viscosity  and  A  =  — 2/3/r  from  Stoke’s  hypothesis. 

The  translational-rotational  and  vibrational-electronic  heat  flnx  vectors  are  given 
by  Fonrier’s  heat  law: 


qj  =  -k 


dT 

dx/ 


Q.v,j 


dTy 
’  dxi 


(6) 


where  k  is  the  mixtnre  translational-rotational  thermal  condnctivity  and  ky  is  the  mix¬ 
tnre  vibrational-electronic  thermal  condnctivity.  The  species  mass  diffnsion  flnxes, 
Js,  are  modeled  nsing  Pick’s  law: 


Js  ps'^sj 


(7) 


where  Dg  is  the  species  diffnsion  coefficient  and  ys  is  the  species  mass  fraction.  The 
species  diffnsion  coefficient,  Dg,  is  replaced  by  a  single  binary  coefficient,  D,  and 
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found  by  assuming  a  constant  Lewis  nnmber,  Le,  by: 


D 


Lek 

pCp 


(8) 


where  k  is  the  mixtnre  translational-rotational  thermal  condnctivity  and  Cp  is  the 
mixtnre  translational-rotational  specihc  heat  at  constant  pressnre.  The  Lewis  nnmber 
is  set  to  1.4. 

There  have  been  a  few  recent  studies  investigating  the  effects  of  different  mass 
diffusion  models  for  planetary  re-entry  vehicles  inclnding  wall  catalytic  effects.  Gosse 
and  Candler  [35]  showed  surface  heating  rates  and  mass  fraction  compositions  com- 
pnted  by  Tick’s  law  began  to  diverge  from  mnlti-component  diffnsion  models  when 
freestream  velocities  were  larger  than  9.7  km/s.  The  main  reasons  for  the  large  dif¬ 
ferences  at  freestream  velocities  greater  than  9.7  km/s  were  due  to  the  flow  becoming 
increasingly  dissociated  and  ionized.  Specihcally,  the  increased  presence  of  N  and  N"*" 
were  cansing  the  snrface  heating  predictions  given  by  Tick’s  law  to  be  larger  than 
those  predicted  by  the  mnlti-component  methods.  Alkandry  et  al.  [5]  fonnd  that 
Tick’s  law  consistently  predicted  larger  stagnation  point  heat  transfer  rates  compared 
to  a  modihed  Tick’s  law  and  other  mnlti-component  diffnsion  models.  The  sonrce  of 
the  problem  was  due  to  Tick’s  law  not  enforcing  the  reqnirement  that  the  diffnsion 
mass  fluxes  sum  to  zero.  Both  studies  used  simple  catalytic  wall  models  and  how  the 
resnlts  wonld  change  with  a  more  sophisticated  snrface  reaction  model  is  not  clear. 
However,  the  freestream  velocity  of  the  experimental  test  condition  is  8.5  km/s,  so 
nsing  Tick’s  law  is  assnmed  to  be  reasonable  for  the  cases  presented  here. 
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3.2.3  Thermodynamic  Properties. 


The  total  pressure  is  found  from  the  equation  of  state  and  is  the  sum  of  the  partial 
pressures: 


p  =  Yl 


JLt 


(9) 


where  R  is  the  universal  gas  constant  and  Ms  is  the  species  molecular  weight.  The 
total  energy  of  the  mixture  is  given  by: 


no  ^  /fccj  na 

E  =  psUiUi  +  pshi + e^  (lo) 

S=1  S=1  S  =  1 

where  h°  is  the  species  heat  of  formation.  This  expression  may  be  inverted  to  deter¬ 
mine  the  translational-rotational  temperature  T.  The  translational-rotational  specihc 
heat  at  constant  volume,  C^^s,  is  given  by: 


Cv^s  —  Cl 


v,tr,s 


+  C 


v.rot.s 


(11) 


where  Cy^tr,s  —  \  C.u^rot,s  =  for  molecules,  and  Cy^rot,s  =  0  for  monatomics.  The 
species  enthalpy  is  dehned  to  be: 


=  Cy^sT 


RgT  -\-  h° 


V^S 


(12) 


where  Rg  =  R/Mg. 

The  species  vibrational-electronic  energy  is  calculated  using  the  NASA  Chemi¬ 
cal  Equilibrium  with  Applications  (CEA)  database.  Gordon  and  McBride  [34]  and 
McBride  et  al.  [70]  developed  a  database  that  can  be  used  to  calculate  thermody¬ 
namic  functions  for  individual  gas  species.  Thermodynamic  data  can  be  calculated 
for  a  temperature  range  from  200  to  20,000  K  with  a  library  of  over  2,000  solid,  liquid, 
and  gaseous  chemical  species.  The  data  are  expressed  as  least-squares  coefhcients  for 
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the  calculation  of  the  specihc  heat  at  constant  pressure  Cp,  enthalpy  H,  and  entropy 
S.  From  these  thermodynamic  quantities,  the  species  vibrational-electronic  energy 
can  be  calculated  from  the  species  enthalpy  (Equation  12),  and  the  species  vibrational 
specihc  heat  at  constant  volume  can  be  calculated  from  the  specihc  heat  at  constant 
pressure  since  Cy  =  Cp  —  R. 

The  thermodynamic  functions  in  dimensionless  form  are  as  follows.  The  specihc 
heat  at  constant  pressure  is: 

+  ^  +  as  +  a^T  +  agT^  +  (13) 

Enthalpy  and  entropy  are  obtained  by  integrating  Cp(T)  and  Cp(T)/T  with  respect 
to  T: 


H{T) 

RT 


Oi  02  hi(T) 


bi 


- '  +  Os  +  2^4^  +  + -oyT  +  — 


(14) 


and 


^  ^  —  ^  +  ^3  In(^)  +  (^4:T  + 

where  oi  —  07  are  the  least-square  coefficients  and  6iand  62  are  integration  constants. 
The  thermodynamic  database  inherently  accounts  for  both  the  vibrational  and  elec¬ 
tronic  internal  energy  modes,  making  it  more  accurate  than  a  harmonic  oscillator 
model  for  high  energy  hows. 


3.2.4  Transport  Properties. 

The  viscosity  and  thermal  conductivities  for  each  energy  mode  are  calculated 
according  to  the  Gupta- Yos  mixing  rule,  which  was  designed  for  weakly  ionized  or 
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non-ionized  flows  [40],  The  formulas  used  to  compute  the  trausport  properties  from 
the  collisiou  cross-sectious  are  obtaiued  from  the  first-order  Chapmau-Euskog  ap- 
proximatiou.  For  a  gas  iu  thermal  equilibrium,  the  mixture  viscosity  is  calculated 
by: 


2  =  1 


M, ; 

N, 


-Xi 


\ 


V 


(16) 


where  the  collisiou  terms,  A^,  are  evaluated  at  the  coutrolliug  temperature,  T,  aud 
the  species  molar  fractiou,  Xs,  is  giveu  by: 


Ps/Ms 
EZi  Ps/Ms 


(17) 


The  trauslatioual  compoueut  of  the  thermal  couductivity  iu  a  mixture  is: 


ktr  —  .  Pb  ^  ^ 


Xi 


(18) 


where  is  the  Boltzmann  constant  and  aij  is  dehned  as: 


C^ij  —  1  “h 


|1  -  (Afj/Mj)]  [0.45  -  2.54  (Mj/Af,)] 
|1  +  {M,/M,p 


(19) 


The  expressiou  for  the  rotatioual  thermal  couductivity  is: 


Xi 


Kot  =  kb  7  —  1 

i=raol  A/ 7=1  Xj^ij 


(20) 


aud  the  vibratioual-electrouic  thermal  couductivity  is: 


kvib  k}y 


a 


v.vib 


Xi 


Rs  ^ 

i=mol  Z—/j=l  Xj^ij 


(21) 
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Note  that  in  the  above  expression  is  evalnated  nsing  the  CEA  database.  The 

mixtnre  translational-rotational  thermal  condnctivity  is  then  given  by  the  snm  of  each 
energy  mode  as: 


(22) 


The  collision  terms  are  dehned  as: 


(23) 


and 


(24) 


where  the  collision  integrals  are  the  weighted  average  cross-section  of  a  collision 


between  species  i  and  j  [40],  The  constant  10“^°  converts  to  sqnare  meters  from 
sqnare  Angstroms.  The  evalnation  of  the  collision  integrals  can  be  obtained  throngh 
a  variety  of  methods.  The  recommended  evalnation  methods  are  given  in  Reference 
[126],  which  provides  a  complete  set  of  collision  integrals  for  the  compntation  of 
transport  properties. 

3.2.5  Source  Terms. 

The  sonrce  terms  appear  in  the  species  mass  (Eqnation  1)  and  vibrational-electronic- 
electron  energy  (Eqnation  4)  conservation  eqnations.  The  formnlation  of  the  species 
mass  prodnction  rates  dne  to  chemical  reactions,  Wg,  and  the  vibrational  energy  sonrce 
term,  Wy,  are  discnssed  here. 

Chemical  Reactions. 

The  species  mass  prodnction  rate  sonrce  terms  are  derived  from  the  chemical 
reactions  that  occnr  in  the  gas.  Dne  to  the  high  freestream  enthalpies  and  velocities 
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of  these  flows,  some  degree  of  chemical  nonequilibrium  was  expected  in  the  post-shock 
conditions.  Hence,  a  20-species,  40  reaction  finite-rate  chemistry  model  for  reacting 
air  that  includes  carbon  species  is  used  in  this  work.  The  species  considered  in  the 
model  are  Na,  O2,  NO,  CO2,  CO,  C2,  C3,  ON,  NO+,  N^,  O^,  CO+,  Ar,  C,  N,  O, 
N"*",  O’*",  C^,  and  e~ .  The  chemical  production  rate  of  species  s  in  reaction  q  can  be 
represented  generally  as  [31,  107]: 


^s,q 


(26) 


where  a  and  (3  are  the  stoichiometric  coefficients  for  each  reaction.  The  source  terms 
in  the  species  mass  conservation  equation  are  given  by: 


Ws  =  Ms'^Ws,q  (26) 

<? 

Each  reaction  is  governed  by  forward  and  backward  reaction  rate  coefficients,  kf^q 
and  fcft  q,  respectively.  The  forward  reaction  rates  are  calculated  using  an  Arrhenius 
equation: 

kf,q  =  exp  i-Df^q/Tf^q)  (27) 

where  the  constants  Af^q,  nf^q  and  Df^q  are  experimentally  determined.  Most  of  the 
chemical  reactions  used  Park  rates  with  reaction  coefficients  determined  in  air  where 
applicable.  However,  a  majority  of  the  rate  coefficients  for  the  reactions  involving 
carbon  species  were  determined  assuming  a  Martian  atmosphere.  The  40  reactions 
considered  are  listed  in  Table  4.  The  reactions  involving  C3  were  not  considered  in 
any  of  the  references  attributable  to  Park,  but  under  strong  sublimation  conditions 
these  reactions  become  important.  Following  the  work  of  Candler  [23],  the  dissocia¬ 
tion,  oxidation,  and  nitridation  reactions  for  C3  were  included.  Candler  showed  the 
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formation  of  C3  does  not  become  appreciable  until  the  surface  temperature  reaches 
approximately  3000  K  for  an  equilibrium  air-carbon  mixture. 

The  backward  reaction  rates  are  obtained  from: 


^b,q 


kf,q 

K 


(28) 


where  Keq^q  is  the  concentration-based  equilibrium  constant  for  reaction  q.  The 
concentration-based  equilibrium  constant  is  calculated  using  the  activity-based  equi¬ 
librium  constant  Ka^q'. 


K  =  K 

^^eq,q 


a,q 


\RTJ 


(29) 


where  the  net  stoichiometric  exponent  Ugq  considers  only  the  gas  species  in  reaction 
q: 

Kg 

iy'kq  -  Y)  (30) 

k=l 

The  activity-based  equilibrium  constant  is  computed  using  Gibbs  free  energy  mini¬ 
mization  and  the  NASA  CEA  curve  hts  for  enthalpy  and  entropy: 


Ka,q  =  exp 


-Agg(T) 

RT 


=  exp 


K 

^kq 

k=l 


Skim 

R  ) 


(31) 


where  G°(T)  is  Gibbs  free  energy  and  Vkq  is  the  net  stoichiometric  for  species  k  in 
reaction  q.  The  activity-based  equilibrium  constant  can  be  calculated  directly  if  the 
necessary  thermodynamic  functions  are  available  for  each  species  in  the  reaction.  The 
US3D  code  uses  the  CEA  thermodynamic  database  with  a  reference  pressure  of  10^  Pa 
to  compute  the  activity  and  concentration-based  equilibrium  constants. 

Eor  dissociation  reactions,  the  forward  rates  should  be  governed  by  the  translational- 
rotational  and  vibrational-electronic  temperatures.  The  exact  form  of  the  tempera¬ 
ture  dependencies  on  the  dissociation  rate  is  unknown,  but  a  few  theories  have  been 
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postulated  in  the  literature.  Hammerling  et  al.  [41]  developed  the  coupled  vibra¬ 
tion  dissociation  (CVD)  model  where  dissociation  rate  expressions  were  derived  for  a 
rotationless  harmonic  oscillator  with  vibrational  levels  populated  according  to  a  Boltz¬ 
mann  distribution.  The  CVD  model  provides  a  dissociation  rate  that  is  a  function  of  T 
and  Ty,  but  proved  to  be  inadequate  because  the  vibrational  energy  removed  by  disso¬ 
ciation  was  neglected.  Park  [92]  suggested  that  the  dissociation  rate  is  governed  by  a 
geometric  average  between  the  translational-rotational  and  the  vibrational-electronic 
temperatures: 


Ta  =  (32) 

This  average  temperature  accounts  for  vibration-dissociation  coupling  [91].  A  more 
general  expression  is  also  typically  used: 

Ta  =  (33) 

where  the  parameter  q  is  taken  to  be  somewhere  between  0.3  and  0.5.  Park  notes  that 
the  CVD  rate  can  be  approximately  reproduced  when  q  =  0.3,  but  this  value  may 
underestimate  the  influence  of  the  vibrational  temperature  [91].  Previous  work  has 
also  shown  that  radiative  heating  calculations  are  extremely  sensitive  to  the  choice 
of  q.  Hartung  [42]  showed  that  a  g  of  0.5  produced  a  peak  radiative  heating  value 
that  was  double  that  of  choosing  q  equal  to  0.3.  The  US3D  code  employs  Equa¬ 
tion  32  to  govern  the  dissociation  rate  and  is  used  for  the  simulations  presented  in 
this  dissertation.  Yet,  the  choice  of  q  is  acknowledged  as  an  uncertain  parameter. 

Additionally,  a  kinetics  model  proposed  by  Johnston  et  al.  [46]  that  mainly  mod- 
ihes  the  Park  CO  and  CO2  dissociation  rates  was  also  used.  Particular  to  this  study, 
the  CO  dissociation  rate  is  increased  by  a  factor  of  13  compared  to  the  Park  rate.  For 
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Table  4.  Reactions  and  rate  coefficients  for  the  gas-phase  chemistry  model. 


q 

Reaction 

^f,q 

Df,q 

Tf,q 

Third  Body,  M 

Ref. 

1 

CO2  +  M  ^  CO  +  0  +  M 

1.4e+19 

-1.50 

6.328e+4 

Ta 

C,  N,  0,  N+,  0+,  C+ 

[96] 

6.9e+17 

-1.50 

6.328e+4 

T 

Ar 

[96] 

6.9e+18 

-1.50 

6.328e+4 

T 

others 

[96] 

2 

CO  +  Me^C  +  O  +  M 

3.4e+17 

-1.00 

1.29e-75 

T 

C,  N,  0,  N+,  0+,  C+ 

[96] 

2.3e+16 

-1.00 

1.29e-55 

T 

Ar 

[96] 

2.3e+17 

-1.00 

1.29e-75 

T 

others 

[96] 

3 

Na  +  M^N  +  N  +  M 

3.0e+19 

-1.60 

1.132e-75 

T 

C,  N,  0,  N+,  0+,  CH- 

[93] 

3.0e+21 

-1.60 

1.132e-75 

T 

e~ 

[93] 

7.0e+18 

-1.60 

1.132e-h5 

T 

others 

[93] 

4 

02  +  M«->0  +  0  +  M 

l.Oe+19 

-1.50 

5.95e+4 

T 

C,  N,  0,  N+,  0+,  C+ 

[93] 

2.0e+18 

-1.50 

5.95e+4 

T 

others 

[93] 

5 

NO  +  Me^N  +  O  +  M 

l.le+14 

0.00 

7.55e+4 

T 

C,  N,  0,  N+,  0+,  C+,  NO,  CO2 

[93] 

5.0e+12 

0.00 

7.55e+4 

T 

others 

[93] 

6 

C2  +  M^C  +  C  +  M 

3.7e+ll 

0.00 

6.99e+4 

T 

All 

[96] 

7 

CN  +  M  ^  C  +  N  +  M 

2.5e+ll 

0.00 

8.774e+4 

T 

All 

[97] 

8 

C3  +  M  C2  +  C  +  M 

3.7e+ll 

0.00 

6.99e-|-4 

T 

All 

[23] 

9 

N  +  <->  N+  + 

2.5e+31 

-3.82 

1.682e-55 

T 

[97] 

10 

0  +  <->  0^  + 

3.9e+30 

-3.78 

1.585e-75 

T 

V 

[93] 

11 

C  +  <->  C+  + 

3.7e+28 

-3.00 

1.307e-75 

T 

[97] 

12 

NO  +  0  ^  O2  +  N 

8.4e+9 

0.00 

1.94e+4 

T 

[93] 

13 

N2  +  0  ^  NO  +  N 

6.4e+14 

-1.00 

3.84e+4 

T 

[93] 

14 

CO  +  0  O2  +  C 

3.9e+10 

-0.18 

6.92e+4 

T 

[96] 

15 

CO  +  C  C2  +  0 

2.0e+14 

-1.00 

5.80e+4 

T 

[96] 

16 

CO  +  N  CN  +  0 

l.Oe+11 

0.00 

3.86e+4 

T 

[96] 

17 

N2  +  C  ^  CN  +  N 

l.le+11 

-0.11 

2.32e+4 

T 

[96] 

18 

CN  +  0  ^  NO  +  C 

1.6e+10 

0.10 

1.46e+4 

T 

[96] 

19 

CN  +  C  C2  +  N 

5.0e+10 

0.00 

1.30e+4 

T 

[96] 

20 

CO2  0  O2  “h  CO 

2.1e+10 

0.00 

2.78e+4 

T 

[96] 

21 

C2  +  C2  C3  +  C 

3.0e+19 

0.00 

0.00 

T 

[23] 

22 

0  +  C3  ^  CO  +  C2 

3.0e+10 

0.00 

0.00 

T 

[23] 

23 

N  +  C3  ^  C2  +  CN 

3.0e+9 

0.00 

0.00 

T 

[23] 

24 

N  +  0  ^  NO+  +  e- 

5.3e+9 

0.00 

3.19e+4 

T 

[92] 

25 

N  +  N  Nj  +  e- 

2.0e+10 

0.00 

6.75e+4 

T 

[92] 

26 

0  +  0  <->  0  J 

l.le+10 

0.00 

8.06e+4 

T 

[92] 

27 

C  +  0^  CO+  +  e- 

8.8e+5 

1.00 

3.31e+4 

T 

[96] 

28 

0+  +  Na  ^  N+  +  0 

9.0e+8 

0.36 

2.28e+4 

T 

[92] 

29 

0+  +  NO  ^  N+  +  O2 

1.4e+2 

1.90 

1.53e+4 

T 

[92] 

30 

NO+  +  O2  0+  +  NO 

2.4e+10 

0.41 

3.26e+4 

T 

[92] 

31 

NO+  +  N  Nj +  0 

7.2e+10 

0.00 

3.55e+4 

T 

[92] 

32 

NO+  +  0  N+  +  O2 

l.Oe+9 

0.50 

7.72e+4 

T 

[92] 

33 

0+  +  N  N+  +  O2 

8.7e+10 

0.14 

2.86e+4 

T 

[92] 

34 

OJ  +  N2  ^  N+  +  O2 

9.9e+9 

0.00 

4.07e+4 

T 

[92] 

35 

Oj  +  0  ^  0+  +  O2 

4.0e+9 

-0.09 

1.80e+4 

T 

[92] 

36 

NO+  +  N  0+  +  Na 

3.4e+10 

-1.08 

1.28e+4 

T 

[92] 

37 

NO+  +  0  OJ  +  N 

7.2e+9 

0.29 

4.86e+4 

T 

[92] 

38 

O2  +  C+  ^  0+  +  C 

l.Oe+10 

0.00 

9.40e+3 

T 

[96] 

39 

CO  +  C+  ^  CO+  +  C 

l.Oe+10 

0.00 

3.14e+4 

T 

[96] 

40 

NO+  +  C  NO  +  C+ 

l.Oe+10 

0.00 

2.32e+4 

T 

[96] 
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Table  5.  Park  gas-phase  chemistry  rate  modifications  made  by  Johnston  et  al.  [46] 


Reaction 

Mnltiplication  Eactor 

CO2  +  M^C0  +  0  +  M 

2 

CO  +  M^C  +  O  +  M 

13 

CN  +  0  ^  NO  +  C 

0.1 

CN  +  C  ^  C2  +  N 

3.2 

CO2  +  0  ^  O2  +  CO 

6 

this  model,  Johnston  et  al.  [46]  made  modihcations  based  on  comparisons  with  shock 
tnbe  radiation  measnrements  nnder  Mars  re-entry  conditions  in  the  NASA  EAST  fa¬ 
cility.  Table  5  lists  a  snmmary  of  the  reaction  rate  modihcations.  Brandis  et  al.  [20] 
fonnd  this  kinetics  model  improved  noneqnilibrinm  radiance  predictions  compared  to 
the  standard  Park  gas-phase  chemistry  model  with  experiments  condncted  in  EAST. 

Vibrational  Energy. 

The  vibrational  energy  sonrce  term  is  given  by: 

Qc—v  T  ^  ^  Qt—v,s 

s 

where  Qc-v  represents  the  amonnt  of  vibrational-electronic  energy  added  or  removed 
dne  to  chemical  reactions.  The  Qc-v  sonrce  term  can  be  dehned  with  a  preferential 
or  non-preferential  model  for  dissociation.  The  non-preferential  model  assnmes  that 
molecnles  are  created  or  destroyed  at  the  average  vibrational  energy.  The  preferential 
model  assnmes  that  molecnles  are  more  likely  to  dissociate  at  higher  vibrational 
energy  levels,  which  tends  to  snppress  the  vibrational  energy  and  dissociation  rate. 
The  amonnt  of  energy  added  or  removed  is  typically  set  to  an  arbitrary  fraction 
(nsnally  0.3)  of  the  dissociation  energy.  However,  the  energy  removal  rate  mnst  be 
consistent  with  the  details  of  the  vibration-dissociation  conpling  model  [84].  The 
two-temperatnre  model  nsed  in  US3D  is  not  based  on  a  detailed  representation  of 
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the  dissociation  process,  so  it  is  impossible  to  derive  an  expression  for  the  vibrational 
energy  added  or  removed  by  dissociation  [23],  Therefore,  the  non-preferential  model 
for  Qc-v  is  nsed  and  dehned  as: 


Q c—v  ^  ^  s^v^s 

s 

Preferential  models  have  been  developed  and  stndied  in  the  past  where  the  energy 
contained  in  the  npper  vibrational  levels  was  acconnted  for  in  the  dissociation  process 
[117,  65,  109],  However,  these  models  are  compntationally  intensive  and  not  often 
implemented  in  CFD  codes. 

The  Qt-v,s  term  is  the  energy  exchange  rate  between  the  vibrational-electronic 
and  translational-rotational  energy  modes.  It  is  assnmed  to  be  modeled  by  a  single 
exchange  rate  dne  to  the  fast  energy  transfer  between  the  translation-rotational  and 
vibrational-electronic  modes  [92].  The  rate  of  change  in  vibrational  state  popnlation 
levels  can  be  described  by  the  Landan-Teller  formnlation  where  it  is  assnmed  that 
molecnles  behave  as  harmonic  oscillators  and  the  vibrational  level  can  only  change 
one  qnantnm  level  at  a  time  [120].  The  single  energy  exchange  rate  for  each  species 
is: 

n  —  n  ~  (ofi\ 

^T—v,s  Ps 

Ts 

where  e„  ^  is  the  vibrational  energy  evalnated  either  at  the  local  translational-rotational 
temperatnre  or  at  the  local  vibrational-electronic  temperatnre.  Again,  the  internal 
vibrational  energy  is  compnted  nsing  the  CEA  database.  The  relaxation  time,  r^,  is 
dehned  as: 

rs=<  Ts,lt  >  +rcs  (37) 

where  <  >  is  the  molar  averaged  Landan-Teller  relaxation  time  [54]  and  t^s  is 


39 


the  collision  limited  relaxation  time.  The  molar  averaged  relaxation  time  is  given  as: 


<  Ts,LT  >  — 


'Yhr  Tr 


X.r/'Ts'. 


(38) 


where  is  the  Landan-Teller  inter-species  relaxation  time  and  is  modeled  nsing 
cnrve  hts  developed  by  Millikan  and  White  [72]  as: 


Tsr  =  -exp[y4sr(T  —  O.OlS/i^/"^)  —  18.42]  ,  p  in  atm, 

A,r  =  1.16  X  (39) 

llsr  =  MMr!  {M,  +  Mr) 


The  relaxation  time  in  Eqnation  36  wonld  calcnlate  a  relaxation  rate  nnrealistically 
large  at  high  temperatnres  dne  to  an  over  prediction  of  the  collision  cross-section  if 
only  the  Landan-Teller  rate  expression  from  Millikan  and  White  was  nsed  [22].  Park 
[92]  corrected  this  problem  by  introdncing  a  collision  limited  relaxation  time,  Tcs, 
which  is  written  as: 

1 

“  arCsN 

where  N  is  the  nnmber  density  of  the  mixtnre,  Cg  is  the  average  molecnlar  speed  of 
species  s  given  as: 

8RT 

ttMs 

and  is  the  limiting  cross  section  given  by: 


/50,000 

V  T 


m 


(42) 


where  ayg  is  the  limiting  cross  section  for  species  s.  This  expression  was  originally 
developed  for  nitrogen,  so  most  molecnles  defanlt  to  the  limiting  cross  section  of 
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Table  6.  Limiting  cross  sections  for  different  gas  species  [93,  96]. 


Molecular  Species  ays  (m^) 


N2 

3 

X 

10- 

-21 

NO 

3 

X 

10- 

-21 

O2 

3 

X 

10- 

-21 

CO 

3 

X 

10- 

-22 

CO2 

1 

X 

10- 

-20 

nitrogen  if  not  specified.  Table  6  provides  the  limiting  cross  sections  used  in  US3D. 

3.2.6  Numerical  Methods. 

The  numerical  methods  utilized  in  US3D  are  similar  to  those  used  in  the  NASA 
Ames  DPLR  Navier-Stokes  solver  [127,  128]  because  US3D  was  originally  intended 
to  be  a  follow-on  to  that  code.  US3D  shares  many  of  the  same  numerics  and  has  also 
been  extensively  validated  against  the  NASA  DPLR  code  on  a  wide  variety  of  test 
problems  [43].  The  compressible  Navier-Stokes  equations  that  also  account  for  inter¬ 
nal  energy  relaxation  and  hnite-rate  chemical  kinetics  are  solved  in  US3D  using  the 
hnite-volume  formulation.  Convective  fluxes  are  calculated  using  the  modihed  Steger- 
Warming  flux  vector  splitting  method  for  steady-state  simulations.  For  unsteady  sim¬ 
ulations,  low-dissipation  centrally  differenced  fluxes  are  available  with  second,  fourth, 
and  sixth  order  spatial  accuracy.  There  are  hrst,  second,  and  third  order  accurate 
explicit  and  implicit  time  integration  methods.  The  MUSCL  scheme  with  a  pressure 
limiter  is  used  to  limit  overshoots  and  prevent  non-physical  values  in  regions  with 
strong  shocks.  Cell-centered  gradients  are  calculated  using  a  weighted-least-squares 
reconstruction  of  the  primitive  variables  while  viscous  fluxes  are  calculated  using  a 
deferred-correction  approach.  The  data  parallel  line  relaxation  numerical  method  is 
used  along  lines  of  cells  normal  to  walls  for  rapid  convergence  to  steady-state.  How¬ 
ever,  the  full  matrix  point  relaxation  numerical  method  is  used  in  regions  where  lines 
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cannot  be  formed.  US3D  incorporates  a  generalized  set  of  bonndary  conditions  in- 
clnding  catalytic  and  partially-catalytic  walls  with  and  withont  radiative  eqnilibrinm, 
wall  blowing  and  snction,  snbsonic  inflow  and  ontflow  conditions,  and  slip  wall  con¬ 
ditions  for  rarefied  flow  applications.  The  next  section  will  discnss  how  an  ablating 
wall  bonndary  condition  is  implemented  in  US3D. 

3.3  Ablation  Modeling 

3.3.1  Past  Work. 

Using  carbon  as  a  TPS  material  for  re-entry  applications  has  been  extensively 
examined.  Most  of  the  early  research  almost  exclnsively  focnsed  on  the  oxidation 
rate  of  carbon  nnder  varions  snrface  temperatnres  and  pressnres.  The  experimental 
data  always  assnmed  that  the  oxidation  process  followed  a  rate  law  in  an  Arrhenins 
form: 

m  =  Ap^  exp  {-Ea/RT)  (43) 

where  A  is  a  pre-exponential  coefficient,  p  is  the  snrface  pressnre,  Ea  is  the  activation 
energy,  R  is  the  nniversal  gas  constant,  and  T  is  the  snrface  temperatnre.  Note  that  m 
is  a  fnnction  of  the  pressnre  to  the  power  n,  the  order  of  the  reaction.  However,  hnge 
scatter  was  seen  in  the  experimental  data,  and  there  were  varying  levels  of  reasoning 
for  why  one  data  set  was  different  from  another.  Scala  [105]  attempted  to  specify 
a  set  of  kinetic  parameters  that  best  bracketed  the  data,  which  spanned  a  reaction 
rate  range  of  fonr  orders  of  magnitnde.  Unfortnnately,  all  of  the  data  nsed  by  Scala 
was  obtained  below  a  snrface  temperatnre  of  1373  K,  and  linear  extrapolation  of 
simple  Arrhenins  expressions  to  higher  temperatnres  is  known  to  predict  oxidation 
rate  probabilities  exceeding  those  theoretically  possible  [62]. 

Maahs  [62]  compiled  high  temperatnre  data  from  many  sonrces  to  demonstrate 
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Figure  12.  Arrhenius  plot  of  carbon  mass  loss  rate  for  a  carbon  filament-oxygen  reaction 
at  p  =  1.33  Pa  [62]. 


the  non-linear  behavior  of  carbon  oxidation.  Fignre  12  shows  the  mass  loss  rate  dne 
to  oxidation  from  a  carbon  hlament-oxygen  reaction  at  p  =  1.33  Pa.  The  dashed  line 
in  the  hgnre  represents  a  kinetic  expression  developed  by  Blyholder  et  al.  [16]  that 
described  the  snrface  as  having  two  types  of  reaction  sites,  one  being  more  reactive 
than  the  other.  Althongh  not  shown  in  Fignre  12,  the  experiments  of  Rosner  and 
Allendorf  [102,  103]  have  also  shown  the  oxidation  rate  in  an  Arrhenins  plot  to  be 
non-linear.  The  mass  loss  rates  differ  by  orders  of  magnitnde,  and  the  cnrves  exhibit 
maxima  and  minima.  Some  of  the  differences  can  be  attribnted  to  the  diversity 
of  carbon  materials  nsed  in  the  experiments  (i.e.  impnrities,  porosities,  strnctnre, 
etc.).  Clearly,  a  simple  linear  Arrhenins  expression  cannot  be  expected  to  adeqnately 
describe  the  oxidation  process  that  wonld  span  a  wide  temperatnre  and  pressnre 
range. 
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Many  researchers  have  attempted  to  theoretically  and  empirically  derive  kinetic 
expressions  for  carbon  oxidation  at  high  temperatnres.  Both  Blyholder  et  ah  [16] 
and  Strickland-Constable  et  ah  [78]  nsed  a  snrface  model  that  assnmed  reaction  sites 
conld  be  classihed  into  a  more  reactive  A-site  and  a  less  reactive  B-site.  At  low 
temperatnres,  the  oxidation  rate  was  dominated  by  the  A-sites.  As  the  temperatnre 
rose,  the  A-sites  thermally  annealed  to  form  less  reactive  B-sites.  The  B-sites  created 
A-sites  as  they  oxidized  and  desorbed  CO.  It  was  the  competition  between  the  genera¬ 
tion  and  deactivation  of  both  these  reaction  sites  that  resnlted  in  observed  maximnms 
in  the  oxidation  rate  [58],  An  important  assnmption  to  note  was  that  CO  was  the 
primary  prodnct.  Ong  [85]  nsed  transition-state  theory  to  dednce  a  set  of  controlling 
kinetic  steps,  calcnlating  reqnired  transition  energies  by  statistical  thermodynamics, 
and  empirically  determining  several  constants  from  high  temperatnre  oxidation  data. 
Ong’s  rate  expression  only  assnmed  one  type  of  active  site  on  the  carbon  snrface. 
A  contrasting  featnre  from  the  Blyholder  and  Strickland-Constable  models  was  that 
Ong  considered  two  prodnct  species  from  the  carbon-oxygen  reaction,  CO2  at  low 
temperatnres  and  CO  at  higher  temperatnres.  A  rate  maximnm  was  also  predicted 
by  Ong’s  model,  which  was  attribnted  to  a  change  in  the  vibration  chemical  potential 
of  molecnlar  oxygen. 

None  of  the  aforementioned  kinetic  models  accnrately  predict  oxidation  rates  be- 
canse  they  rely  on  linear  Arrhenins  expressions.  A  gas-snrface  reaction  can  be  very 
complex  considering  there  may  be  a  large  nnmber  of  separate,  competing  kinetic 
steps  involved.  These  inclnde  adsorption,  snrface  migration,  chemical  bonding  with 
a  carbon  atom,  and  desorption  of  the  gaseons  prodnct.  One  of  these  steps  may  be 
of  greater  or  lesser  importance  becanse  the  apparent  activation  energy  and/or  steric 
factor  may  change  as  the  snrface  temperatnre  and  pressnre  changes.  Therefore,  while 
it  may  be  possible  to  correlate  kinetic  data  with  a  simple  Arrhenins  expression  over 
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a  narrow  range  of  temperatnres  and  pressnres,  snch  an  expression  cannot  be  applied 
with  reliability  ontside  that  range  [62], 

Mnch  nncertainty  remains  in  predicting  mass  loss  rates  dne  to  snblimation  pro¬ 
cesses.  Snblimation  is  known  to  be  the  primary  mass  loss  mechanism  for  snrface 
temperatnres  above  approximately  3000  K.  The  snblimation  regime  has  typically 
been  dehned  as  the  range  of  conditions  where  the  mass  loss  dne  to  vaporization  ex¬ 
ceeds  the  diffnsion-controlled  oxidation  mass  loss  rate  [106],  In  order  to  predict  these 
rates,  thermodynamic  databases  are  commonly  nsed  to  calcnlate  transport  properties 
of  individnal  chemical  species  at  a  given  temperatnre.  However,  previons  research  has 
shown  that  these  databases  predict  large  variations  in  carbon  snblimation  rates.  Mi¬ 
los  and  Chen  [73]  compared  the  JANAF  [1]  and  CEA  [70]  thermodynamic  databases 
and  fonnd  that  C3  and  C5  contained  large  discrepancies.  The  C3  disagreement  was 
tronblesome  becanse  this  is  the  primary  species  present  for  carbon  snblimation  and  is 
critical  for  the  prediction  of  ablation  rates.  Havstad  and  Ferencz  [44]  recommended 
the  addition  of  C5  and  C7  snblimation  into  the  noneqnilibrinm  chemistry  model  for 
a  more  accnrate  ablation  rate  prediction  when  snrface  temperatnres  are  greater  than 
approximately  3900  K. 

Gosse  et  al.  [36]  compared  eqnilibrinm  ablation  prohles  from  the  CEA  database 
and  other  researchers  who  have  attempted  to  improve  the  eqnilibrinm  vapor  pressnres 
of  carbon  gas  species.  Eignre  13  shows  eqnilibrinm  ablation,  or  B'  cnrves,  for  differ¬ 
ent  thermochemical  models  in  the  snblimation  regime.  The  divergence  in  predicted 
ablation  rates  becomes  more  severe  as  snrface  temperatnres  rise  above  3400  K,  which 
highlights  the  discrepancy  between  the  different  thermochemical  models.  Addition¬ 
ally,  Gosse  et  al.  [36]  compared  the  pressnre-temperatnre  phase  change  diagram  for  a 
graphite  system  again  nsing  valnes  from  CEA  and  eqnilibrinm  vapor  pressnres  given 
by  noneqnilibrinm  snblimation  models  for  carbon.  Eignre  14  shows  a  large  variation 
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Figure  13.  Equilibrium  ablation  rate  curves  from  the  NASA  CEA  database  and  equi¬ 
librium  vapor  pressures  from  selected  finite-rate  chemistry  models  [36]. 


Figure  14.  Graphite  pressure-temperature  phase  change  diagram  comparing  different 
thermodynamic  data  [36]. 
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between  the  predicted  pressure-temperature  at  which  carbon  would  sublimate.  The 
CEA  database,  which  is  currently  assumed  to  give  the  best  thermodynamic  values, 
does  not  line  up  with  the  graphite  to  vapor  phase  change. 


3.3.2  Equilibrium:  B'  Curve  for  Air-Carbon. 

A  common  approximation  to  account  for  the  effects  of  ablation  products  being 
released  into  the  gaseous  flow  held,  and  accompanying  surface  blowing  that  occurs,  is 
to  employ  a  blowing  correction  term,  S',  that  is  based  on  boundary  layer  theory  and 
mass  balance  at  the  surface.  Consider  a  mixture  of  air  and  carbon  in  the  gas-phase 
interacting  with  a  surface  of  solid  carbon.  The  mass  hux  of  carbon  gas,  rhc,  from  the 
wall  is: 


PvjDc^  TJc 


T  nZc  Vcyw'^c 


W 


(44) 


where  r/c  is  the  total  mass  fraction  of  gaseous  carbon  and  the  subscript  w  represents 
the  wall  state. 

Now  assume  that  the  heat/mass  transfer  analogy  holds  such  that  Cm  =  Ch  with 
unity  Prandtl/Lewis  numbers  and  equal  species  diffusion.  Using  the  dehnition  of  the 
blowing  and  heat  transfer  coefficients,  it  can  be  shown  that  B'  is  given  by: 


B' 


m  m 

_  _ 

Pe^eC M  Pe'^eC'H 


yc,w 

1  Uc^w 


yc,w 

ya^w 


(46) 


where  the  subscript  e  represents  boundary  layer  edge  conditions  and  is  the  mass 
fraction  of  air  species  at  the  wall  conditions.  Thus,  B'  is  the  ratio  of  the  mass  fraction 
of  gaseous  carbon  species  to  air  species  at  the  local  wall  conditions  [23].  MacLean 
[63]  showed  that  the  blowing  correction  can  be  adjusted  to  match  experimental  data 
using  the  expression: 

Ch  _  2\B' 

^  _  1 
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where  Ch,o  is  the  non-blowing  heat  transfer  coefficient  and  A  is  an  empirical  parameter 
that  can  be  tailored  for  the  expected  ontpnt. 

To  obtain  the  eqnilibrinm  satnrated  state  for  an  air-carbon  mixtnre  at  a  given 
pressnre  and  temperatnre,  bnlk  carbon  is  added  to  air  nntil  the  carbon  starts  to 
condense  to  a  solid.  That  condensation  point  gives  the  eqnilibrinm  valne  for  B' . 
Fignre  15a  plots  B'  for  an  air-carbon  mixtnre  as  fnnction  of  temperatnre  for  several 
pressnres.  At  low  temperatnres,  air  reacts  with  carbon  to  prodnce  CO2,  which  resnlts 
in  a  platean  at  B'  =  0.089.  In  the  next  temperatnre  regime,  CO  is  the  dominant 
reaction  prodnct  dne  to  reactions  of  atomic  oxygen  with  carbon.  The  rate  of  ablation 
increases  nntil  all  of  the  oxygen  is  consnmed  resnlting  in  a  valne  of  B'  =  0.178.  As  the 
snrface  temperatnre  continnes  to  increase,  the  solid  carbon  begins  to  snblime  directly 
to  gaseons  carbon.  The  B'  of  this  process  increases  exponentially  as  the  snrface 
temperatnre  approaches  and  passes  the  eqnilibrinm  vapor  pressnre  of  the  carbon  gas 
mixtnre. 

Fignre  15b  plots  the  composition  of  the  satnrated  air-carbon  mixtnre  at  a  pressnre 
of  1  atm.  For  snrface  temperatnres  below  2500  K,  the  gas  is  a  simple  mixtnre  of 
N2,  CO,  and  CO2.  Above  2500  K,  appreciable  qnantities  of  CN  start  to  form.  As 
snblimation  becomes  stronger,  the  snrface  chemical  state  becomes  complicated  as 
many  polyatomic  carbon  and  nitrogen  species  are  formed.  The  primary  species  formed 
in  order  of  concentration  level  are  C3,  C5,  CN,  C2,  C4,  and  then  C. 

Many  ablation  modeling  approaches  nse  B'  to  obtain  the  gas-solid  interface  bonnd- 
ary  condition  [76,  8,  26,  73,  63].  For  each  gas-snrface  bonndary  of  the  CFD  grid,  the 
snrface  temperatnre,  pressnre,  and  existing  gas-phase  composition  are  nsed  as  inpnt 
variables.  The  resnlting  satnrated  air-carbon  composition  provides  mass  fractions  of 
all  gaseons  species  that  are  injected  into  the  flow.  The  new  mass  fractions  are  set 
as  explicit  bonndary  conditions  [63].  The  gas  blowing  rate  m  is  then  fonnd  nsing 
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(a)  B'  for  air-carbon  as  function  of  temperature  and  pressure. 


(b)  Mass  fractions  and  B'  at  1  atm  pressure. 

Figure  15.  Equilibrium  saturated  air-carbon  mixture  composition  and  dimensionless 
blowing  rate. 
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Equation  44.  The  heat  flux  can  also  be  reduced  using  B'  through  Equation  46  to 
account  for  the  effect  of  ablation. 

3.3.3  Finite-Rate  Surface  Reaction  Formulation. 

The  nonequilibrium,  finite-rate  gas-surface  interaction  system  model  utilized  in 
this  work  follows  the  formulation  developed  by  Marschall  and  MacLean  [67]  and 
MacLean  et  al.  [64].  A  summary  of  the  derivation  is  presented  to  provide  informa¬ 
tion  necessary  to  follow  the  implementation  discussion  in  the  next  section.  Please 
refer  to  the  identified  references  for  more  information.  The  model  consists  of  three 
environments  that  can  exist  at  the  gas-surface  interface: 

1.  Gas  environment 

2.  Surface  environment 

3.  Bulk  environment 

The  gas  environment  is  a  single  phase  containing  all  the  gas  species.  The  surface 
environment  can  have  multiple  phases,  each  occupying  a  fraction  of  the  total  surface 
area.  Each  surface  phase  consists  of  one  or  more  sets  of  “active  sites"  where  surface  re¬ 
actions  can  occur.  Each  set  of  active  sites  has  an  associated  set  of  chemically-distinct 
species.  The  surface  species  can  include  open  or  available  sites  and  chemically  ad¬ 
sorbed  species.  The  total  number  of  open  and  occupied  sites  for  each  set  is  conserved. 
The  bulk  environment  can  also  have  multiple  phases  with  each  representing  a  differ¬ 
ent  material  or  constituent.  Each  bulk  phase  (i.e.  solid  phase)  is  associated  with 
a  unique  set  of  chemical  species  that  can  participate  in  surface  reaction  processes 
as  would  occur  with  a  re-entry  vehicle  undergoing  ablation  [64].  All  reactions  must 
occur  at  active  sites  on  the  surface  and  involve  open  and/or  occupied  sites.  Note  that 
the  same  atom  or  molecule  is  treated  as  a  distinct  chemical  species  when  in  the  gas 
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or  surface  phase.  The  different  environments  recognize  that  surface  reactions  occur 
with  chemical  species  in  different  chemical  states  (i.e.  gas,  adsorbed,  bulk),  which 
provides  the  flexibility  to  specify  rate  parameters  in  an  appropriate  manner. 

A  surface  reaction  should  involve  surface  species  to  include  interactions  with  each 
other  and/or  species  from  the  gas  and  bulk  environments.  Additionally,  surface  reac¬ 
tions  must  involve  empty  and/or  hlled  sites.  The  following  surface  reaction  types  can 
be  dehned  where  (s)  denotes  an  open  active  site,  A(s)  denotes  an  adsorbed  species 
on  the  surface,  and  A{b)  denotes  a  bulk  phase  species: 


•  Adsorption/desorption 

•  Eley-Rideal  (ER) 

•  Langmuir-Hinshelwood  (LH) 

•  Oxidation/reduction 

•  Sublimation/condensation 


A  +  (s)  ^  A(s) 

A  +  B{s)  ^  AB  +  (s) 

A(s)  +  B{s)  AB  +  2{s) 
A+{s)  +  B{b)  ^  AB  +  {s) 
(s)  +  A(b)  ^  A  +  (s) 


The  only  bulk  phase  material  considered  in  this  work  was  that  of  carbon,  so  the 
number  of  bulk  species  is  equal  to  one.  The  production  or  loss  of  any  species  comes 
from  the  summation  of  the  forward  and  backward  rate  processes  for  each  reaction. 
The  procedure  is  analogous  to  the  hnite-rate  formulation  for  the  gas  phase  chemistry 
discussed  in  Section  3.2.5,  except  species  production  occurs  at  the  surface  and  the 
rate  is  per  unit  surface  area  per  unit  time  (e.g.  kmol  m“^  The  forward  and 

backward  rate  coefficients  are  related  through  an  equilibrium  constant,  which  can  be 
obtained  in  three  ways: 


1.  Specifying  individual  Gibbs  free  energies  for  each  species 


2.  Explicitly  specifying  the  equilibrium  constant  by  use  of  a  curve  ht  (i.e.  Equa¬ 
tion  31)  or  a  direct  expression 
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Figure  16.  Surface  mass  balance. 


3.  Specifying  forward  and  backward  reaction  rate  coefficients 
The  reaction  flnx,  r*,  for  reaction  i  is  given  by: 

K  ,  K  „ 

n^k„Y[Xl‘“-kuY[Xl‘“  (47) 

k=l  k=l 

where  kfi  and  kbi  are  the  forward  and  backward  reaction  rate  coefficients  for  reaction 
i  at  a  given  temperatnre,  is  the  reactant  stoichiometric  coefficient  for  species  /c, 
is  the  prodnct  stoichiometric  coefficient  for  species  k,  and  Xk  is  a  generalized 
concentration  of  species  k.  The  concentration  of  species  k  is  nsnally  calcnlated  as 
Pk/Mk  as  shown  in  Eqnation  25.  The  reaction  flnx  is  given  in  nnits  of  (kmol  m“^ 
s“^).  The  net  prodnction  rate  of  species  k  is: 

NR 

WZc  =  ^  {l^kiTi)  (48) 

i=l 

where  z/^j  =  (z/(!-  —  z/(,J  and  NR  represents  the  total  nnmber  of  reactions  involving 
species  k.  The  chemical  sonrce  term  in  Eqnation  48  applies  to  any  species  at  the 
gas-snrface  interface. 
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Mass  balance  at  the  surface  is  used  to  set  the  boundary  condition  for  each  species 
in  the  gas  phase  by  considering  the  fluxes  of  mass  entering  and  leaving  an  infinitely 
thin  control  volume  fixed  at  the  surface  as  shown  in  Figure  16.  The  surface  mass 
balance  for  species  k,  with  the  assumption  that  no  material  is  being  removed  due  to 
mechanical  erosion  (spallation),  can  be  written  as: 


PwDk^  Vk 


TTffciC/j  pw'^wVkjV 


(49) 


where  yk  is  the  species  mass  fraction  and  is  the  surface  normal  velocity.  The  terms 
on  the  left-hand  side  of  Equation  49  represent  the  mass  fluxes  entering  the  surface 
due  to  diffusion  and  species  production  from  surface  reactions.  These  mass  fluxes 
are  exactly  balanced  by  the  terms  on  the  right-hand  side  that  represent  the  rate  at 
which  the  mixture  is  convected  away  from  the  surface  due  to  ablation.  Equation  49 
represents  the  rate  at  which  the  mass  of  each  gas  phase  species  changes  at  the  surface 
|64|. 


3.3.4  Surface  Equilibrium  Constants. 

Calculation  of  the  reaction  flux  r*  in  Equation  47  requires  the  specification  of  two  of 
the  three  quantities:  kfi,  kbi,  or  iCeg-  Typically,  the  forward  rate  coefficient  is  specified, 
and  the  backward  rate  coefficient  is  calculated  using  the  equilibrium  constant.  The 
calculation  of  equilibrium  constants  for  surface  reactions  requires  thermodynamic 
data  for  each  surface  and  gas  phase  species  in  the  reaction.  If  the  backward  rate 
coefficient  is  specified,  the  concentration-based  equilibrium  constant  can  be  calculated 
directly  using  Equation  28.  Otherwise,  it  can  be  computed  using  the  activity-based 
equilibrium  constant  as  shown  in  Equation  29.  In  the  calculation  of  the  activity- 
based  equilibrium  constant,  the  enthalpy  and  entropy  are  not  available  for  species  in 
the  surface  phase,  so  equilibrium  constant  expressions  are  specified  directly  for  the 
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adsorption/desorption  reactions. 

The  adsorption/desorption  equilibrium  constants  are  expressed  in  dimensionless 
activity-based  form  and  depend  on  the  type  of  adsorption  (mobile  or  immobile).  For 
mobile  adsorption,  the  adsorbed  species  remains  in  a  gaseous  state  and  acts  as  a  two- 
dimensional  gas  over  the  surface.  It  is  common  to  assume  that  mobile  adsorption  has 
no  activation  energy  since  there  is  a  very  weak  van  der  Waals  bond  to  the  adsorption 
sites.  The  activity-based  equilibrium  constant  is  given  by: 


K 

a,m 


-1/2 


exp 


^  des 

~f~ 


(50) 


where  po  is  a  reference  pressure  of  10®  Pa,  rrii  =  mo  or  m^r,  h  is  the  Planck  con¬ 
stant,  and  Tdes  is  related  to  the  desorption  energy  by  Edes  =  RTdes-  The  surface  site 
concentration  B  represents  the  total  number  of  available  reaction  sites  given  in  m“^. 
Substituting  k  =  R/A^  and  =  Mi/ Ay  gives: 


_  T,po  f27rM,RTY^^\_  fTdes\ 
BRT  V  Alh^  ) 


(51) 


where  Ay  is  Avogadro’s  number.  For  immobile  adsorption,  the  adsorbed  species  stays 
wherever  it  was  adsorbed  until  it  desorbs  or  reacts  with  the  bulk  surface.  This  process 
requires  an  activation  energy  for  the  gas-phase  species  to  overcome  in  order  to  bond 
to  an  adsorption  site.  The  activity-based  equilibrium  constant  is  given  by: 

..  _  f2nM,RTy^/^  fTdes\ 

^a,im  ^2^2  J  [  T  J 


The  activity-based  equilibrium  constant  is  converted  to  a  concentration-based  equi¬ 
librium  using  Equation  29.  For  an  adsorption  reaction,  Ug  equals  -1  and  gives  the 
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concentration-based  eqnilibrinm  constant  for  mobile  adsorption  as: 


K 


eq,m 


Ay 

l3 


2irM,RTY''^ 


(B3) 


and  for  immobile  adsorption: 


K. 


eq^ira 


Ay 


2irMiRTY'‘'‘‘ 


(B4) 


Marschall  and  MacLean  [67]  provided  a  general  eqnation  for  the  adsorption/desorp¬ 
tion  reactions  that  can  be  implemented  into  a  CFD  code: 


Keq  =  AegT^Ko  exp 


Edes  E, 


ad 


RT 


(BB) 


where  A^g  is  a  constant,  f3  and  T'  are  dimensionless  {T'  =  T/IK),  and  Edes  and  Ead 
are  the  activation  energies  for  desorption  and  adsorption,  respectively.  Relating  the 
above  constants,  it  is  shown  that: 


A 


eq  — 


mobile 

immobile 


/27rM,RTy^/^ 

/27rMiRTY^^^ 


mobile 

immobile 


Therefore,  the  type  of  adsorption  process  is  critical  to  characterizing  the  overall  snr- 
face  kinetic  model.  The  sensitivity  of  the  ablation  process  to  the  adsorption  process 
is  discnssed  in  Chapter  IV. 
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Table  7.  Forward  reaction  rate  formulas. 


Type 

Arrhenius 

Adsorption 

Eley-Rideal 


Rate  Formula 


kf  = 
kf  = 


r  1 

SoT'^exp 

fill's 

■  Vs  - 

JerT'^  exp  ( 

Ead\ 

W  ) 

Eer  \ 

'WJ 


Specified 

Parameters 

A,f3,E 
*5*0)  -Fad 

^eri  Pi  Egr 


3.3.5  Forward  Reaction  Rates. 

Forward  reaction  rates  must  be  specified  for  each  surface  reaction  in  the  model. 
Please  refer  to  Reference  [67]  for  the  rate  formulas  available  for  the  different  reaction 
types.  The  implemented  surface  models  only  use  the  Arrhenius,  adsorption,  and  ER 
reaction  types  and  the  corresponding  rate  formulas  are  specified  in  Table  7. 

The  units  of  A  in  the  Arrhenius  formula  vary  depending  on  the  type  of  reactants 
and  stoichiometries  involved.  For  example,  the  units  of  A  would  be  (m^  krnol”^  s“^) 
for  a  reaction  involving  one  gas  and  one  surface  species.  For  a  reaction  involving 
two  surface  species,  the  units  of  A  would  be  (m^  krnol”^  s“^)  .  Finally,  the  units 
of  A  would  be  (s“^)  for  a  reaction  involving  one  surface  and  one  bulk  species.  The 
Arrhenius  expression  is  the  most  general  way  to  specify  a  forward  reaction  rate. 
The  disadvantage  of  using  the  basic  Arrhenius  formulation  is  that  the  magnitudes 
of  A  and  P  are  difficult  to  relate  to  physical,  chemical,  and  kinetic  processes  [67]. 
However,  kinetics-based  formulations  (i.e.  adsorption,  ER)  provide  a  more  intuitive 
way  to  specify  surface  reaction  rates  because  they  are  describing  physical  processes 
occurring  at  the  gas-surface  interface. 

The  leading  term  in  the  brackets  for  the  kinetics-based  rates  contains  the  following 
variables:  is  the  total  active  site  density,  is  the  site  density  exponent  equal  to  the 
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sum  of  stoichiometric  coefficients  for  all  surface  reactants,  and  Ug  is  the  mean  thermal 
velocity  of  the  incoming  reactant  gas  species: 


(B6) 


The  sticking  coefficient  So  and  ER  reaction  efficiency  coefficient  7er  are  dimensionless, 
and  their  values  should  he  between  zero  and  one. 

3.3.6  Park  Model. 

Two  of  the  more  common  air-carbon  surface  models  contained  in  the  literature  are 
from  Zhluktov  and  Abe  (ZA)  [136]  and  Park.  As  MacLean  et  al.  [64]  point  out,  the 
simple  Park  model  was  actually  hrst  published  by  Chen  and  Milos  [25].  This  “Park 
model”  is  more  or  less  a  compilation  of  Park’s  published  work  on  carbon  ablation 
models  [89,  90,  92,  95].  The  surface  reactions  involved  in  the  Park  model  are  given 
as: 


(1)  O  +  (s)  +  C(6)  ^  CO  +  (s) 

(2)  O2  +  2(s)  +  2C(6)  ^  2CO  +  2(s) 


(3)  N  +  (s)  +  C(6)  ^CN  +  (s) 


(4)  3(s)  +  3C(6)  ^C3  +  3(s) 

(5)  C3  +  3(s)  ^3(s)  +  3C(6) 


The  above  reaction  set  contains  irreversible  oxidation  of  both  atomic  and  molec¬ 
ular  oxygen,  irreversible  nitridation  of  atomic  nitrogen,  and  C3  sublimation.  The 
reaction  rates  are  based  on  kinetic  theory  [92]  with  experimentally  determined  re¬ 
action  probabilities  assigned  for  each  surface  reaction.  The  fluxes  for  each  reaction 
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are: 


ri 

__  Me 

—  ProVO^Oll  ,  , 

Mo 

(B7) 

r2 

_  o  -  Me 

^Pwyo2^02'y2 

1VI02 

(58) 

rs 

^  Me 
—  PwVnt^nIz  , , 

IVIn 

(59) 

r4 

=  pwyc3,E^C3'yA 

(60) 

r5 

=  PwyC3^C3l5 

(61) 

Here,  Ug  is  defined  as  and  7j  is  the  efficiency  of  snrface  reaction  i.  The 

reaction  efficiencies  are  taken  from  the  literatnre  to  be  [89,  25]: 


7l 

-  0.63  exp  ^  ) 

(62) 

72 

=  0.5 

(63) 

74 

=  1 

(64) 

75 

=  0.1 

(65) 

There  has  been  a  lot  of  recent  work  in  determining  the  reaction  probability  or 
efficiency  of  carbon  nitridation  becanse  it  is  considered  as  an  important  reaction  for 
carbon-based  TPS  materials.  The  nitridation  reaction  contribntes  directly  to  the 
snrface  mass  loss  and  radiative  heating.  The  efficiency  of  this  reaction  pnblished  by 
Park  and  Bogdanoff  [95]  states  that  30%  of  N  atoms  impacting  the  snrface  form  CN. 
However,  Driver  and  MacLean  [28]  fonnd  that  an  efficiency  of  1  x  10“^  provides  the 
best  fit  for  recession  data  in  arc  jet  tests.  Zhang  et  al.  [135]  provides  an  excellent 
comparison  of  experimental  work  determining  the  nitridation  efficiency,  inclnding 
their  own,  and  fonnd  the  efficiency  to  be  abont  two  orders  of  magnitnde  lower  than 
that  obtained  by  Park  and  Bogdanoff  [95].  The  experimental  efficiencies  generated 
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Figure  17.  Experimentally  determined  reaction  efficiency  for  carbon  nitridation  by 
atomic  nitrogen  from  different  investigators.  Reproduced  from  Zhang  et  al.  [135]  and 
references  cited  in  the  plot  refer  to  those  contained  in  that  article.  The  “o”  symbols 
refer  to  the  work  performed  by  Zhang  et  al. 


by  Goldstein  [33],  Suzuki  et  al.  [114,  116]  and  Zhang  et  al.  [135]  were  all  around 
10“^  and  showed  a  weak  temperature  dependence  as  shown  in  Figure  17.  Suzuki  et 
al.  [115]  proposed  a  nitridation  efficiency  based  on  their  results  as: 


jnit  =  8.441  X  10  ^  exp 


(66) 


The  above  efficiency  expression  hts  nitridation  data  with  surface  temperatures  rang¬ 
ing  from  1351  to  2184  K.  In  this  dissertation,  results  will  be  presented  using  the 
Park  ablation  model  with  a  nitridation  efficiency  of  0.3  and  with  that  specihed  by 
Equation  66. 

The  sublimation  reactions  are  based  on  a  Knudsen-Langmuir  formulation  for 
nonequilibrium  surface  evaporation  where  the  C3  concentration  in  Equation  60  is 
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computed  from  the  saturated  vapor  pressure  [8,  48],  The  reactiou  flux  theu  becomes: 


r4  =  Py^4 


M, 


Cs 


27tRTu 


(67) 


with  the  saturated  carbou  vapor  pressure  as: 


Pv 


5.19  X  10^^  exp 


-93310\ 

) 


(68) 


Species  mass  couservatiou  at  the  surface  is  theu  solved  usiug  Equatiou  49  with: 

ns 

PwVw  =  ^  rhs  (69) 

The  source  terms  for  each  species  are  dehued  as: 


mco  = 

,  ,  ^  Mco 

Me 

(70) 

rhcN  = 

Mcn 

Me 

(71) 

mcs  = 

r4  -  rs 

(72) 

rriN  = 

Mn 

"‘"'Mc 

(73) 

rho  = 

1 

o  o 

(74) 

mo2  = 

M02 

—To - 

2Mc 

(75) 

aud  Tils  =  0  for  all  other  species  cousidered. 


3.3.7  Zhluktov  and  Abe  Model. 

Cousidered  oue  of  the  most  advauced  kiuetic  models  for  air-carbou  gas-surface 
chemistry,  the  huite-rate  Zhluktov  aud  Abe  [136]  model  iuvolves  12  surface  reactious 
with  dehued  forward  aud  backward  rates  that  are  thermodyuamically  coustraiued. 
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Table  8.  Reaction  set  for  the  Zhluktov  and  Abe  gas-snrface  chemistry  model. 


Reaction 

Type 

So  /  7er.  /  A 

/3 

E  (kJ  mol  ^) 

1.)  0  +  (s)  ^O(s) 

Ads 

1 

0 

0 

2.)  N  +  {s)^  N{s) 

Ads 

1 

0 

0 

3.)  20{s)  ^  O2  +  2(s) 

Arrh 

3.58  X  10^3 

1 

256.07 

4.)  02  +  (s)  ^0  +  0(s) 

ER 

1 

0 

118.06 

5.)  CO2  +  (s)  ^CO  +  0{s) 

ER 

0.9 

0 

0 

6.)  0(s)  +  C'(6)  ^C'0  +  (s) 

Arrh 

2.08  X  10^ 

1 

332.56 

7.)  0  +  0{s)  +  C{b)  ^  CO2  +  (s) 

ER 

0.8 

0 

16.63 

8.)  20{s)  +  C{b)^C02  +  2{s) 

Arrh 

3.58  X  10^^ 

0 

332.56 

9.)  ^  +  (5)  ^  {s)  +  Cib) 

ER 

0.24 

0 

0 

10.)  C2  +  2{s)  ^  2{s)  +  2C{b) 

ER 

0.5 

0 

0 

11.)  (^3  +  3(5)  ^3(s)  + 3(7(6) 

ER 

0.023 

0 

0 

12.)  N2  +  {s)^N  +  N{s) 

ER 

1 

0 

636.85 

The  implementation  into  US3D  was  first  described  by  Candler  [23],  The  original 
implementation  had  to  be  modified  in  this  work  due  to  errors  found  through  initial 
testing.  Verification  of  the  implementation  in  US3D  is  supplied  in  Chapter  IV.  The 
approach  involves  solving  for  the  surface  coverage  (number  of  adsorbed  N  and  O 
atoms  on  the  available  surface  sites),  and  then  computing  the  rate  of  formation  of 
gas-phase  species  through  competing  kinetics-based  processes  [23].  The  reactions  and 
rate  coefficients  that  comprise  the  ZA  model  are  taken  from  Reference  [64]  and  are 
shown  in  Table  8. 

The  first  two  reactions  represent  the  adsorption  and  desorption  of  atomic  oxygen 
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and  nitrogen.  Reactions  3-8  and  12  represent  the  recombination  of  adsorbed  atoms, 
recombination  of  gas-phase  species  with  adsorbed  atoms,  and  oxidation  of  bnlk  car¬ 
bon.  Reactions  9-11  represent  the  snblimation  of  bnlk  carbon  or  the  condensation  of 
atomic  and  molecnlar  carbon  gas.  A  major  difference  from  the  Park  model  is  that 
snrface  reactions  are  allowed  to  take  place  in  both  directions  with  the  preferred  direc¬ 
tion  dependent  on  the  kinetic  rates  and  eqnilibrinm  constants  of  each  reaction.  The 
original  ZA  model  lacks  a  carbon  nitridation  reaction,  bnt  a  method  to  implement 
it  is  discnssed  in  Section  3.4.  The  importance  of  nitridation  will  be  discnssed  in  the 
next  few  chapters. 

Corresponding  net  reaction  flnxes  are  as  follows: 

ri  = 

r2  = 
rs  = 

r4  = 

r5  = 
re  = 
ri  = 
rg  = 
rg  = 
rio  = 
rn  = 
ri2  = 

where  <he,  d’oi  and  are  the  concentrations  of  empty  snrface  sites,  sites  with  ad- 


Mo 

kf2T^^e  —  kb2^N 


‘M. 


N 


kfs  (*ho)^  -  kbgj^  (*he)^ 
M02 

kfA-^^e  -  kbA-^^O 


M02 
PCO2 


o 


kfei^^e-kbb-^^o 


Mco2  Mco 

Pco 


kfe^o  ~  kbe 


M, 


CO 


1  Po  ^  J  PCO2  ^ 
kfj  . ,  d>r)  —  kh7^^ — ‘hf. 


M, 


■O  ~  i^bl 


O 


M, 


CO2 


{■ief 

^vico2 

Pc 

kfgj^^e  —  kbg^e 
Me 


k 


fl2 


Me, 

PN2 


-  k 


612 


N2 


Pn 

W 


<!> 


N 


N 
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Table  9.  Adsorption/desorption  equilibrium  rate  constants  for  the  Zhluktov  and  Abe 
model. 


Reaction 

^eq 

p 

Edes  (kj  mol  ^) 

lb.) 

0 

+ 

(s)  ^  0(s) 

1.72  X  10^ 

0 

374.13 

2b.) 

N 

+ 

(s)  ^  N{s) 

1.72  X  10^ 

0 

304.29 

sorbed  atomic  oxygen,  and  sites  with  adsorbed  atomic  nitrogen,  respectively.  The 
calcnlation  of  the  eqnilibrinm  constants  begins  by  nsing  known  eqnilibrinm  expres¬ 
sions  for  reactions  one  and  two.  The  eqnilibrinm  constant  for  reaction  one  is  a  fnnction 
of  the  concentration  of  oxygen  atoms  Xq  in  the  gas  phase  on  the  reactants  side  and 
the  amonnt  of  adsorbed  O  atoms  on  the  prodncts  side: 


kbl  Xo^e 


(77) 


where  X^  represents  the  concentration  of  species  k.  Solving  for  <he/ gives: 


*ho 


1  $o 

K^Xo  <he 


KiXo 


(78) 


The  <ho/*he  fraction  is  common  to  many  of  the  eqnilibrinm  constant  expressions  for 
the  reactions  listed  in  Table  8.  Similarly  for  reaction  two: 


<h7V 


^  _  kf2  _ 
kb2 

1 

K2Xn  ^ 


K2Xm 


(79) 


The  eqnilibrinm  constants  for  both  adsorption/desorption  reactions  {Ki  and  K2)  are 
fonnd  nsing  Eqnation  55  with  their  respective  rate  constants  and  activation  energies 
fonnd  in  Table  9.  Mobile  adsorption  was  assnmed  with  the  snrface  site  concentration 
set  to  B  =  3.5  X  10^^  m“^,  which  gives  =  5.8  x  10“®  kmol  m“^  [136]. 
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Next,  the  equilibrium  constant  for  reaction  three  can  be  written  as: 


Ks 


X02  1 


(80) 


The  ratio  Xq^I^’o  the  equilibrium  constant  for  the  gas  phase  reaction: 


0  +  0  — >  O2 


(81) 


The  equilibrium  constant  can  be  computed  using  thermodynamic  data  from  the  CEA 
database  and  the  Gibbs  free  energy  equation.  Hence,  is  easily  computed  from 
known  quantities.  Following  this  same  procedure,  the  rest  of  the  equilibrium  constants 
are  calculated  as  follows: 


K4 

K, 

K, 

Ks 

K9 

Kn 

Ku 


Xo^o 
Xo2^e  " 
Xco^o 
Xc02^e 
Xco^e 
Xc{b)^0 


Xo, 

XcoXq 

Xco2 

Xco 


Ki 


XoXcib)  K, 


Xc02^e  _  Xc02  1 
XoXcib)^o  ~  K, 

Xco2^l  ^  XcOj  1 
Xc(,)d>S  “  XlXcib)  Kf 

Xc{b)^e  _  Xc(b) 

Xc^e  ~ 


Xc2n 

x°u+l 

XcM 

Xn^n  _ 
XN2^e 


■^2 

■  Xc, 

Xc, 

X,. 


X 


N2 


(82) 


The  backward  rate  coefficients  are  simply  found  by  evaluating  ku  =  kfi/Ki.  It  is 
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important  to  stress  that  the  surface  equilibrium  constants  are  dependent  on  gas- 
phase  equilibrium  constants,  and  thus,  there  is  a  coupling  between  the  gas-phase 
kinetic  model  and  the  surface  boundary  condition  [23],  There  are  hve  reactions  in  the 
above  equilibrium  relations  that  involve  gas  phase  species  and  carbon  bulk  material 
(e.g.  Xc{b)/^c',  C  —>■  C{b)).  Fortunately,  the  CEA  database  includes  thermodynamic 
properties  for  certain  bulk  materials  such  as  carbon,  and  the  equilibrium  constants 
can  be  computed. 

The  rates  of  species  production  on  the  surface  are: 


mo  = 

(-ri  +  ri-  rr)  Mq 

rhco  = 

(rs  re)  Mco 

mco2  = 

(-rs  -f  ry  -k  rg)  Mcoa 

me  = 

-rgMc 

mc2  = 

-rioMc2 

mc3  = 

-riiMcg 

TTIn  = 

(-ra  +  ria) 

mo2  = 

(rg  -  r4)  M02 

171^2  = 

-riaMATj 

The  surface  production  rate  for  any  species  included  in  the  gas  phase  chemistry, 
including  charged  particles,  that  does  not  contribute  to  any  surface  reaction  is  set  to 
zero.  The  mass  fluxes  for  the  amount  of  O  and  N  absorbed  on  the  surface  are: 

mo(s)/Mo  =  ri  -  2r3 r4 -t- rs  -  re  -  ry  -  2r8  =  0 

mN{s)/MN  =  r2  +  ri2  =  0  (84) 
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Using  these  expressions,  the  total  snrface  mass  blowing  rate  is: 


m  =  ^  rhj  =  Me  (rg  +  ry  +  rg  -  rg  -  2rio  -  Srn)  (85) 

The  solntion  of  Eqnations  84  together  with  the  snm  of  snrface  coverage  concentrations: 

To  +  (86) 


yields  To,  and  Te  valnes. 

Newton’s  method  for  non-linear  systems  in  two-dimensions  is  nsed  for  determining 
the  snrface  coverage  concentrations,  which  is  written  in  vector  notation  as: 


P„+l  =  Pn-  /■'  (A)  F  (F„) 


(87) 


where  J  represents  the  Jacobian.  The  vectors  P  and  F  are  given  by: 


To 

fi  (<^>o) 

ri  -  2r3  +  r4  -1-  rs  -  re  -  r7  -  2r8 

,  F  = 

Tat 

/2  (*hAr) 

r2  +  ri2 

(88) 


The  Jacobian  and  its  inverse  are  written  as: 


r  dh 

dh  1 

dh 

dh  1 

J  = 

d^o 
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dh 
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Last,  the  partial  derivatives  of  fi  and  /2  need  to  be  calcnlated  nsing  the  chain  rnle 
by  invoking  the  constraint: 


Te  =  T 
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aT, 


aT 
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=  -l, 


aT, 


aT 


=  -1 


N 


(90) 
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The  solution  procedure  is  iterated  until  we  have  converged  solutions  for  $0  and  $7v 
By  substituting  $0,  and  <he  into  Equations  76,  the  necessary  mass  production 
rates  can  be  obtained  from  Equations  83.  Next,  the  surface  mass  balance  equation 
is  solved  where  is  obtained  from  Equation  85  since  rh  =  p^v^.  Equation  49  is 
non-linear  that  is  solved  iteratively  for  each  species  until  the  surface  normal  velocity 
converges.  The  approach  assumes  that  Tick’s  law  for  diffusion  holds. 

As  an  example,  the  surface  mass  balance  equation  can  be  solved  for  the  N2  surface 
mass  fraction.  The  discretized  form  of  Equation  49  is  given  as: 


,D 


(l/jVa.l  -  yN2,w) 

An 


+  PwVwyN2,w  —  —  'I^N2 


(91) 


where  7/^2, 1  is  the  mass  fraction  of  N2  in  the  hrst  cell  away  from  the  wall  and  An  is  the 
distance  from  the  hrst  cell  center  to  the  wall.  Substituting  for  171^2  from  Equations  83 
and  76  and  rearranging  terms,  results  in: 


D 

^  )  yN2,z 

D 

^  )  yN2,z 


-^yN2,i  4 - i—ri2MN  ) 

An  p^ 

D  ,  M7V2  Pn  ^  , 

+  kbl2^ - —  k 


Mjv  Pu 


fl2^—^^e  (92) 

Pn 


where  Pn/Pw  =  yN,w  and  PN2/ Pw  =  yN2,w  Solving  for  yN2,w  gives: 


yN2,w 


D  M]\f^ 

N—yN2,l  +  kbl2N-;—yN,w^N 
An  M]s[ 

D 

+  Vw  +  kfi2^e 


(93) 


Equation  93  is  an  approximate  solution  for  yN2,w  because  we  are  assuming  yN,w  is 
known  when  it  actually  has  its  own  mass  balance  equation  to  solve.  However,  the 
entire  solution  procedure  is  iterated  until  the  wall  normal  velocity  converges.  All  mass 
fractions  at  the  wall  should  balance  when  the  convergence  criterion  for  is  met  and 
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then  the  proper  gradients  for  the  wall  bonndary  condition  can  be  set. 

3.3.8  Comparison  of  Park  and  ZA  Models. 

Many  anthors  have  applied  and  compared  the  Park  and  ZA  models  to  different 
re-entry  trajectories  and  arcjet  tnnnel  data.  Havstad  and  Ferencz  [44]  fonnd  the 
models  predicted  similar  ablative  mass  flnxes  nnder  snblimation  conditions.  The 
primary  species  resnlting  from  carbon  snblimation  is  C3,  and  both  models  inclnde 
a  C3  reaction  with  similar  forward  rates.  Chen  and  Milos  [27]  showed  that  the  ZA 
model  predicted  snbstantially  lower  stagnation  point  ablation  rates  compared  to  the 
Park  model  bnt  larger  convective  heat  flnxes.  Another  interesting  discovery  by  Chen 
and  Milos  was  the  effect  of  the  nitridation  reaction  inclnded  in  the  Park  model,  which 
the  ZA  model  lacks.  The  removal  of  the  nitridation  reaction  from  the  Park  model 
signihcantly  lowered  the  ablation  rate  to  valnes  near  the  ZA  prediction  bnt  had  a 
relatively  minor  impact  on  other  parameters  snch  as  the  wall  heating  rate. 

Snznki  et  al.  [115,  116]  calcnlated  the  nitridation  rate  of  graphite  test  models  to 
be  100  times  smaller  than  the  valne  obtained  by  Park  and  Bogdanoff  [94]  in  plasma 
wind  tnnnel  tests  at  snrface  temperatnres  ranging  from  1351-1723  K.  However  even 
with  a  mnch  lower  nitridation  efficiency,  they  predicted  the  snrface  recession  of  the 
HAYABUSA  capsnle  to  be  25%  larger  than  that  calcnlated  withont  nitridation.  Beer- 
man  et  al.  [12]  nsed  the  Park  model  with  and  withont  nitridation  to  analyze  a  carbon 
ablative  material  at  a  Stardnst  peak  heating  condition.  They  fonnd  that  nitridation 
had  a  very  minimal  impact  on  the  stagnation  point  ablation  and  convective  heating 
rates,  which  is  in  stark  contrast  to  the  resnlts  of  Chen  and  Milos  [27]  and  Snznki  et  al. 
[115,  116]  The  rationale  was  that  nitridation  was  an  insignfficant  snrface  reaction,  or 
did  not  occnr,  and  it  conld  be  replaced  by  a  snrface  nitrogen  catalytic  process.  This 
idea  was  carried  over  into  the  work  of  Chen  and  Gokcen  [24]  who  compared  the  two 
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models  against  arcjet  data  but  replaced  Park’s  nitridation  reaction  with  a  nitrogen 
surface  recombination  reaction.  Note  that  the  arcjet  tests  had  stagnation  pressures 
ranging  from  approximately  0.023-0.83  atm  and  surface  temperatures  ranging  from 
approximately  1600-3000  K.  It  will  be  shown  in  this  research  that  the  level  of  nitrida¬ 
tion  is  dependent  on  the  wall  pressure  and  temperature.  The  ablation  rates  predicted 
by  the  modihed  Park  model  matched  arcjet  data  extremely  well  for  a  majority  of 
the  tests  while  the  ZA  model  consistently  under  predicted  the  ablation  rates.  The 
inclusion  of  a  nitridation  reaction  into  the  ZA  model  could  possibly  improve  these 
predictions. 

Recognize  that  both  ablation  models  specify  forward  rates  using  relatively  linear 
Arrhenius  expressions.  As  previously  discussed  and  shown  in  Figure  12,  measured 
oxidation  rates  displayed  non-linear  Arrhenius  behavior.  The  ZA  model  attempts  to 
account  for  this  non-linearity  by  specifying  “linear”  expressions  for  many  separate, 
competing  kinetic  steps.  The  surface  reactions  specihed  in  the  legacy  models,  includ¬ 
ing  the  Park  model,  only  describe  one-step  kinetic  processes  and  cannot  adequately 
capture  this  non-linear  behavior. 

3.4  Carbon  Nitridation  Reaction 

Carbon  nitridation  is  recognized  as  an  important  surface  reaction  to  characterize 
in  order  to  properly  compute  convective  and  radiative  heat  fluxes  [108,  38].  Currently, 
the  ZA  model  can  only  be  considered  an  oxidation  and  sublimation  surface  reaction 
model  because  it  does  not  include  a  carbon  nitridation  mechanism.  The  addition 
of  a  carbon  nitridation  reaction  will  make  it  a  more  comprehensive  and  physically 
authentic  model.  The  development  of  the  nitridation  reaction  rates  for  inclusion  in 
the  ZA  model  shown  in  this  section  represents  new  and  original  work. 

Zhang  et  al.  [135]  presented  a  review  experimental  carbon  nitridation  reaction 
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efficiencies  determined  from  measnred  carbon  mass  loss  rates.  Additionally,  Lntz  et 
al.  [60]  is  attempting  to  measnre  nitridation  and  oxidation  efficiencies  nsing  laser- 
indnced  flnorescence.  The  nitridation  reaction  defined  in  these  works  is: 

N  +  C(b)  CN  (94) 

There  have  been  attempts  to  implement  this  reaction  into  compntational  ablation 
models,  however  it  was  nsnally  dehned  as  an  irreversible  process  [95,  25,  6],  Following 
the  snrface  chemistry  model  developed  in  the  previons  section,  the  carbon  nitridation 
reaction  can  be  represented  nsing  an  ER  recombination  process  or  an  Arrhenins- 
type  expression.  All  dehned  snrface  reactions  are  assnmed  to  proceed  throngh  a 
single  active  site  type  (i.e.  bnlk  carbon).  The  hrst  path  is  simply  dehned  as  direct 
nitridation: 

N  +  (s)  +  C(b)  ^  CN  +  (s)  (95) 

The  second  path  is  dehned  as  a  snrface  participating  reaction: 

N(s)  +  C(b)  ^  CN  +  (s)  (96) 

These  formnlations  allow  for  the  snrface  reactions  to  occnr  in  forward  and  backward 
directions  with  the  preferred  direction  dependent  on  the  kinetic  rates  and  eqnilibrinm 
constants  of  each  reaction. 

The  net  reaction  hnx  for  each  pathway  can  be  expressed  in  terms  of  forward  and 
backward  reaction  rate  coefficients  (kf  and  kb),  stoichiometric  coefficients  (z/),  and 
the  generalized  concentrations  of  gas,  snrface,  and  bnlk  species  (X^)  as  shown  in 
Eqnation  47.  Eor  a  gas  phase  species  Xg  =  Cg  =  Ps/Mg  in  kmol  m“^,  for  a  snrface 
species  Xg  =  in  kmol  m“^,  and  for  a  pnre  bnlk  species  X*  =  1  (dimensionless). 
The  net  reaction  hnx  is  different  between  the  reaction  paths.  Eor  direct  nitridation. 
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the  reaction  flux  is: 


r  —  kfCN^e  —  kbCcN^e  (97) 

The  units  of  the  reaction  flux  are  kmol  s“^,  which  gives  the  forward  reaction  rate 
in  units  of  kmol“^  s“^  for  consistency.  The  net  flux  for  the  surface  participating 
reaction  is  given  by: 

r  =  kf^N  —  kbCcN^e  (98) 

For  consistency,  the  units  of  the  forward  reaction  rate  are  s“^. 

The  direct  nitridation  equilibrium  constant  can  be  written  as: 


^N^C{b)^e  ^N^C{b) 

which  can  be  computed  using  the  CEA  database  [70]  and  the  minimization  of  Gibbs 
free  energy  (Equation  31).  As  before,  CEA  computes  the  activity-based  equilibrium 
constant  which  needs  to  be  converted  to  the  concentration-based  equilibrium  constant. 
The  surface  participating  equilibrium  constant  is  given  by: 


^  _  XcN^e 
^NXc{b) 

which  is  a  function  of  the  nitrogen  adsorption  reaction.  Substituting  in  Equation  79 
for  ^e/^N,  the  equilibrium  constant  can  now  be  specified  as: 


XNXc{h)  K2 

The  solution  process  is  the  same  as  described  previously,  but  only  the  flux  for  the 
surface  participating  nitridation  reaction  will  affect  the  surface  coverage  equation  for 
<h7v.  Because  the  ZA  model  does  not  include  any  other  CN  production  mechanism. 
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the  CN  mass  production  rate  is  simply: 


^CN  —  r  Mqn  (102) 

The  total  surface  mass  blowing  rate  will  now  also  include  contribution  from  this 
reaction  flux  computed  by  either  nitridation  mechanism.  Therefore,  to  close  the 
system  of  equations,  the  forward  reaction  rates  must  be  dehned  for  each  path. 

3.4.1  Determination  of  Forward  Rates. 

Even  though  the  purpose  of  Zhang  et  al.  [135]  was  to  publish  carbon  nitridation 
reaction  efhciencies,  the  most  fundamental  information  contained  in  the  article  is  the 
carbon  mass  loss  rate  at  different  surface  temperatures,  pressures,  and  nitrogen  con¬ 
centrations.  The  experiments  were  performed  in  N-N2  mixtures  with  carbon  mass  loss 
rates  measured  over  a  pure  graphite  rod.  The  reaction  efficiencies  simply  parametrize 
the  observed  mass  loss  within  a  very  simple  model  that  only  considers  the  reaction 
N  +  C(6)  — CN,  with  no  competing  surface  processes  like  adsorption  or  catalysis  nor 
the  limitation  of  a  hnite  number  of  active  sites.  The  tabulated  carbon  mass  loss  rates, 
which  are  averaged  over  time  and  surface  location,  can  be  matched  to  an  analytical 
expression  developed  for  a  simple  N-N2  system.  The  experimental  carbon  mass  loss 
rate  due  to  nitridation  can  be  computed  from  values  in  Tables  2  and  4  from  Zhang 
et  al.  [135]: 

_  ^^^C,con 

^  AtndsLs 

where  Amc,tot  is  the  total  carbon  mass  loss,  Amc,con  is  the  control  mass  loss.  At 
is  the  test  time,  dg  is  the  diameter  of  the  graphite  rod,  and  Lg  is  the  rod  length. 
The  diameter  of  all  the  rods  was  3.175  mm.  The  corresponding  test  conditions  are 
contained  in  Tables  1  and  3  from  Zhang  et  al.  [135] 
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The  experimental  carbon  mass  loss  rate  should  match  Equation  102,  but  with 
Me  instead  of  Mcn,  and  with  the  reaction  flux  specihed  by  either  Equation  97  or 
98,  depending  on  the  path.  The  fraction  of  available  surface  sites  covered  by  atomic 
nitrogen  must  be  considered,  which  may  be  different  under  different  environmental 
conditions,  because  it  will  influence  the  nitridation  rate  for  both  paths.  The  goal  is  to 
obtain  a  set  of  experimentally-based  kf  values  as  a  function  of  surface  temperatures 
that  reproduce  the  set  of  carbon  mass  loss  measurements.  Then  a  ht  can  be  computed 
based  on  an  ER  or  Arrhenius  form,  which  can  be  implemented  into  the  ZA  model. 

Consider  a  simple  N-N2  system  that  accounts  for  adsorption  and  ER  recombination 
at  the  surface: 


N  +  (s)  ^N(s) 

(104) 

N2  +  (s)  ^  N  +  N(s) 

(106) 

with  parameters  taken  from  Table  8.  The  concentration  of  empty  and  occupied  sites 
are  related  to  the  total  site  concentration  as: 


(|)  =  (|)g  +  $^  (106) 

where  <l>  =  5.8  x  10“®  kmol  m“^  [136].  The  next  subsections  develop  the  analytical 
expressions  for  the  forward  reaction  rates  based  on  either  nitridation  path. 

3.4.2  Direct  Nitridation. 

The  hrst  step  in  hnding  an  expression  for  fc/  is  to  determine  the  concentration  of 
empty  sites  from  the  concentration  of  sites  occupied  by  atomic  nitrogen.  The  rate  of 
change  of  the  adsorbed  atomic  nitrogen  under  steady-state  conditions  is: 
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d$7V 

dt 


0  =  n  +  r2 


(107) 


where  the  subscripts  1  and  2  refer  to  the  adsorption  and  ER  reactions,  respectively, 
shown  in  Equations  104  and  105.  Substituting  for  the  reaction  fluxes  and  empty 
surface  site  concentrations  gives: 


kflCjsi  (<h  —  ^n)  —  +  kf2Cjq^  (<l>  —  <l>Ar)  —  kb2C]y^N  —  0  (108) 


Solving  for  gives: 


= 


{kfiCN  +  kf2CN2)  d* 
kfiCN  +  kbi  +  kf2CN2  +  kb2CN 


(109) 


All  forward  and  backward  reaction  rates  are  taken  from  the  ZA  model,  and  the 
concentrations  of  N  and  N2  are  measured  in  the  experiment.  The  concentration  of 
empty  sites  can  then  be  calculated  by  <he  =  d*  —  d’Ar-  The  experimental  carbon  mass 
loss  rate  should  match  the  analytic  rate: 


The  —  r^Mc  —  {kf^CN^e  —  kbsCcN^e)  Me  (110) 


where  the  subscript  3  refers  to  the  direct  nitridation  reaction.  The  reaction  flux  for 
direct  nitridation  is  given  in  Equation  97.  Solving  for  kf^\ 


kf3  = 


kne  kbsCcN 

^cCnMc  Cn 


(111) 
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The  equilibrium  constant  for  the  nitridation  reaction  could  also  be  substituted  giving 
kfs  as: 

,  _  _ me _ 

^ 

The  concentration  of  CN  at  the  surface  was  not  measured  in  the  experiment,  so  it  is 
an  unknown.  However,  an  assumption  may  be  made  that  Ccn  is  negligible  compared 
to  Cn  and  0^2  {Ccn  -C  Cn,Cn2)  ,  which  gives  fc/3  as: 


^  me 


(113) 


In  fact,  this  assumption  is  consistent  with  the  experiments  because  the  mole  fractions 
were  computed  based  on  assuming  that  N  and  N2  were  the  only  gas  species  in  the  flow. 
Equation  113  provides  a  simple  relation  to  calculate  forward  reaction  rates  based  on 
experimental  carbon  mass  loss  rates  at  different  surface  temperatures. 

The  ER  forward  rate  can  be  used  to  dehne  the  direct  nitridation  reaction  as  shown 
in  Table  7.  The  required  inputs  for  the  ER  rate  are  jer,  P,  and  E^r-  Eor  a  pure  kinetic 
formulation,  P  equals  zero  and  7er  is  within  the  range  of  zero  to  one  [67].  The  reaction 
efficiency  is  dimensionless,  so  the  bracketed  term  provides  the  units  for  kf,  which  is 
consistent  with  the  net  reaction  flux  shown  in  Equation  97. 

The  ER  coefficients  can  be  determined  by  a  linear  least  squares  ht  on  the  log  of 
the  forward  reaction  rate.  The  resulting  linear  equation  is: 


In  (1^/1)  =  In 


+  In  (l7er|)  + 


{-Eer) 

RT 


(114) 


Equation  114  can  be  represented  by  the  matrix  equation  Ax  =  b  with  a  known  linear 
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least  squares  solution  of  x 


^  A^b  [4],  The  system 


in  matrix  form  is: 


1  ■ 

^  ~BT 

ln(|7er|) 

todl^/l)  '"(44) 

{  —  Eer) 

1  0 

0 

(116) 


where  the  constraint  0  <  7er  <  1  is  also  folded  into  the  system  of  equations.  The 
method  described  here  will  provide  an  ER  forward  reaction  rate  that  can  be  imple¬ 
mented  into  the  ZA  model. 


3.4.3  Surface  Participating  Reaction. 

For  the  surface  participating  reaction,  the  situation  is  more  complicated  because 
the  surface  coverage  is  coupled  to  the  nitridation  rate.  Now  the  rate  of  change  of  the 
adsorbed  atomic  nitrogen  under  steady-state  conditions  is: 

-— =  0  =  ri -Fra -rg  (116) 

where  ri  and  ra  represent  the  same  reactions  as  dehned  for  direction  nitridation  and 
r3  is  dehned  by  Equation  98.  Substituting  for  the  reaction  huxes  and  empty  surface 
site  concentrations  results  in: 


kflCN  —  ^N)  —  kbl^N  +  kf2CN2  ~  ^N)  —  kb2CN^N  —  kf3^N  +  kb3CcN  ($  “  ^n)  —  0 

(117) 


Solving  for  gives: 


— 


(/C/iC'tv  +  kf2C]Sl2  +  kb3CcN)  ® 
kfiCjq  -F  -F  kf2Ciq^  +  kb2C]\[  +  kf3  -F  kb3CcN 


(118) 
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The  analytic  solution  for  now  contains  forward  and  backward  rates  for  the  surface 
participating  nitridation  reaction.  Solving  for  the  experimental  carbon  mass  loss  rate 
shows: 

rhc  =  —  kbaCcN^e)  Me  (119) 

Substituting  for  and  using  the  known  equilibrium  constant  instead  of  kbs  gives: 


me  =  kfsMe^ 


kfiCjsf  +  kf2CN2  +  kbsCeN 


kfiC^  +  kbi  +  kf2Cpf2  +  kb2CM  +  fc/3  +  kbsCeN 


a 


CN 


1  - 


kfiCpf  +  kf2Cjq^  +  kbzCeN 


kfiC^  +  kbi  +  kf2CM2  +  kb2CN  +  kf3  +  kb^CeN 


(120) 


Assuming  again  that  CeN  CnjCn^^  carbon  mass  loss  rate  can  be  simplihed  to: 


kne  =  kfsMe^ 


kf^Civ  +  kfo.C]' 


(121) 


Solving  the  above  equation  for  kf^  gives: 


me 

,  _  Mc<h 

fc/3  -  - 


{kfiC^i  +  kbi  +  kf2CM2  +  kb2CN) 


kflCN  +  kf2CN2  — 


me 


(122) 


Again,  all  parameters  above  are  known  from  either  the  ZA  model  or  experiment. 
Accordingly,  the  forward  reaction  rate  values  can  be  tabulated  as  a  function  of  surface 
temperature. 

An  Arrhenius  expression  will  be  used  for  the  forward  rate  of  the  surface  partic¬ 
ipating  reaction  because  the  units  are  in  s“^.  The  coefficients  A,  /3,  and  Ea  can  be 
found  using  a  linear  least  squares  ht  as  was  done  previously.  For  consistency  with  the 
ZA  oxidation  reactions  that  are  speeffied  in  an  Arrhenius  form,  (3  is  set  equal  to  one. 
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Figure  18.  Forward  rate  fits  for  both  considered  carbon  nitridation  reactions. 

The  linear  equation  is: 


ln(|A:/|)  =  ln(|A|)  +  ln 


,  i-Eg) 
RT 


(123) 


The  linear  least  squares  solution  can  be  found,  and  the  surface  participating  nitrida¬ 
tion  reaction  can  also  be  added  to  the  ZA  model. 


3.4.4  Forward  Rate  Fits. 

The  results  of  the  linear  least  squares  solution  for  both  nitridation  paths  are  shown 
in  Table  10.  As  expected,  there  is  excellent  agreement  between  the  proposed  forward 
reaction  rates  and  experimental  data  as  shown  in  Figure  18.  Additional  experimental 
data  is  desired  to  further  validate  the  forward  rate  expressions,  especially  at  higher 
surface  temperatures.  However,  there  is  no  known  experimental  data  that  attempts 
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Table  10.  Forward  rate  fits  for  nitridation  reactions. 


Reaction 

Type 

7er  or  A  /3 

E  (k.J  mol 

N  +  (s)  +  C(6)  ^CN  +  (s) 

ER 

0.36  0 

36.86 

N(s)  +  C(6)  ^  CN  +  (s) 

Arrh 

0.57  1 

69.46 

to  characterize  the  nitrogen-carbon  snrface  interaction,  mnch  less  consider  the  steps 
of  the  snrface  reaction  process. 

The  direct  nitridation  reaction  has  a  lower  activation  energy  compared  to  the 
snrface  participating  reaction,  which  will  likely  give  the  occnrrence  of  this  reaction  a 
higher  probability.  Fnrthermore,  direct  nitridation  is  a  one-step  reaction  that  is  only 
dependent  on  the  availability  of  open  snrface  sites.  In  comparison,  the  snrface  par¬ 
ticipating  reaction  follows  a  general  snrface  reaction  process  similar  to  the  oxidation 
mechanisms  fonnd  in  the  ZA  model,  where  atomic  nitrogen  needs  to  be  in  an  adsorbed 
state  before  it  will  react  with  the  bnlk  carbon.  Ma  et  al.  [61]  demonstrated  adsorption 
and  diffnsion  of  N-atoms  on  single  layers  of  graphene  nsing  density-fnnctional  theory. 
Therefore,  the  reaction  that  represents  the  physical  process  may  be  the  Arrhenins- 
type  mechanism,  bnt  this  needs  to  be  experimentally  verihed.  It  is  expected  that  if 
both  reaction  types  predicted  CN  prodnction  at  the  snrface,  the  direct  nitridation 
mechanism  wonld  prodnce  larger  qnantities  of  CN. 

3.5  Radiative  Transitions  of  Cyanogen  (CN) 

Atoms  and  molecnles  emit  radiation  in  high  temperatnre  hypersonic  flows  that 
can  be  measnred  to  determine  radiative  heating  contribntions.  Radiative  transitions 
occnr  when  energy  is  released  resnlting  from  the  internal  energy  state  of  an  atom 
or  molecnle  going  from  an  npper  to  a  lower  level.  The  transitions  from  the  hrst 
two  excited  electronic  states  to  the  gronnd  state  of  CN  are  known  as  the  CN  red 
(A^n  — >  and  CN  violet  bands.  The  notation  nsed  to  indicate 
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Table  11.  Radiative  transition  wavelengths  for  the  CN  violet  band.  The  v’  and  v” 
colnmns  represent  different  vibrational  qnantnm  states  [52]. 


v’  v”  A  (nm) 

0  0  387.629 

0  1  420.956 

0  2  459.995 

1  0  358.145 

1  1  386.410 

1  2  419.057 

1  3  457.177 

2  0  333.274 

2  1  357.617 

2  2  385.404 

2  3  417.414 

2  4  454.677 

3  1  333.291 

3  2  357.299 

3  3  384.645 


v’  v”  A  (nm) 

3  4  416.067 

3  5  452.537 

4  0  293.755 

4  2  333.518 

4  3  357.224 

4  4  384.168 

4  5  415.053 

4  6  450.800 

5  3  333.985 

5  4  357.423 

5  5  384.009 

5  6  414.412 

5  7  449.508 

6  2  295.608 

6  3  314.082 


v’  v”  A  (nm) 

6  4  334.724 

6  5  357.930 

6  6  384.203 

6  7  414.184 

6  8  448.705 

7  3  296.885 

7  4  315.262 

7  5  335.765 

7  6  358.780 

7  7  384.791 

7  8  414.411 

7  9  448.439 

8  4  298.433 

8  5  316.743 

8  7  360.012 


v’  v”  A  (nm) 


8 

8 

385.813 

8 

9 

415.141 

8 

10 

448.759 

9 

4 

283.778 

9 

5 

300.284 

9 

6 

318.560 

9 

8 

361.667 

9 

9 

387.317 

9 

10 

416.422 

10 

5 

285.947 

10 

6 

302.472 

10 

7 

320.750 

10 

9 

363.791 

10 

10 

389.351 

the  type  of  transition  (i.e.  ^  X^S+)  comes  from  the  description  of  molecnlar 

strnctnre  nsing  qnantnm  physics.  The  different  electronic  states  are  represented  by 
X,  A,  B,  etc,  with  X  designating  the  gronnd  electronic  state.  The  hrst  three  total 
angnlar  momentnm  qnantnm  (energy)  states  are  designated  as  S,  11,  or  A.  The  pre- 
snperscript  and  snperscript  to  the  total  angnlar  momentnm  qnantnm  state  identihes 
the  total  spin  mnltiplicity  of  all  the  electrons  in  the  ontermost  orbit  and  the  symmetry 
of  the  electronic  wave  fnnctions,  respectively.  The  radiative  transition  wavelengths 
for  the  CN  violet  bands  are  given  in  Table  11  [52]. 


3.6  Radiation  Solver 


The  radiation  solver  nsed  is  the  Noneqnilibrinm  Air  Radiation  (NEQAIR)  code 
version  13.2  [125].  NEQAIR  is  a  spectral  high-resolntion  line-by-line  code  that  com- 
pntes  the  radiative  emission  and  adsorption  of  atomic  species  and  molecnlar  electronic 
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and  infrared  band  systems  along  a  line-of-sight.  The  solver  reqnires  temperatnres 
and  species  nnmber  densities  along  that  line-of-sight  to  calcnlate  qnantities  snch  as 
spectral  radiance  or  radiative  heating.  The  code  can  model  bonnd-free  and  free-free 
continnnm  radiation.  NEQAIR  has  been  nsed  extensively  to  compare  against  experi¬ 
mental  measnrements  from  a  nnmber  of  different  facilities  and  simnlated  atmospheres 
[21,  87,  68,  39].  As  noted  previonsly,  the  radiation  solver  is  rnn  nnconpled  from  the 
CFD  solver. 

The  popnlations  of  the  excited  electronic  states  can  be  fonnd  by  nsing  either  Boltz¬ 
mann  or  non-Boltzmann  distribntions.  The  Boltzmann  distribntion  is  appropriate 
nnder  eqnilibrinm  conditions  and  is  evalnated  at  the  specihed  electronic  temperatnre, 
which  is  eqnal  to  the  vibrational  temperatnre  nnder  the  assnmed  two-temperatnre 
model.  The  non-Boltzmann  distribntion  is  appropriate  for  thermal  noneqnilibrinm 
conditions  and  follows  the  qnasi-steady-state  (QSS)  method  [92,  125].  The  QSS 
method  assnmes  the  rate  of  change  in  the  popnlation  of  any  electronic  energy  level  is 
slow  compared  to  the  excitation/de-excitation  rates  into  and  ont  of  that  energy  level. 
Under  this  condition,  the  electronic  state  popnlations  are  determined  by  solving  a 
simple  set  of  linear  algebraic  eqnations  instead  of  a  set  of  differential  eqnations. 

All  NEQAIR  resnlts  nsed  a  non-Boltzmann  distribntion  for  the  popnlation  of 
excited  electronic  states.  Previons  stndies  [19,  130,  87,  71]  have  shown  that  Boltz¬ 
mann  distribntions  over  predict  radiation  intensities.  However,  these  over  predictions 
were  fonnd  in  highly  noneqnilibrinm  regimes,  and  in  particnlar,  for  carbon  containing 
environments  like  Mars  and  Titan.  The  CN  violet  bands  were  the  only  atomic  or 
molecnlar  systems  considered  for  the  NEQAIR  calcnlations  becanse  these  were  the 
primary  radiative  transitions  measnred  in  the  experiments.  The  radiative  transition 
wavelengths  for  the  CN  violet  bands  span  the  spectral  range  from  approximately  283- 
460  nm  [52] .  The  strongest  bands  are  contained  within  an  approximate  spectral  range 


81 


of  350-390  nm,  which  were  targeted  by  the  spectrometer  during  experimental  testing 
and  used  as  inputs  into  NEQAIR.  The  only  other  relevant  species  that  radiates  near 
this  wavelength  range  is  N^,  which  has  a  transition  wavelength  at  391.1  nm  [52], 
However,  Figure  10  from  Chapter  II  showed  there  were  no  substantial  spectra  mea¬ 
sured  near  this  wavelength  throughout  the  shock  layer  for  each  surface  temperature. 
Therefore,  was  not  contributing  to  the  measured  radiation  for  these  conditions. 

Radiance  comparisons  are  made  between  the  spectrometer  measurements  and  the 
NEQAIR  results.  The  radiance  is  dehned  as  the  radiative  flux  emitted  per  unit 
of  observed  area  and  is  calculated  by  integrating  the  spectral  radiance  between  the 
wavelength  range  of  353-391  nm.  The  radiance  is  measured  in  units  of  W/(cm^-sr). 
The  post-shock  flow  quantities  computed  by  US3D  are  extracted  from  the  volumetric 
solution  corresponding  to  the  experimental  line-of-sight.  The  experimental  line-of- 
sight  is  normal  to  the  stagnation  streamline  and  captures  the  radiation  emitted  from 
the  surface  out  to  the  shock.  Therefore,  multiple  slice  extractions  were  performed 
along  the  stagnation  streamline  that  provided  radiance  predictions  as  a  function  of 
wall  normal  distance.  The  length  of  the  line-of-sight  used  for  the  NEQAIR  simulations 
corresponded  to  the  X-2  diameter  of  8.5  cm. 

The  wavelength  calibration  of  the  UV  spectrometer  tended  to  drift  slightly,  so  the 
experimental  spectra  were  shifted  to  match  the  computed  CN  transition  at  approxi¬ 
mately  388.4  nm.  The  slit  function  of  the  experimental  setup  has  been  characterized 
as  a  Voigt  prohle  with  Gaussian  and  Lorentzian  line  widths  of  1.1  nm  and  0.3  nm, 
respectively.  To  verify  this  slit  function.  Figure  19  plots  the  scaled  CN  spectral  ra¬ 
diance  0.2  mm  away  from  the  wall  along  the  stagnation  streamline  for  Ty;  =  2410  K. 
The  plot  displays  the  results  from  the  experiment  (X-2)  and  from  using  the  ZA  and 
Park  ablation  models.  Each  data  set  was  scaled  by  its  maximum  spectral  radiance 
value.  After  convolution  with  the  Voigt  prohle,  the  numerical  spectra  calculated  by 
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Figure  19.  Comparison  of  the  scaled  CN  spectral  radiance  0.2  mm  away  from  the  wall 
along  the  stagnation  streamline  for  the  =  2410  K  case. 

NEQAIR  matched  the  experimental  spectra  well.  The  same  Gaussian  and  Lorentzian 
line  widths  that  comprised  the  Voigt  prohle  were  used  for  all  solutions. 


3.6.1  Spatial  Resolution  Function. 

A  direct  comparison  of  the  radiance  measured  in  the  experiment  and  that  pre¬ 
dicted  by  NEQAIR  has  limited  accuracy  because  there  are  instrument  constraints. 
The  spatial  resolution  of  the  collection  and  spectrometer  optics  is  limited  by  how 
sharply  the  shock  layer  can  be  resolved.  The  resolution  of  the  image  focused  on  the 
spectrometer  camera  is  determined  by  the  quality  of  optical  alignment.  Also,  the 
spectrometer  camera  is  subject  to  charge  smearing,  which  causes  adjacent  pixels  to 
share  intensity  measurements.  The  net  effect  of  these  resolution  limitations  results  in 
a  spatial  prohle  that  is  broadened  in  comparison  to  the  actual  result. 

Eollowing  the  work  of  Brandis  et  al.  [20],  the  broadening  may  be  accounted  for 
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by  convolving  the  predicted  radiance  compnted  by  NEQAIR  with  a  spatial  resolntion 
function  (SRF).  The  SRF  is  a  convolntion  of  the  optical  and  camera  resolntion  fnnc- 
tions  that  is  determined  by  the  specihc  experimental  setnp.  The  convolntion  fnnction 
nsed  for  this  work  was  based  on  similar  experiments  at  NASA  FAST,  which  matched 
the  same  UV  spectrometer  camera  settings.  The  camera  pixel  to  wall  normal  distance 
correlation  was  mapped  to  the  X-2  setnp  to  calcnlate  appropriate  convolntions.  Note 
that  this  is  an  estimated  SRF  calcnlation,  bnt  it  shonld  provide  a  better  comparison 
to  the  experimental  resnlts. 

3.7  ZA  Rate  Coefficients  Screening 

Compntational  analyses  are  sensitive  to  the  physical,  chemical,  and  nnmerical 
models  nsed  for  aeroheating  predictions  [32],  The  inpnt  parameters  of  the  models 
generally  nse  a  single  nominal  valne,  bnt  there  is  a  range  of  nncertainty  associated  with 
each  valne  [88],  To  qnantify  this  nncertainty,  Monte  Carlo  methods  are  widely  nsed 
dne  to  their  robnstness  and  convergence  rate  properties.  A  Monte  Carlo  nncertainty 
analysis  provides  statistical  data  that  determines  the  fractional  contribntion  of  each 
inpnt  parameter  to  the  overall  variability  of  a  selected  ontpnt  variable  [88].  However 
to  reach  statistical  convergence,  2000-3000  CFD  solntions  are  reqnired  for  each  inpnt 
parameter,  which  makes  the  method  very  compntationally  expensive. 

A  more  recently  developed  method,  called  a  dehnitive  screening  design,  has  the 
ability  to  address  the  sensitivity  of  inpnt  variables  on  qnantities  of  interest.  Jones 
and  Nachtsheim  [47]  proposed  a  three-level  dehnitive  screening  design  that  assesses 
the  relative  impact  of  a  large  nnmber  of  factors  with  the  ability  to  identify  main 
(linear)  and  second-order  (non-linear)  effects  independent  of  two-factor  interactions. 
The  design  nses  three  factor  levels  with  the  nnmber  of  rnns  eqnal  to  2m  -|-  1  for  m 
factors  with  one  center  point  rnn  (baseline).  The  general  strnctnre  for  m  factors  is 
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Table  12.  General  structure  of  a  definitive  screen  design  with  m  factors  [47] 


Run  Factor  Levels 


(*) 

Xi^l 

Xi^2 

Xi,^ 

•  •  • 

1 

0 

±1 

±1 

...  ±1 

2 

0 

±1 

±1 

...  ±1 

3 

±1 

0 

±1 

...  ±1 

4 

±1 

0 

±1 

...  ±1 

5 

±1 

±1 

0 

...  ±1 

6 

±1 

±1 

0 

...  ±1 

2m  —  1 

±1 

±1 

±1 

0 

2m 

±1 

±1 

±1 

0 

2m  ±  1 

0 

0 

0 

0 

shown  in  Table  12.  Each  run  has  exactly  one  factor  at  its  center  point  and  sets  all 
other  factors  to  high  or  low  levels.  In  addition,  two  runs  were  added  where  no  factors 
were  set  at  their  center  point. 

The  13  pre-exponential  reaction  rate  coefEcients  of  the  ZA  model  (with  nitrida- 
tion)  were  screened  for  their  relative  signihcance  in  determining  the  total  CN  radia¬ 
tive  heat  flux.  In  the  absence  of  precise  uncertainty  estimates  for  each  coefficient, 
the  screening  process  was  performed  by  varying  the  values  within  a  conservative  un¬ 
certainty  range  guided  by  the  work  of  Palmer  et  al.  [88]  The  adsorption  sticking 
coefficients,  Sq,  and  ER  reaction  efficiencies,  7er,  were  varied  by  ±10%,  while  the  Ar¬ 
rhenius  reaction  rate  coefficients  were  varied  by  ±1  order  of  magnitude.  The  total  CN 
radiance  from  each  combination  of  factors  was  compared  to  the  case  with  all  coeffi¬ 
cients  set  to  their  nominal  values  (i.e.  baseline  case).  The  quantitative  metric  was  the 
absolute  percent  difference  from  the  baseline  case.  The  goal  of  the  screening  process 
is  to  show  which  surface  reactions  have  the  greatest  influence  on  the  overall  variabil¬ 
ity  of  the  CN  radiance.  Specifically,  identification  of  those  reactions  that  increase 
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the  radiance  were  of  most  interest.  It  shonld  be  evident  that  an  efficiency  greater 
than  one  is  not  physically  meaningfnl  in  the  context  of  the  snrface  reaction  model 
framework.  However  for  the  pnrpose  of  this  screening  design,  it  will  provide  insight 
into  the  sensitivity  of  the  efficiency  parameters.  Before  the  Park  and  ZA  models  are 
applied  to  the  X-2  experiments,  the  next  chapter  will  characterize  the  performance 
of  each  model  nnder  representative  hypersonic  re-entry  conditions. 
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IV.  Characterization  of  Ablation  Models 


In  order  to  better  understand  the  behavior  of  the  different  ablation  models,  a 
simple  set  of  test  cases  is  used  that  is  representative  of  a  typical  re-entry  trajectory. 
The  ablation  models  will  be  applied  to  the  simulations  described  by  Candler  [23],  The 
results  presented  below  examine  the  influence  of  nitridation,  the  effect  of  immobile 
or  mobile  adsorption,  and  the  sensitivity  of  the  total  surface  site  concentration. 

4.1  Simulation  Parameters 

The  test  cases  involved  a  10  cm  radius  sphere,  8°  cone  geometry  with  freestream 
conditions  corresponding  to  altitudes  of  20,  25,  30,  and  40  km  at  a  velocity  of  7  km 
sec“^.  The  surface  temperatures  and  pressures  for  each  altitude  are  shown  in  Fig¬ 
ure  20.  The  flow  fields  were  computed  on  the  same  grid  used  in  Candler  [23],  which 
consisted  of  48,000  hexahedral  elements  (160  axial  and  300  wall-normal).  The  gas 
is  considered  to  be  composed  of  the  following  species:  N2,  O2,  NO,  CO2,  CO,  C2, 
C3,  CN,  C,  N,  and  O.  Candler’s  results  compared  different  gas-phase  models  cou¬ 
pled  with  the  ZA  ablation  boundary  condition.  The  results  presented  here  used  a 
gas-phase  chemistry  model  essentially  equivalent  to  that  shown  in  Table  4,  but  with 
the  elimination  of  charged  species  chemistry  and  the  addition  of  important  reactions 
identified  by  Martin  and  Boyd  [68].  The  three  additional  Martin  and  Boyd  reactions 
shown  in  Arrhenius  form  (Equation  27)  are  given  in  Table  13. 

4.2  ZA  Model  Verification 

First,  a  code  verification  study  was  performed  because  the  implementation  of 
the  ZA  model  was  relatively  complex.  The  ZA  model  implemented  in  US3D  was 
verified  against  the  NASA  DPLR  code  [129,  64].  The  verification  used  the  30  km  test 
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(a)  Surface  temperature 


(b)  Surface  pressure 


Figure  20.  Surface  temperature  and  pressure  distributions  for  each  representative  re¬ 
entry  test  case. 
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Table  13.  Martin  and  Boyd  reactions  and  rate  coefficients  [68]. 


q 

Reaction 

Tf 

1 

C  +  NO  ^  CO  +  N 

2.3e+10 

0.00 

0.00 

T 

2 

C2  +  Na  ^  CN  +  CN 

1.5e+10 

0.00 

2.10e+4 

T 

3 

N  +  CO2  ^  NO  +  CO 

3.0e+9 

0.00 

5.69e+3 

T 

Fignre  21.  Comparison  of  species  snrface  mass  flnxes  predicted  by  the  US3D  and  DPLR 
codes  nsing  the  Zhlnktov  and  Abe  ablation  model. 

conditions.  Fignre  21  shows  the  wall  temperatnre  distribntion  along  with  the  species 
snrface  mass  flnxes  compnted  by  the  US3D  and  DPLR  codes.  There  is  excellent 
agreement  between  the  codes  and  verihes  the  implementation  of  the  ZA  model  in 
US3D. 


4.3  Influence  of  Nitridation 

The  following  hgnres  will  display  resnlts  for  the  fonr  different  cases  considered  at 
the  fonr  freestream  conditions.  The  original  ZA  model  withont  nitridation  is  iden- 
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tified  as  the  baseline  model.  The  modihed  ZA  model  includes  the  direct  nitridation 


reaction  and  the  surface  participating  nitridation  reaction  which  are  identihed  as  “ER 
nitridation”  and  “Arrh  nitridation”,  respectively,  due  to  the  reaction  types. 

4.3.1  Surface  Mass  Flux  Comparison. 

First,  the  total  mass  flux  from  the  surface  is  compared  as  a  function  of  surface  dis¬ 
tance  measured  from  the  stagnation  point.  Figure  22  plots  this  quantity  and  shows 
that  at  strong  sublimation  conditions  (20  km),  all  results  gave  similar  total  mass 
fluxes  near  the  stagnation  point,  with  the  ZA  model  (with  or  without  nitridation) 
being  about  10%  higher.  Further  along  the  surface,  as  the  wall  temperature  cools, 
the  effect  of  direct  nitridation  is  seen  where  there  is  an  increase  in  mass  loss  rate  over 
a  short  surface  distance  compared  to  the  baseline  model.  The  effect  of  direct  nitrida¬ 
tion  becomes  more  signihcant  as  the  wall  temperatures  and  pressures  are  lowered  at 
the  higher  altitude  conditions  and  as  sublimation  becomes  more  of  an  inconsequential 
reaction.  The  largest  impact  of  direct  nitridation  occurred  at  the  30  km  condition 
where  it  increased  the  surface  mass  flux  at  the  stagnation  point  by  about  154%  com¬ 
pared  to  the  baseline  model.  In  comparison,  the  Park  model  had  an  approximate  46% 
higher  mass  flux  than  the  baseline  model.  Even  though  there  are  substantial  increases 
to  the  mass  fluxes  caused  by  direct  nitridation,  it  does  not  necessarily  translate  to 
signihcant  increases  in  gaseous  CN  concentrations.  For  most  cases,  the  modihed  ZA 
mass  loss  rates  relax  to  the  baseline  values  beyond  the  nose  of  the  sphere-cone  geom¬ 
etry.  Clearly,  the  surface  participating  nitridation  reaction  had  a  very  minimal  effect 
compared  to  the  baseline  model  and  did  not  alter  any  mass  loss  rates  or  gas-phase 
species  compositions  as  will  be  shown  in  subsequent  hgures. 

The  carbon-bearing  species  surface  mass  huxes  are  plotted  in  Figures  23-29  for 
each  test  condition.  At  the  20  km  condition,  the  baseline  and  modihed  ZA  model 
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Surface  Mass  Flux  (kg/m^  sec)  Surface  Mass  Flux  (kg/m^  sec) 


(a)  20  km  (b)  25  km 


(c)  30  km  (d)  40  km 

Figure  22.  Predicted  total  surface  mass  flux  as  a  function  of  normalized  distance  from 
stagnation  point. 
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shows  that  C3  formation  is  the  dominant  surface  reaction  at  the  stagnation  point, 
but  the  Park  model  has  relatively  equal  mass  flux  contributions  from  C3  sublimation 
and  nitridation.  Furthermore,  the  Park  model  produces  less  C3  than  the  ZA  models. 
The  ZA  models  have  C2  and  C  recombining  at  the  surface,  whereas  the  Park  model 
has  zero  mass  flux  for  both  of  these  species  at  the  stagnation  point.  However,  as 
noted  above,  the  Park  model  has  no  mechanism  for  C2  formation  or  C  recombination 
on  the  surface,  which  will  become  an  important  distinction  at  the  higher  altitude 
conditions.  The  Park  model  shows  signihcant  CN  production  around  the  entire  nose 
of  the  geometry,  but  the  modihed  ZA  model  with  direct  nitridation  only  has  CN 
forming  around  the  sphere-cone  juncture. 

For  the  25  km  case,  oxidation  and  nitridation  become  more  important  closer  to  the 
stagnation  point  while  the  level  of  sublimation  decreases  for  the  modihed  ZA  model. 
The  CN  formation  rate  has  increased  near  the  stagnation  region,  but  it  has  also  caused 
the  amount  of  oxidation  to  decrease  compared  to  the  baseline  model.  Sublimation  and 
nitridation  remain  the  dominant  reactions  for  the  Park  model.  Also,  the  Park  model 
predicts  approximately  zero  CO  production.  At  the  stagnation  point,  the  Park  model 
predicts  four  times  as  much  CN  formation  compared  to  the  modihed  ZA  model.  The 
level  of  nitridation  becomes  similar  between  the  Park  and  modihed  ZA  models  at  the 
30  km  condition,  which  is  shown  in  Figure  27.  Interestingly,  the  production  of  CN 
is  now  15%  higher  at  the  stagnation  point  for  the  modihed  ZA  model  compared  to 
the  Park  model.  The  amount  of  oxidation  predicted  by  the  modihed  ZA  model  near 
the  stagnation  point  is  also  substantially  reduced  compared  to  the  baseline  model. 
There  are  relatively  large  levels  of  C  recombination  at  the  wall  with  direct  nitridation, 
which  will  affect  the  CN  concentrations  in  the  near-wall  boundary  layer.  In  contrast, 
the  Park  model  shows  nitridation  as  the  only  relevant  surface  reaction  for  this  test 
condition. 
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Surface  Mass  Flux  (kg/m^  sec)  Surface  Mass  Flux  (kg/m^  sec) 


(a)  CO2 


(b)  CO 


(c)  C3  (d)  C2 

Figure  23.  CO2,  CO,  C3,  and  C2  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  stagnation  point  at  20  km  conditions. 


93 


(a)  C  (b)  CN 

Figure  24.  C  and  CN  surface  mass  fluxes  as  a  function  of  normalized  distance  from 
stagnation  point  at  20  km  conditions. 

The  competition  between  the  oxidation  and  nitridation  processes  is  further  high¬ 
lighted  in  Figure  29,  which  shows  the  40  km  species  surface  mass  fluxes.  Now,  CO2 
formation  is  the  dominant  surface  reaction  with  nitridation  still  contributing  to  the 
total  mass  flux  for  the  modified  ZA  model.  The  ZA  models  predict  CO  to  recombine 
at  the  surface,  which  is  opposed  to  the  Park  model  that  predicts  CO  production.  The 
CN  production  rate  for  the  Park  model  has  again  overtaken  the  modified  ZA  model 
being  60%  higher  at  the  stagnation  point.  Thus,  even  though  the  total  mass  flux 
at  the  stagnation  point  for  the  Park  and  modified  ZA  models  is  fairly  close  (within 
10%),  the  mechanisms  are  completely  different. 

The  previous  analysis  has  shown  that  nitridation  can  be  as  or  more  important 
of  a  surface  reaction  as  oxidation  and  sublimation  under  certain  simulated  re-entry 
conditions.  The  implemented  direct  nitridation  mechanism  proved  to  have  a  more 
significant  impact  on  surface  mass  fluxes  than  the  surface  participating  nitridation 
mechanism.  The  reason  is  due  to  a  lack  of  N  adsorbing  to  the  surface,  which  is  a 
precursor  step  for  the  nitridation  reaction.  For  all  of  the  test  conditions,  the  mass  flux 
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Figure  25.  CO2,  CO,  C3,  and  C2  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  stagnation  point  at  25  km  conditions. 
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Figure  26.  C  and  CN  surface  mass  fluxes  as  a  function  of  normalized  distance  from 
stagnation  point  at  25  km  conditions. 


of  N  either  away  or  towards  the  surface  was  essentially  zero.  The  effect  of  nitridation 
can  be  further  assessed  by  considering  Figures  30-33,  which  plot  select  gas-phase 
species  concentrations  in  the  boundary  layer.  The  hgures  will  focus  on  how  the  CN 
mass  fraction  changes  between  the  different  models. 


4.3.2  Species  Mass  Fraction  Comparison. 

Figure  30  plots  the  species  mass  fractions  as  a  function  of  distance  from  the  surface 
at  the  stagnation  point  for  the  20  km  case.  Note  that  the  same  gas-phase  chemistry 
model  was  used  with  each  surface  reaction  model.  The  initial  mass  fractions  at 
the  surface  reflect  the  surface  mass  flux  results.  Therefore,  C2  and  CN  have  higher 
fractions  at  the  wall  for  the  Park  model  and  C3  is  higher  at  the  wall  for  the  ZA  models. 
C2  and  CN  form  from  C3  immediately  away  from  the  surface.  Then,  as  the  amount  of 
C3  continually  diminishes,  the  concentrations  of  C2  and  CN  also  rapidly  decline.  The 
ZA  and  Park  models  predict  similar  CN  mass  fraction  prohles  due  to  the  differences 
in  C3  and  CN  production  rates.  Even  though  the  amount  of  CN  formation  at  the 
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27.  Species  surface  mass  fluxes  as  a  function  of  normalized  distance  from  stag- 
point  at  30  km  conditions. 
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(a)  C  (b)  CN 

Figure  28.  C  and  CN  surface  mass  fluxes  as  a  function  of  normalized  distance  from 
stagnation  point  at  30  km  conditions. 

stagnation  point  is  mnch  lower  for  the  ZA  model  than  the  Park  model,  the  higher 
levels  of  C3  predicted  by  the  ZA  models  accelerate  CN  formation  in  the  gas-phase.  In 
Fignre  31  for  the  25  km  case,  the  CN  mass  fraction  prohles  are  dne  to  the  amonnt  of 
nitridation  that  occnrred  at  the  snrface  for  either  ablation  model  and  the  dissociation 
rate  of  CN  determined  from  the  gas-phase  chemistry.  Note  the  large  difference  in 
C  concentration  between  the  Park  and  ZA  models,  which  is  a  common  resnlt  at  all 
freestream  conditions  and  will  have  a  strong  inflnence  on  CN  concentrations  at  the 
higher  altitndes. 

At  the  two  highest  altitnde  conditions,  the  ability  of  C  to  recombine  at  the  snrface 
is  the  primary  discrepancy  between  the  resnlting  CN  mass  fraction  prohles.  Fignres  32 
and  33  show  the  species  mass  fractions  for  the  30  km  and  40  km  cases.  Again,  for 
these  cases  the  oxidation  and  nitridation  mechanisms  are  the  most  dominant  snrface 
reaction  processes.  It  is  interesting  to  observe  the  concentration  of  CN  is  higher  for 
the  Park  model  than  the  ZA  model  with  direct  nitridation  becanse  the  CN  mass  hnx 
at  the  stagnation  point  was  higher  for  the  ZA  model.  However,  this  can  be  explained 
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Figure  29.  Species  surface  mass  fluxes  as  a  function  of  normalized  distance  from  stag¬ 
nation  point  at  40  km  conditions.  At  these  conditions,  there  is  no  flux  of  C3  or  C2  for 
any  model  considered. 
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Figure  30.  Species  mass  fractions  as  a  function  of  surface-normal  distance  at  the  stag¬ 
nation  point  at  20  km  conditions. 
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31.  Species  mass  fractions  as  a  function  of  surface- normal  distance  at  the  stag- 
point  at  25  km  conditions. 
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32.  Species  mass  fractions  as  a  function  of  surface-normal  distance  at  the  stag- 
point  at  30  km  conditions. 
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Mass  Fraction  Mass  Fraction 


(a)  CO2  (b)  CO 


(c)  C  (d)  CN 

Figure  33.  Species  mass  fractions  as  a  function  of  surface-normal  distance  at  the  stag¬ 
nation  point  at  40  km  conditions. 
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by  comparing  the  C  mass  fraction.  The  Park  model  predicts  relatively  more  C  in 
the  boundary  layer  than  the  ZA  model,  which  in  turn  produces  more  CN  due  to  the 
gas-phase  reaction  N2  +  C  <->■  CN  +  N.  At  these  freestream  conditions,  only  a  small 
amount  of  N2  dissociates  across  the  shock,  so  it  is  readily  available  to  react  with 
C  in  the  boundary  layer.  A  similar  conclusion  can  be  drawn  for  the  results  at  the 
40  km  condition  shown  in  Figure  33.  The  Park  model  consistently  predicts  a  higher 
CN  mass  fraction  in  the  boundary  layer  than  the  modihed  ZA  models.  The  larger 
C  mass  fraction  is  again  a  main  contributor  to  this  difference,  but  additionally,  the 
CN  surface  mass  flux  at  the  stagnation  point  for  the  Park  model  is  also  higher  than 
the  direct  nitridation  ZA  model.  All  of  these  results  stem  from  the  species  that  are 
allowed  to  react  at  the  surface  and  the  preferred  direction  of  reactions,  dehned  in  the 
ablation  models. 

4.3.3  Nitridation  Comparison. 

The  surface  temperatures  and  pressures  at  which  nitridation  is  or  is  not  active 
for  the  direct  nitridation  ZA  model  are  shown  in  Figure  34  to  further  differentiate 
the  nitridation  mechanisms  from  the  Park  model.  By  comparing  all  four  plots,  it 
is  fair  to  say  that  the  nitridation  rate  has  a  stronger  temperature  dependence  than 
pressure  dependence.  Direct  nitridation  is  most  active  between  surface  temperatures 
of  approximately  2000-4000  K,  with  a  consistent  maximum  around  3600  K.  As  the 
temperature  approaches  and  goes  below  2000  K,  the  nitridation  rate  goes  to  zero.  The 
rate  also  appears  to  go  to  zero  as  the  temperatures  reach  4600  K.  In  contrast,  the 
Park  model  nitridation  rate  remains  relatively  constant  at  very  high  wall  temperatures 
above  4000  K.  The  rate  does  decrease  as  the  temperature  decreases  but  with  a  smaller 
slope  than  the  direct  nitridation  ZA  model.  Even  as  temperatures  decrease  below 
2000  K,  the  Park  model  predicts  a  CN  mass  flux  greater  than  zero,  whereas  the 
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direct  nitridation  rate  goes  to  zero. 


4.4  Surface  Adsorption  Process 

The  effect  of  assuming  the  surface  adsorption  process  as  mobile  or  immobile  is 
analyzed  for  each  test  condition.  The  baseline  ZA  model  assumes  mobile  adsorp¬ 
tion,  but  there  is  no  direct  experimental  evidence  that  it  dominates  at  all  conditions. 
Therefore,  the  limiting  case  of  immobile  adsorption  is  studied  to  understand  the  sensi¬ 
tivity  of  the  ablation  process  to  the  surface  mobility.  Figures  22-33  display  the  surface 
mass  fluxes  and  species  mass  fractions  under  the  assumption  of  immobile  adsorption 
at  each  test  condition. 

At  the  lower  altitude  conditions,  the  immobile  adsorption  solution  predicts  a  sim¬ 
ilar  surface  mass  flux  (total  and  species  specihc)  as  the  solution  predicted  with  a 
mobile  adsorption  assumption  near  the  nose  of  the  sphere-cone  geometry.  There  is  a 
signihcant  reduction  in  the  total  surface  mass  flux  for  the  higher  altitude  cases  with 
the  assumption  of  immobile  adsorption.  In  fact,  immobile  adsorption  reduces  the 
mass  loss  rate  to  almost  zero  over  the  entire  surface  for  the  40  km  case.  The  inability 
of  oxygen  and  nitrogen  atoms  to  diffuse  across  the  surface  reduced  reaction  probabil¬ 
ities  and,  when  combined  with  the  much  lower  surface  temperatures  and  pressures  of 
the  40  km  case,  generated  very  little  mass  loss. 

The  oxidation  reactions  are  affected  the  most  under  an  immobile  assumption  for 
all  test  conditions,  while  the  sublimation  reactions  are  altered  slightly  for  the  30  km 
test  condition.  When  comparing  the  mobile  and  immobile  CO2  and  CO  mass  fluxes, 
immobile  adsorption  has  the  tendency  to  increase  CO2  levels  and  causes  CO  to  recom¬ 
bine  at  the  surface.  Mobile  adsorption  predicts  the  opposite  and  shows  CO  formation 
at  similar  surface  locations  as  shown  in  Figures  23b,  25b,  and  27b.  The  CO2  and  CO 
mass  fractions  along  the  stagnation  streamline  also  highlight  these  differences. 
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(a)  20  km  (b)  25  km 


(c)  30  km  (d)  40  km 

Figure  34.  CN  surface  mass  fluxes  as  a  function  of  surface-normal  distance  at  the 
stagnation  point  for  each  test  case  using  the  modifled  Zhluktov  and  Abe  model  with  the 
Eley-Rideal  nitridation  mechanism.  Displaying  the  surface  temperatures  and  pressures 
at  speciflc  locations. 
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Figure  35.  Surface  pressure  distributions  for  representative  re-entry  cases. 

As  expected,  the  sublimation  dominated  test  conditions  at  20  and  25  km  have  no 
changes  to  any  of  the  predicted  species  mass  fractions  since  the  stagnation  point  mass 
loss  rates  are  identical  for  mobile  and  immobile  adsorption.  Since  immobile  adsorption 
predicts  strong  CO  recombination  at  the  stagnation  point  for  the  30  km  case,  the  CO 
mass  fraction  at  the  wall  is  reduced  by  73%  compared  to  the  mobile  adsorption 
prediction.  The  amount  of  CN  is  also  reduced  due  to  the  lower  concentrations  of  CO. 
For  the  40  km  case,  there  are  signihcantly  rednced  levels  of  both  CO2  and  CO  for 
immobile  adsorption,  and  the  snrface  is  essentially  a  non-reactive,  as  compared  to  the 
mobile  adsorption  predictions. 

The  conditions  at  which  mobile  and  immobile  adsorption  would  be  expected  to 
provide  similar  results  are  correlated  with  the  surface  pressure  and  possibly  also  with 
the  snrface  temperatnre.  Fignre  35  plots  the  snrface  pressnre  distribntions  for  each 
test  condition.  Comparing  this  plot  with  Fignres  22-29,  any  snrface  pressnre  above 
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approximately  20  atm  produces  similar  surface  reaction  processes  independent  of  a 
mobile  or  immobile  adsorption  assumption.  Conversely,  any  surface  pressure  below 
approximately  16  atm  produces  different  oxidation  mechanisms  dependent  on  a  mobile 
or  immobile  assumption.  A  transition  zone  occurs  in  between  these  pressures  where 
the  surface  temperature  may  have  an  influence  on  the  relevant  oxidation  processes. 
When  the  surface  pressure  is  in  this  transition  region  and  the  surface  temperature 
is  at  sublimation  conditions,  immobile  adsorption  tends  to  predict  higher  levels  of 
CO2  formation  and  CO  recombination  as  compared  to  mobile  adsorption.  However, 
when  the  surface  pressure  is  in  this  transition  region  and  the  surface  temperature  is 
near  oxygen  saturation  conditions,  the  predicted  species  mass  loss  rates  are  similar 
for  mobile  and  immobile  adsorption.  Hence,  the  differences  in  the  species  mass  fluxes 
are  correlated  with  the  surface  pressure  and  temperature,  but  the  total  surface  mass 
flux  appears  to  be  directly  correlated  with  the  surface  pressure,  as  seen  in  Figure  22. 

4.5  Surface  Site  Concentration  Sensitivity 

One  of  the  key  attributes  of  the  ZA  model  that  differentiates  it  from  other  ab¬ 
lation  models  is  the  specihcation  of  a  hnite  number  of  surface  reaction  sites.  The 
surface  reaction  site  concentration  is  also  called  the  “active  site  density”  and  is  usu¬ 
ally  dehned  with  the  parameter  $.  Different  grades  of  graphite  exhibit  large  property 
variations  because  of  different  impurity  levels,  grain  sizes,  grain  boundaries,  and  struc¬ 
ture  (amorphous  or  highly-oriented) ,  which  complicate  gas-surface  interaction  studies 
[134].  Determining  the  sensitivity  of  <h,  or  characterizing  an  ablative  surface  in  gen¬ 
eral,  is  essential  to  any  surface  reaction  study  because  it  could  have  a  direct  impact 
on  the  level  of  oxidation  predicted  by  the  ZA  model.  In  fact,  this  is  probably  the 
most  difficult  parameter  to  determine  for  any  system.  The  surface  site  concentration 
parameter  could  also  serve  as  a  means  to  control  the  reactivity  of  the  graphite  model. 
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The  active  site  density  of  3.5  x  10^®  atoms  (i.e.  =  5.8  x  10“®  kmol  m“^) 

published  by  Zhluktov  and  Abe  [136]  was  actually  taken  from  a  study  performed  by 
Blyholder  and  Eyring  [17].  Blyholder  and  Eyring  specihed  the  active  site  density  as 
the  number  of  carbon  atoms  per  cm^  of  graphite  lattice.  It  is  unknown  how  this  value 
was  determined,  but  it  can  be  related  to  the  physical  locations  where  oxygen  atoms 
tend  to  adsorb  on  a  graphene  sheet.  Numerous  studies  have  shown  that  adsorbed 
atomic  oxygen  mostly  forms  epoxide-like  structures  (C-O-C)  meaning  a  majority  of 
the  reaction  sites  are  located  between  C-C  bonds  [110,  11,  99]. 

As  a  simplihcation  to  this  observation,  several  assumptions  are  made.  Eirst,  as¬ 
sume  a  graphene  sheet  is  arranged  in  a  honeycomb  lattice  (hexagonal  pattern)  with 
a  C-C  atom  separation  distance  of  0.142  nm.  Eor  this  structure,  the  number  of  C-C 
bonds  is  equal  to  the  number  C-atoms  for  one  cell.  Therefore,  it  is  assumed  active 
sites  are  equivalent  to  C-atom  locations  to  take  advantage  of  known  geometrical  rela¬ 
tionships.  The  average  distance  between  active  sites  is  given  by  Marschall  et  al.  [67] 
as: 

A  =  — ^ 

If  the  average  distance  is  assumed  to  be  0.142  nm,  the  active  site  density  equals 
8.2  X  10“®  kmol  m“^.  If  the  Blyholder  and  Eyring  active  site  density  is  used,  the 
average  distance  is  equal  to  0.169  nm.  A  more  accurate  average  distance  between 
active  sites  is  found  by  calculating  the  distance  between  one  C-atom  and  the  other  12 
C-atom  nearest  neighbors  by  considering  a  three-cell  hexagonal  structure.  Using  this 
method,  the  average  distance  is  found  to  be  0.229  nm,  which  equates  to  an  active  site 
density  of  3.2  x  10“®  kmol  m“^.  Therefore,  it  has  been  shown  that  a  range  of  active 
site  densities  could  be  considered  depending  on  how  reaction  sites  are  dehned. 

A  surface  reaction  site  concentration  sensitivity  study  was  performed  using  the 
baseline  ZA  model  for  the  30  km  test  case.  Eive  surface  site  concentrations  were 
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Table  14.  Surface  site  concentration  sensitivity  study  parameters. 


(kmol  m 

B  (m 

Aeq  (m^  kmol 

1.0  X  10"® 

6.022  X  10^® 

1.0  X  10® 

3.0  X  10"® 

1.807  X  10^9 

3.33  X  10^ 

5.0  X  10"® 

3.011  X  10^9 

2.0  X  10^ 

7.0  X  10"® 

4.215  X  10^9 

1.43  X  10^ 

9.0  X  10"® 

5.420  X  10^9 

1.11  X  10^ 

considered  that  encompass  the  range  of  densities  derived  above  given  in  Table  14. 
The  atomic  oxygen  and  nitrogen  adsorption  eqnilibrinm  constant  coefficients  (Aeg) 
are  also  provided  in  Table  14  becanse  the  valnes  change  based  on  the  valne  of  <h. 

Fignres  36-39  show  the  snrface  mass  flnxes  and  mass  fractions  along  the  stagnation 
streamline.  The  ablation  predictions  are  insensitive  to  the  active  site  density.  Even 
with  wide  variations  in  snrface  temper atnre  and  pressnre,  the  predicted  snrface  mass 
flnxes  and  mass  fractions  are  essentially  identical  for  each  valne  of  <h.  Therefore,  $ 
will  be  set  to  5.8  x  10“®  kmol  m“^  as  originally  pnblished  for  the  simnlations  of  the 
X-2  experiments  presented  in  the  next  chapter.  In  addition,  the  assnmption  of  mobile 
adsorption  will  be  nsed  in  the  next  chapter  with  the  nnderstanding  that,  in  reality, 
some  adsorbed  atoms  conld  behave  as  immobile.  The  inclnsion  of  a  direct  nitridation 
reaction  in  the  ZA  model  will  also  be  considered  in  the  next  chapter  as  a  means  to 
improving  the  ablation  predictions  when  compared  to  the  X-2  experimental  resnlts. 
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Figure  36.  Total  surface  mass  flux  as  a  function  of  normalized  distance  from  stagnation 
point  at  different  surface  reaction  site  concentrations  for  the  30  km  conditions. 
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Figure  37.  CO2,  CO,  C3,  and  C2  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  stagnation  point  at  different  surface  reaction  site  concentrations  for  the  30  km 
conditions. 
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Figure  38.  C  surface  mass  fluxes  as  a  function  of  normalized  distance  from  stagnation 
point  at  different  surface  reaction  site  concentrations  for  the  30  km  conditions. 
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(a)  CO2  (b)  CO 


(c)  CN 


Figure  39.  Species  mass  fraction  as  a  function  of  surface-normal  distance  at  the  stag¬ 
nation  point  at  different  surface  reaction  site  concentrations  for  the  30  km  conditions. 
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V.  Results  and  Analysis 


The  results  shown  in  this  chapter  are  presented  in  a  progressive  order  of  the  work 
performed.  Simulations  were  hrst  performed  for  the  four  wall  temperature  cases  from 
the  hrst  X-2  campaign  shown  in  Table  3  using  the  Park,  ZA,  and  equilibrium  (i.e. 
B')  ablation  models.  The  Johnston  et  al.  [46]  kinetics  model  was  used  with  the 
equilibrium,  ZA  model,  and  the  Park  model  modihed  with  the  nitridation  efficiency 
proposed  by  Suzuki  et  al.  [115].  Based  on  those  results,  modihcations  to  the  ZA 
model  were  proposed  to  improve  CN  radiance  comparisons.  The  development  and 
implementation  of  those  modihcations  are  presented  below.  Finally,  simulations  of 
the  four  wall  temperature  cases  from  the  second  X-2  campaign  shown  in  Table  3  were 
completed  and  a  discussion  of  those  radiation  comparisons  is  provided. 

5.1  Grid  Convergence 

Initially,  a  grid  convergence  study  was  completed  to  ensure  accurate,  stable  re¬ 
sults  that  minimized  the  number  of  grid  points  and  computational  resources.  High 
resolution  of  the  mean  how  was  desired  in  the  front  surface  of  the  model  near  the 
stagnation  region.  Hence,  the  number  of  wall  normal  grid  points  was  increased  in  this 
region  until  mean  how  parameters  stabilized.  Specihcally  the  temperature,  CN  mass 
fraction,  and  total  surface  mass  hux  were  observed  with  increasing  grid  resolution. 
The  results  of  the  grid  convergence  study  are  shown  in  Figure  40.  The  number  of 
grid  points  in  the  wall  normal  direction  was  increased  from  150  to  200  and  250  points 
corresponding  to  the  total  number  of  cells  identihed  in  the  plots  as  4.6,  5.5,  and  6.4 
million  cells,  respectively.  Even  though  all  of  the  grids  generated  similar  temperature 
and  CN  mass  fraction  prohles,  there  were  differences  in  the  total  surface  mass  hux 
near  the  stagnation  point.  The  grid  with  5.5  million  cells  provided  similar  mass  huxes 
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(a)  Temperature 


(b)  CN  mass  fraction 


(c)  Surface  mass  flux 


(d)  Surface  mass  flux  near  stagnation  point 


Figure  40.  Convergence  of  mean  flow  parameters  with  increasing  grid  resolution  in  the 
wall  normal  direction.  Profiles  are  shown  along  the  stagnation  line. 
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Figure  41.  Unsteady  CN  mass  fraction  profiles  along  the  stagnation  line  at  different 
times  after  shock  arrival  for  =  2410  K. 


as  the  grid  with  6.4  million  cells,  so  the  following  results  used  a  grid  with  5.5  million 
cells  to  conserve  computational  time. 

Steady-state  simulations  were  run  because  of  the  process  in  which  the  experimen¬ 
tal  data  was  captured.  The  spectrometer  camera  delayed  for  at  least  20  /is  after  shock 
arrival  before  it  started  recording  measurements,  which  allowed  for  the  initial  tran¬ 
sient  behavior  to  relax  to  steady-state.  Conhrmation  of  steady-state  conditions  was 
achieved  by  analyzing  the  high-speed  video  and  by  performing  an  unsteady,  time- 
accurate  simulation.  Figure  41  displays  the  CN  mass  fraction  prohles  at  multiple 
instances  in  time  after  shock  arrival  for  the  case  with  a  wall  temperature  of  2410  K. 
The  flow  exhibits  steady  behavior  after  about  15  /rs  and  advancing  the  solution  further 
in  time  will  only  compute  the  same  CN  mass  fraction  prohle.  Therefore,  calculation 
of  steady-state  solutions  was  deemed  acceptable  for  these  simulations. 
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5.2  First  Campaign  Results 


5.2.1  Steady-State  Flow  Analysis. 

The  steady-state  solution  profiles  for  all  wall  temperatures  are  shown  in  Fig¬ 
ure  42.  The  shock  location  is  approximately  1.80  mm  from  the  wall,  and  the  thermal 
nonequilibrium  region  has  a  thickness  of  roughly  1  mm.  Although  thermal  equilib¬ 
rium  is  present  within  a  distance  of  about  0.8  mm  from  the  wall,  the  populations  of 
internal  energy  states  do  not  necessarily  follow  a  Boltzmann  distribution.  Laux  et 
ah  [53]  and  Johnston  [45]  demonstrated  that  the  excited  electronic  states  of  atoms 
are  underpopulated  by  a  factor  of  10  in  equilibrium  regions.  Due  to  these  reasons, 
a  non-Boltzmann  distribution  for  the  excited  electronic  states  was  assumed  for  all 
NEQAIR  simulations. 

The  total  surface  mass  flux  at  each  wall  temperature  is  plotted  in  Figure  43  as  a 
function  of  the  normalized  surface  distance  measured  from  the  stagnation  point.  The 
surface  distance  follows  the  model’s  horizontal  centerline  with  zero  corresponding  to 
the  stagnation  point  and  is  normalized  by  the  outer  diameter  of  the  model.  Overall, 
each  respective  ablation  model  predicted  similar  mass  fluxes  across  this  wall  tempera¬ 
ture  range.  For  example,  the  ZA  model  predicted  a  15%  drop  in  the  stagnation  point 
total  mass  flux  from  =  2410  K  to  1770  K,  while  the  Park  model  only  predicted 
a  1%  drop.  Reducing  the  nitridation  efficiency  in  the  Park  model  had  a  significant 
effect  on  the  total  mass  flux.  The  =  2410  K  case  showed  a  40%  stagnation  point 
mass  flux  reduction  between  the  original  and  modified  Park  models,  and  similar  re¬ 
sults  are  seen  at  the  other  wall  temperatures.  As  will  be  shown,  the  Park  model 
predicts  nitridation  to  be  the  dominant  surface  reaction  at  each  wall  temperature, 
which  is  why  there  is  a  significant  reduction  in  surface  mass  flux  with  the  modified 
Park  model.  Incorporating  the  Johnston  gas-phase  chemistry  model  had  a  minimal 
effect  and  had  a  tendency  to  slightly  reduce  the  mass  fluxes  when  using  the  ZA  or 
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Figure  42.  Steady-state  solution  profiles  for  the  X-2  cases  along  the  stagnation  line. 
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Surface  Mass  Flux  (kg/m^-s)  Surface  Mass  Flux  (kg/m^-s) 


Non-DImenslonal  Distance  on  Surface  (x/r„) 


(a)  =  2410  K 


(b)  =  2170  K 


Non-DImenslonal  Distance  on  Surface  (x/r„) 


(c)  =  1920  K 


(d)  =  1770  K 


Figure  43.  Predicted  total  surface  mass  flux  as  a  function  of  normalized  distance  from 
the  stagnation  point  that  follows  the  model’s  horizontal  centerline. 
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modified  Park  ablation  model.  Considering  the  eqnilibrinm  air-carbon  mixtnre  resnlt 
shown  in  Fignre  15b,  the  identical  total  snrface  mass  finx  profiles  for  the  eqnilibrinm 
ablation  assnmption  was  as  expected  for  this  wall  temperatnre  range.  The  eqnilib- 
rinm  total  snrface  mass  finxes  fell  between  the  ZA  and  modified  Park  ablation  model 
predictions. 

There  are  notable  differences  in  the  oxidation  and  nitridation  mechanisms  between 
the  ablation  models.  The  primary  ablation  prodncts  predicted  by  the  snrface  models 
were  CO2,  CO,  and  CN.  Fignres  44-47  plot  the  snrface  mass  finxes  for  these  species. 
Within  this  wall  temperatnre  and  pressnre  range,  there  was  no  predicted  finx  of  C3, 
C2,  and  C  for  any  model  considered.  A  positive  mass  finx  indicates  that  the  species 
is  being  formed  at  the  snrface  and  diffnses  away.  A  negative  mass  finx  means  the 
species  is  diffnsing  to  the  snrface  where  it  is  adsorbed  or  reacts  with  the  bnlk  carbon 
snrface  or  other  adsorbed  snrface  species.  For  all  wall  temperatnres  considered,  the 
ZA  model  predicted  significant  CO2  prodnction  and  had  CO  consnmption  at  the 
snrface.  In  contrast,  the  Park  and  eqnilibrinm  models  predicted  CO  prodnction  as 
the  only  oxidation  mechanism.  The  Park  model  predicted  CN  to  be  the  primary 
species  prodnced  at  the  snrface,  while  the  ZA  model  has  no  nitridation  mechanism, 
so  it  had  zero  CN  mass  finx.  The  eqnilibrinm  model  also  predicted  zero  CN  mass  finx 
at  all  snrface  temperatnres.  Applying  the  Snznki  nitridation  efficiency  to  the  Park 
model  changes  the  relative  importance  of  the  nitridation  and  oxidation  reactions.  The 
modified  Park  model  dropped  the  CN  finx  by  86%  and  more  than  donbled  the  CO 
prodnction  rate  for  the  =  2410  K  case.  The  .Johnston  chemistry  model  had  a 
negligible  effect  on  the  modified  Park  model  resnlts,  bnt  it  did  increase  the  CO2  mass 
finx  predicted  by  the  ZA  model  at  all  wall  temperatnres. 

The  species  mass  fractions  in  the  near  wall  bonndary  layer  refiect  the  snrface 
mass  finx  resnlts  and  will  have  a  direct  impact  on  the  predicted  radiative  heat  finxes. 
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Figure  44.  Predicted  species  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  the  stagnation  point  for  =  2410  K. 
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Figure  45.  Predicted  species  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  the  stagnation  point  for  =  2170  K. 
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Figure  46.  Predicted  species  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  the  stagnation  point  for  =  1920  K. 
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Figure  47.  Predicted  species  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  the  stagnation  point  for  =  1770  K. 
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Figure  48.  Predicted  species  mass  fractions  along  the  stagnation  streamline  for 

2410  K. 
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Figure  49.  Predicted  species  mass  fractions  along  the  stagnation  streamline  for 

2170  K. 
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Figure  50.  Predicted  species  mass  fractions  along  the  stagnation  streamline  for 

1920  K. 
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Figure  51.  Predicted  species  mass  fractions  along  the  stagnation  streamline  for 
1770  K. 
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Figures  48-51  plot  the  species  mass  fractions  for  CO2,  CO,  and  CN  as  a  function 
of  wall  normal  distance  along  the  stagnation  streamline.  The  modihed  Park  and 
equilibrium  models  gave  very  similar  mass  fraction  prohles  at  all  surface  temperatures. 
Within  the  considered  surface  temperature  range  where  oxidation  processes  dominate, 
the  modihed  Park  and  equilibrium  models  have  comparable  behavior.  Both  hnite- 
rate  surface  models  have  mechanisms  to  produce  CO  but  do  so  through  different 
kinetic  processes.  The  Park  model  produced  CO  and  predicted  it  to  diffuse  away 
from  the  surface  resulting  in  a  larger  mass  fraction  in  the  near-wall  boundary  layer.  In 
comparison,  the  ZA  model  predicted  CO  to  recombine  at  the  surface  resulting  in  lower 
CO  boundary  layer  concentrations.  Additionally,  the  predicted  CO2  mass  fractions 
were  substantially  different  in  the  near-wall  boundary  layer  region  between  the  hnite- 
rate  surface  models.  The  ZA  model  predicted  over  three  orders  of  magnitude  more 
CO2  at  the  stagnation  point  compared  to  the  Park  model  for  all  wall  temperatures. 
Interestingly,  using  the  Johnston  chemistry  with  the  ZA  model  increased  the  rate  at 
which  CO2  was  destroyed  along  the  stagnation  streamline,  but  it  actually  increased 
the  concentration  of  CO  near  the  wall.  The  increase  in  gaseous  CO  concentration  is 
non-intuitive  because  the  Johnston  chemistry  model  has  a  CO  dissociation  rate  13 
times  higher  than  the  Park  gas-phase  chemistry  model  rate.  Therefore,  it  would  be 
expected  the  amount  of  CO  in  the  boundary  layer  would  decrease  using  the  Johnston 
chemistry  model.  However,  as  noted  in  Section  3.3.7,  the  ZA  surface  reactions  are  non- 
linearly  coupled  to  the  gas  chemistry  model,  so  results  displaying  linear  relationships 
should  not  be  expected. 

Similar  CN  mass  fraction  prohles  were  achieved  for  all  cases  at  wall  normal  dis¬ 
tances  greater  than  approximately  1  mm,  which  coincided  with  the  region  of  thermal 
nonequilibrium.  The  CN  concentration  for  the  Park  model  is  strongly  affected  by  the 
amount  of  CN  produced  at  the  wall.  Consequently,  the  Park  model  predicted  the 
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highest  boundary  layer  CN  concentration  across  the  wall  temperature  range  consid¬ 
ered.  Due  to  the  much  lower  nitridation  rate  of  the  modihed  Park  model,  the  peak 
CN  concentration  in  the  boundary  layer  was  approximately  cut  in  half  compared  to 
the  original  Park  model  for  each  wall  temperature.  The  ZA  model,  which  has  no 
nitridation  mechanism,  predicted  the  lowest  peak  CN  concentrations.  The  amount  of 
CN  in  the  boundary  layer  was  entirely  controlled  by  the  gas-phase  chemistry  model 
for  the  ZA  results.  As  the  wall  temperature  was  lowered,  the  peak  boundary  layer 
CN  concentration  predicted  by  the  ZA  model  was  reduced  by  18%  from  =  2410  K 
to  1770  K.  Comparing  the  mass  fraction  prohles  for  all  cases  considered,  there  ap¬ 
pears  to  be  a  correlation  between  the  CO  and  CN  boundary  layer  concentrations.  In 
general,  higher  CO  concentrations  led  to  higher  CN  concentrations. 

5.2.2  Radiation  Comparisons. 

Good  agreement  was  found  for  the  locations  of  peak  radiation  intensity  measured 
in  the  experiment  and  the  predicted  maximum  CN  mass  fraction  in  the  boundary 
layer.  Figure  52  demonstrates  this  agreement  with  a  series  of  plots.  The  predicted 
CN  mass  fraction  shown  in  this  hgure  was  computed  using  the  ZA  model.  The 
spectrometer  measured  the  strongest  radiative  intensity  at  the  stagnation  point  to 
occur  approximately  between  0.10-0.30  mm  from  the  wall.  The  spectral  radiance 
was  extracted  along  a  line  of  constant  wavelength  at  388.4  nm  and  plotted  with  the 
CN  mass  fraction  in  Figure  52c.  As  will  be  shown  below,  the  SRF  has  a  tendency 
to  broaden  and  shift  the  predicted  radiance  prohle  closer  to  the  wall  which  further 
improves  this  comparison.  There  is  also  good  agreement  for  the  shock  location,  which 
is  approximately  located  1.80  mm  away  from  the  wall.  Note  that  the  predicted  local 
maximum  of  CN  in  the  boundary  layer  was  also  at  a  similar  wall  normal  distance  when 
computed  using  the  Park  and  equilibrium  models.  Hence,  all  surface  ablation  models 
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(a)  Predicted  CN  mass  fraction. 


(b)  Experimental  spectral  radiance. 


(c)  Comparison  of  peak  and  shock  locations. 


Figure  52.  Comparison  of  peak  experimental  spectral  radiance  and  predicted  maximum 
CN  concentration  using  the  ZA  model  for  the  =  2410  K  case. 
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are  consistent  with  the  experiment  in  this  regard,  bnt  it  is  the  actnal  concentration 
level  that  will  have  the  most  signihcant  impact  on  the  radiation  resnlts. 

The  radiance  resnlts  as  a  fnnction  of  wall  normal  distance  along  the  stagnation 
streamline  are  compared  in  Fignre  53.  The  X-2  radiance  nncertainty  shown  in  this 
hgnre  varied  between  31%-33%  as  dehned  by  Lewis  et  ah  [57],  Note  that  the  same 
estimated  SRF  was  applied  to  each  predicted  radiance  prohle  calcnlated  by  NEQAIR, 
which  broadened  and  slightly  shifted  the  prohles  toward  the  snrface.  The  resnlt  was 
a  mnch  better  agreement  to  the  shape  and  location  of  the  bonndary  layer  peak  valne 
compared  to  the  measnred  radiance  prohle. 

Overall,  the  nnmerical  simnlations  demonstrated  a  better  ability  to  predict  radi¬ 
ance  valnes  near  the  snrface  than  in  the  shock  layer.  The  experimental  measnrements 
displayed  a  more  signihcant  amonnt  of  CN  violet  radiation  emitted  from  the  shock 
layer  than  what  the  simnlations  indicated.  The  reasons  for  this  are  cnrrently  nn- 
known,  bnt  there  is  likely  some  sonrce  of  facility  contamination  present.  It  was  pecn- 
liar  that  the  shock  layer  CN  radiation  appeared  to  increase  with  snrface  temperatnre, 
and  Lewis  et  al.  [57]  specnlated  that  this  conld  be  dne  to  spallation.  The  high-speed 
video  from  the  =  2410  K  case  possibly  displayed  a  small  amonnt  of  carbon  particles 
ejected  from  the  snrface.  The  theory  is  that  spalled  carbon  particles  entered  the  invis- 
cid,  high  temperatnre  shock  layer  and  vaporized  allowing  carbonaceons  gas  species, 
snch  as  CN,  to  form  and  radiate.  These  simnlations  can  not  predict  the  mechanical 
breakdown  of  the  bnlk  carbon,  nor  track  spalled  particles  as  they  diffnse  throngh  the 
shock  layer.  Another  consideration  may  be  that  there  are  trace  amonnts  of  CN  in 
the  freestream  that  emits  as  soon  as  it  passes  throngh  the  shock.  Possible  sonrces 
conld  be  from  the  facility’s  vacnnm  pnmp  system  or  diaphragm  materials.  More  in¬ 
vestigation  into  this  phenomenon  is  reqnired  to  determine  the  discrepancy  between 
the  experiments  and  simnlations.  Therefore,  the  following  comparison  discnssion  will 
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Radiance  (W/cm^-sr)  Radiance  (W/cm^-sr) 


(a)  =  2410  K 


(b)  =  2170  K 


(c)  T„  =  1920K 


(d)  =  1770  K 


Figure  53.  Comparison  of  radiance  values  as  a  function  of  wall  normal  distance  along 
the  stagnation  streamline. 
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focus  on  the  radiance  profiles  near  the  model  surface. 

The  following  trends  observed  between  the  two  ablation  models  indicate  that 
nitridation  may  be  an  important  surface  reaction  to  consider  in  these  flows.  In  ac¬ 
cordance  with  the  relatively  higher  CN  concentrations  in  the  boundary  layer,  the 
Park  model  predicted  radiance  values  above  those  measured  in  the  experiment  for 
all  considered  surface  temperatures.  However,  the  conclusion  is  not  that  nitridation 
was  an  insignificant  surface  reaction.  The  modified  Park  model  displayed  reasonably 
good  agreement  at  =  2410  K  with  marginal  improvement  when  combined  with 
the  Johnston  chemistry  model.  The  ablation  of  CN  into  the  boundary  layer  increased 
the  CN  mass  fraction  to  concentrations  found  in  the  experiment  for  the  =  2410  K 
case.  Unfortunately,  the  comparison  got  worse  for  the  modified  Park  model  as  the 
wall  temperature  was  lowered,  which  does  lead  to  the  conclusion  that  nitridation 
becomes  less  significant.  The  ZA  model  was  shown  to  exhibit  the  opposite  trend 
and  had  the  best  agreement  at  the  lowest  wall  temperature  considered.  Although 
the  applied  SRF  is  more  or  less  arbitrary  for  this  X-2  data  set,  the  ZA  model  with 
Johnston  chemistry  had  excellent  agreement  for  the  near  surface  radiance  profile  at 

=  1770  K.  The  ZA  model  then  began  to  under  predict  the  radiance  values  at 
the  higher  wall  temperatures.  Based  off  these  observations,  the  nitridation  reaction 
appears  to  be  activated  above  a  surface  temperature  of  approximately  1800  K. 

The  equilibrium  result  gave  the  same  radiance  profile  for  each  surface  temperature. 
Clearly,  the  radiation  emitted  from  CN  was  sensitive  to  a  varying  surface  temperature 
as  measured  in  the  experiment.  Therefore,  the  use  of  an  equilibrium  ablation  model  is 
not  appropriate  to  simulate  these  specific  flows  even  though  there  is  the  appearance  of 
good  agreement  at  a  few  of  the  surface  temperatures.  At  =  2170  K,  for  example, 
the  equilibrium  model  provided  the  best  prediction  compared  to  the  experiment,  but 
it  would  be  presumptuous  to  conclude  that  the  surface  and  gas  were  indeed  in  chemical 
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equilibrium. 

To  quantify  these  differences,  the  normalized  absolute  error  with  respect  to  the 
experimental  values  was  calculated  for  each  predicted  result.  The  error  is  only  shown 
out  to  a  wall  normal  distance  of  0.4  mm  due  to  the  uncertainties  in  the  shock  layer. 
The  Park  model  results  are  not  shown  at  a  wall  temperature  of  1770  K  because 
its  error  was  around  400%.  These  error  results  reinforce  the  observed  trends.  At 
Tw  =  2410  K,  the  modihed  Park  model  with  .Johnston  chemistry  had  a  maximum 
error  of  approximately  8%,  while  the  ZA  and  original  Park  models  hovered  around  an 
error  between  60-70%.  The  ZA  and  modihed  Park  models  both  had  comparable  error 
estimates  with  a  mean  around  50%  for  wall  temperatures  of  2170  K  and  1920  K.  At 
the  lowest  wall  temperature  of  1770  K,  the  ZA  model  with  the  .Johnston  chemistry 
had  errors  of  1-2%  that  coincided  with  the  location  of  peak  intensity.  The  modihed 
Park  models  began  to  signihcantly  diverge  from  the  experimental  results  with  errors 
above  180%. 

From  this  analysis,  it  is  concluded  that  the  ZA  model  is  preferred  over  the  original 
or  modihed  Park  model  because  it  is  a  mechanism-based  approach  with  a  physically 
meaningful  methodology.  The  Park  models  are  too  simplistic  with  highly  empirical 
rate  parameters  that  only  include  a  few  of  the  relevant  ablative  species  for  an  air- 
carbon  system.  The  ZA  model  makes  use  of  gas-kinetic  theory  and  accounts  for 
physical  mechanisms  such  as  the  hux  of  atoms  impinging  on  the  surface  and  adsorption 
processes.  Furthermore,  many  more  reaction  pathways  are  dehned  with  the  ZA  model 
that  includes  more  of  the  relevant  ablative  species.  However,  parameterizing  the  ZA 
model  remains  an  issue  because  the  predicted  species  mass  fluxes  were  nearly  constant 
between  surface  temperatures  of  1770-2410  K,  which  contradicted  the  trends  observed 
in  the  experiment.  Additionally,  there  may  be  a  dehciency  in  the  baseline  ZA  model 
because  it  does  not  include  any  CN  surface  reactions.  The  following  sections  describe 
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Figure  54.  Normalized  absolute  error  of  predicted  radiance  values  with  respect  to  the 
experimental  results. 
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how  a  new  model  derived  from  the  baseline  ZA  model  is  developed  to  fnrther  improve 
the  CN  radiance  comparisons. 

5.3  Screening  Design  Results 

5.3.1  Statistical  Correlations. 

The  ZA  model  rate  coefficients  screening  design  described  in  Section  3.7  was  exe- 
cnted  nsing  the  30  km  condition  from  Section  4.1  becanse  that  condition  was  shown 
to  inclnde  all  three  ablation  mechanisms  of  oxidation,  nitridation,  and  snblimation. 
For  a  more  complete  screening  design,  the  direct  nitridation  reaction  was  inclnded  in 
the  ZA  model  becanse  it  was  shown  in  Section  4.3  that  the  snrface  participating  nitri¬ 
dation  reaction  had  a  negligible  effect  on  the  resnlts.  The  baseline  ZA  rate  coefficient 
valnes  were  shown  in  Tables  8  and  10.  To  compnte  the  CN  radiative  heat  flnxes,  the 
NEQAIR  code  was  nsed  again.  For  the  pnrposes  of  this  screening  design,  the  popnla- 
tions  of  the  excited  CN  electronic  states  were  fonnd  nsing  a  Boltzmann  distribntion. 
The  radiance  was  calcnlated  by  integrating  the  spectral  radiance  between  the  wave¬ 
length  range  of  353-391  nm.  The  post-shock  flow  qnantities  were  extracted  from  the 
volnmetric  solntion  corresponding  to  a  line-of-sight  that  was  normal  to  the  stagnation 
streamline.  Mnltiple  slice  extractions  were  performed  along  the  stagnation  streamline 
that  provided  radiance  predictions  as  a  fnnction  of  wall  normal  distance.  To  obtain 
a  single  total  CN  radiance  valne,  the  radiance  along  the  stagnation  streamline  was 
integrated  from  the  snrface  ont  to  a  wall  normal  distance  of  0.4  mm. 

There  were  13  screening  design  factors  that  represented  each  pre-exponential  re¬ 
action  rate  coefficient  of  the  ZA  model  inclnding  nitridation.  Table  15  shows  the 
screening  design  matrix  and  resnlting  metrics.  Again,  the  qnantitative  metric  was 
the  absolnte  percent  difference  of  the  CN  radiance  from  the  baseline  case  (Rnn  1). 
For  simplicity,  the  13  factor  levels  were  denoted  by  the  symbol  Aj  with  the  snbscript 
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Table  15.  The  Zhluktov  and  Abe  rate  coefficient  screening  design  and  CN  radiance 
metrics. 


Factors  ^  Levels 

Metrics 

Run 

Al 

A2 

A3 

A4 

A5 

Ae 

A7 

As 

Ag 

Aio 

All 

Ai2 

Ai3 

CN  Radiance 

Difference 

1 

1.0 

1.0 

3.58 

1.0 

0.90 

2.08  X  10* 

0.80 

3.58 

X  10“ 

0.240 

0.50 

0.0230 

1.0 

0.360 

0.84667 

0.0000 

2 

1.0 

1.1 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.264 

0.55 

0.0253 

1.1 

0.396 

0.85440 

0.0091 

3 

1.0 

0.9 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.216 

0.45 

0.0207 

0.9 

0.324 

0.82614 

-0.0242 

4 

1.1 

1.0 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.264 

0.55 

0.0207 

0.9 

0.396 

0.83630 

-0.0122 

5 

0.9 

1.0 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.216 

0.45 

0.0253 

1.1 

0.324 

0.85355 

0.0081 

6 

1.1 

0.9 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.264 

0.55 

0.0253 

0.9 

0.324 

0.83324 

-0.0159 

7 

0.9 

1.1 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.216 

0.45 

0.0207 

1.1 

0.396 

0.85495 

0.0098 

8 

1.1 

1.1 

3.58 

X  10“ 

1.0 

0.81 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.264 

0.55 

0.0253 

1.1 

0.324 

0.85013 

0.0041 

9 

0.9 

0.9 

3.58 

X  10“ 

1.0 

0.99 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.216 

0.45 

0.0207 

0.9 

0.396 

0.84084 

-0.0069 

10 

1.1 

0.9 

3.58 

X  10“ 

0.9 

0.90 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.216 

0.55 

0.0253 

1.1 

0.396 

0.84017 

-0.0077 

11 

0.9 

1.1 

3.58 

X  10“ 

1.1 

0.90 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.264 

0.45 

0.0207 

0.9 

0.324 

0.84595 

-0.0009 

12 

1.1 

0.9 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.216 

0.45 

0.0253 

1.1 

0.396 

0.85367 

0.0083 

13 

0.9 

1.1 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.264 

0.55 

0.0207 

0.9 

0.324 

0.84037 

-0.0074 

14 

1.1 

1.1 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10* 

0.80 

3.58 

X  10“ 

0.264 

0.45 

0.0207 

1.1 

0.396 

0.82980 

-0.0199 

15 

0.9 

0.9 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10“ 

0.80 

3.58 

X  10“ 

0.216 

0.55 

0.0253 

0.9 

0.324 

0.85431 

0.0090 

16 

1.1 

1.1 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.216 

0.55 

0.0207 

0.9 

0.396 

0.85592 

0.0109 

17 

0.9 

0.9 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.264 

0.45 

0.0253 

1.1 

0.324 

0.82584 

-0.0246 

18 

1.1 

1.1 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.240 

0.45 

0.0253 

0.9 

0.324 

0.82702 

-0.0232 

19 

0.9 

0.9 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.240 

0.55 

0.0207 

1.1 

0.396 

0.85303 

0.0075 

20 

1.1 

1.1 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.216 

0.50 

0.0207 

1.1 

0.324 

0.83419 

-0.0147 

21 

0.9 

0.9 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.264 

0.50 

0.0253 

0.9 

0.396 

0.85400 

0.0087 

22 

1.1 

0.9 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10“ 

0.72 

3.58 

X  10“ 

0.264 

0.45 

0.0230 

0.9 

0.396 

0.85227 

0.0066 

23 

0.9 

1.1 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10* 

0.88 

3.58 

X  10“ 

0.216 

0.55 

0.0230 

1.1 

0.324 

0.83626 

-0.0123 

24 

1.1 

0.9 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.216 

0.55 

0.0207 

1.0 

0.324 

0.85130 

0.0055 

25 

0.9 

1.1 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.264 

0.45 

0.0253 

1.0 

0.396 

0.83912 

-0.0089 

26 

1.1 

1.1 

3.58 

X  10“ 

0.9 

0.99 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.216 

0.45 

0.0253 

0.9 

0.360 

0.85622 

0.0113 

27 

0.9 

0.9 

3.58 

X  10“ 

1.1 

0.81 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.264 

0.55 

0.0207 

1.1 

0.360 

0.82612 

-0.0243 

28 

1.1 

0.9 

3.58 

X  10“ 

0.9 

0.81 

2.08  X  10“ 

0.88 

3.58 

X  10“ 

0.264 

0.45 

0.0207 

1.1 

0.324 

0.84619 

-0.0006 

29 

0.9 

1.1 

3.58 

X  10“ 

1.1 

0.99 

2.08  X  10* 

0.72 

3.58 

X  10“ 

0.216 

0.55 

0.0253 

0.9 

0.396 

0.84026 

-0.0076 

Table  16.  Correlations  for  Zhluktov  and  Abe  rate  coefficient  screening.  Significant 
correlations  are  highlighted. 


Al 

A2 

A3 

A4 

A5 

Ae 

A7 

Ag 

Ag 

Alo 

A„ 

Ai2 

Ai3 

CN  Radiance 

0.0065 

0.0267 

0.0060 

-0.0265 

0.0185 

0.8222 

0.0010 

0.1348 

-0.2285 

0.0722 

0.1457 

-0.0162 

0.2863 

Difference 

0.0070 

0.0269 

0.0052 

-0.0266 

0.0181 

0.8219 

0.0015 

0.1344 

-0.2287 

0.0719 

0.1453 

-0.0163 

0.2865 

i  referring  to  the  reaction  number.  The  metrics  were  evaluated  in  a  JMP®  [104] 
multivariate  analysis  report.  The  report  provided  independent  correlation  values  for 
each  factor  and  response,  which  are  shown  in  Table  16.  The  pairwise  correlations 
gave  an  estimate  of  factor  signihcance.  Lehman  et  al.  [55]  proposed  that  correlations 
between  0.2  and  0.5  are  weakly  significant,  correlations  between  0.5  and  0.8  are  mod¬ 
erately  signihcant,  and  correlations  greater  than  0.8  are  strongly  significant.  These 
correlations  are  also  applicable  for  negative  values  where  they  would  be  dehned  as  an 
inverse  relationship. 

The  only  statistically  significant  factor  that  had  an  impact  on  the  CN  radiance 
was  Ag,  the  Arrhenius  coefficient  for  the  oxidation  reaction  0(s)  +C{b)  ^  CO  +  (s). 
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The  nitridation  efficiency  and  atomic  carbon  condensation  efficiency,  factors  A13  and 
Ag,  were  the  only  other  parameters  that  conld  possibly  be  identihed  as  somewhat 
signihcant.  However,  the  correlations  indicate  that  any  changes  to  these  coefficients 
wonld  only  have  minor  impacts.  Note  that  the  atomic  carbon  condensation  efficiency 
is  negatively  correlated  meaning  a  decrease  in  efficiency  wonld  resnlt  in  an  increase 
in  CN  radiance  valnes.  Fnrthermore,  these  resnlts  show  that  more  CO  needs  to  be 
prodnced  at  the  snrface  if  the  goal  is  to  increase  CN  radiance  valnes. 

5.3.2  Surface  Reaction  Rate  Modifications. 

Based  on  the  rate  coefficient  screening  resnlts,  rate  modffications  shonld  consider 
the  snrface  reactions  that  had  the  most  inflnence  on  increasing  CN  radiance  predic¬ 
tions.  The  initial  set  of  modihcations  applied  to  the  ZA  model  are: 

•  Increase  the  Arrhenins  coefficient  by  at  least  one  order  of  magnitnde  for  the 
oxidation  reaction  0(s)  +  C{b)  ^  CO  +  (s). 

•  Increase  the  ER  nitridation  efficiency  by  no  more  than  10%  becanse  there  is  still 
nncertainty  abont  extrapolating  the  derived  rate  parameters  to  higher  snrface 
temperatnres. 

•  Decrease  the  atomic  carbon  condensation  efficiency  by  10%. 

The  inflnence  of  CO2  shonld  also  be  considered  becanse  it  conld  indirectly  effect  CN 
concentrations.  The  stagnation  point  reaction  flnxes  from  Rnn  1  of  the  screening 
design  (baseline  case)  are  shown  in  Table  17  to  highlight  preferred  reaction  direc¬ 
tions.  The  snrface  temperatnre  was  approximately  3600  K.  Under  these  conditions, 
reaction  5  prodnces  CO2  becanse  it  proceeds  in  the  backward  direction.  Addition¬ 
ally,  Fignres  44-47  showed  that  CO2  was  the  dominant  ablative  species  for  each 
snrface  temperatnre  considered.  Therefore,  a  modihcation  to  the  snrface  reaction 
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Table  17.  Stagnation  point  reaction  flnxes  from  Rnn  1  of  the  Zhlnktov  and  Abe  rate 
coefficient  screening  design  (T^  «  3600  K). 


Reaction 

r  (kmol  m  ^  s 

(1) 

0  ±  (s)  <->•  0(s) 

4.53  X  10-3 

(2) 

N  ±  (s)  ^  N(s) 

1.16  X  10“^ 

(3) 

20(s)  ^02±2(s) 

2.16  X  10-3 

(4) 

O2  ±  (s)  ^  0  ±  0(s) 

-2.17  X  10-6 

(5) 

C02±(s)  ^CO±0(s) 

-3.47  X  10-^ 

(6) 

0(s)  ±  c[h)  ^  CO  ±  (s) 

4.12  X  10-3 

(7) 

0  ±0(s)  ±C(6)  ^  CO2  ±  (s) 

4.29  X  10-6 

(8) 

20(s)  ±C(6)  ^  CO2  ±2(s) 

2.59  X  10-6 

(9) 

C  ±  (s)  <—>  (s)  ±  C(b) 

2.10  X  10-2 

(10) 

C2  ±  2(s)  ^  2(s)  ±  2C(&) 

-1.77  X  10-3 

(11) 

C3  ±  3(s)  ^  3(s)  ±  3C(&) 

-1.42  X  10-3 

(12) 

N2  ±  (s)  ^  N  ±  N(s) 

-1.16  X  10-^ 

(13) 

N±  (s)  ±C(6)  ^  CN±  (s) 

2.76  X  10-2 

CO2  +  (s)  ^  CO  +  0(s)  is  recommended  because  it  appears  to  have  a  tendency  to 
generate  CO2  rather  than  CO.  The  screening  design  results  showed  that  the  forward 
rate  ER  efficiency  for  this  reaction  had  no  correlation  to  predicted  CN  radiance  val¬ 
ues.  However,  this  reaction  prefers  to  proceed  in  the  reverse  direction  at  a  surface 
temperature  of  approximately  3600  K,  so  a  ±10%  forward  rate  perturbation  is  not 
going  to  change  the  reaction  direction.  Hence,  the  statistical  insignfficance  of  this 
reaction  makes  sense  because  any  slight  increase  or  decrease  in  CO2  mass  flux  is  not 
going  to  alter  CN  levels. 

For  the  following  reasons,  and  without  additional  detailed  reaction  rate  data,  it 
is  proposed  to  remove  the  CO2  ±  (s)  ■«-»•  CO  ±  0(s)  reaction  from  the  ZA  model  in 
hopes  of  increasing  CN  radiance  values.  A  review  of  the  ZA  model  revealed  that 
the  incorporated  oxidation  reactions  were  influenced  by  the  work  found  in  Refer¬ 
ences  [123,  111,  85,  50,  79].  References  [50]  and  [79]  only  describe  oxygen  and 
nitrogen  surface  recombination  processes  and  have  no  reactions  to  produce  CO  or 
CO2.  References  [123,  111,  85]  discuss  several  possible  mechanisms  for  CO  and  CO2 
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Figure  55.  Species  mass  fractions  of  an  equilibrium  carbon-air  mixture  as  a  function  of 
temperature  at  1  atm  pressure. 


production,  but  show  it  as  only  possible  through  reactions  of  adsorbed  O  or  O2 
with  the  bulk  carbon.  In  fact,  Zhluktov  and  Abe  [136]  identify  the  work  of  Walker 
et  ah  [123]  as  the  preferred  mechanism  for  CO2  formation,  of  which  the  reaction 
20(s)  +  C(6)  ^  CO2  +  2(s)  is  primarily  responsible.  Walker  et  al.  [123]  does  list  the 
CO  +  0(s)  <->■  CO2  +  (s)  reaction  as  a  possible  CO2  mechanism  but  concluded  that  it 
was  not  likely  to  occur  based  on  their  experiments.  None  of  the  other  references  list 
any  possible  path  similar  to  this  reaction.  Expanding  the  review  to  other  experimen¬ 
tal  work  analyzing  carbon  oxidation  discovered  that  CO  was  consistently  found  as 
the  principal  reaction  product  across  a  range  of  surface  pressures  and  temperatures 
[78,  124,  83,  82,  58,  100,  77].  Additionally,  MacLean  [63]  developed  a  tool  to  compute 
the  composition  of  an  equilibrium  saturated  air-carbon  mixture  at  selected  pressures 
and  temperatures.  The  species  mass  fractions  of  an  equilibrium  air-carbon  mixture 
as  a  function  of  temperature  at  1  atm  pressure  is  plotted  in  Figure  55,  which  was  the 
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stagnation  pressnre  in  the  X-2  experiments  shown  in  Fignre  42.  Note  that  the  gas  is 
a  simple  mixtnre  of  primarily  N2  and  CO  within  an  approximate  temperatnre  range 
of  1500-2500  K. 

The  proposed  modihcations  were  implemented  into  the  ZA  model  and  applied 
to  the  30  km  case.  A  comparison  of  the  species  mass  flnxes  with  the  Park  and 
modihed  ZA  models  is  shown  in  Fignre  56.  The  C2  and  C3  mass  flnx  prohles  were 
nnchanged  becanse  none  of  the  proposed  modihcations  affected  those  snblimation 
reactions.  Additionally,  the  C  and  CN  mass  hnxes  remained  relatively  the  same.  The 
major  improvement  with  this  proposed  model  compared  to  the  other  model  variations 
is  that  CO  is  now  the  main  oxidation  prodnct  instead  of  CO2  over  a  majority  of 
the  snrface.  However,  despite  increased  CO  snrface  prodnction,  the  calcnlated  total 
CN  radiance  was  not  that  mnch  higher  compared  to  the  baseline  case  with  only  a 
0.001  difference.  The  small  difference  is  likely  related  to  the  gas-phase  chemistry 
model,  which  has  a  strong  inhnence  on  species  prodnction  rates  in  the  bonndary 
layer.  The  Johnston  gas  chemistry  model  nsed  in  the  X-2  simnlations  increases  the  CO 
dissociation  rate  by  a  factor  of  13  compared  to  the  model  nsed  for  the  screening  design 
stndy.  Combining  these  ZA  model  modihcations  with  the  Johnston  gas  chemistry 
model  may  now  prodnce  drastically  different  resnlts. 

5.4  ZA  Model  Modifications  Applied  to  X-2  Experiments 

The  fonr  wall  temperatnre  cases  from  the  hrst  X-2  campaign  were  rernn  with  the 
proposed  modihcations  described  in  the  previons  section.  Fignres  57  and  58  plot 
the  predicted  species  snrface  mass  hnxes  and  stagnation  streamline  mass  fractions 
at  Tw  =  2410  K.  The  ZA  model  with  nitridation  now  has  a  CN  mass  hnx  that 
is  approximately  three  times  larger  than  the  modihed  Park  model.  However,  CN 
dissociates  immediately  away  from  the  snrface,  and  there  is  only  a  slight  increase  in 
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Figure  56.  Species  surface  mass  fluxes  as  a  function 
nation  point  at  30  km  conditions  with  proposed  ZA 


of  normalized  distance 
model  modiflcations. 


from  stag- 
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the  peak  boundary  layer  concentration  compared  to  the  ZA  model  without  nitridation 
as  shown  in  Figure  58.  The  inclusion  of  a  nitridation  reaction  does  not  signihcantly 
alter  the  oxidation  reactions,  and  CO2  is  still  the  main  product. 

There  was  a  signihcant  change  when  all  of  the  proposed  nitridation  and  oxidation 
reaction  modihcations  were  implemented.  The  main  oxidation  product  becomes  CO, 
and  very  little  CO2  was  produced  at  the  surface.  The  CO  mass  fractions  at  the  surface 
were  an  order  of  magnitude  higher  and  were  similar  to  the  modihed  Park  model 
predictions.  It  was  previously  noted  that  higher  CO  concentrations  led  to  higher  CN 
concentrations,  which  is  again  reflected  in  Figure  58.  The  peak  CN  boundary  layer 
concentration  was  twice  as  high  with  the  ZA  model  with  all  proposed  modihcations 
compared  to  the  baseline  ZA  model.  The  increase  in  CN  levels  is  likely  due  to  a 
combination  of  these  three  gas-phase  reactions: 

CO  +  M  C  +  O  +  M 
N2  +  C  ^  CN  +  N 
CO  +  N  CN  +  0 

The  Johnston  [46]  gas  chemistry  model  dissociates  CO  at  a  much  higher  rate  than 
the  traditional  Park  [96]  model,  which  will  make  C  more  readily  available  to  react 
with  N2.  Any  leftover  CO  can  react  with  N  to  also  form  CN. 

The  ZA  model  augmented  with  direct  nitridation  only  provided  a  minimal  im¬ 
provement  in  radiance  predictions  at  the  higher  considered  surface  temperatures. 
Updated  CN  radiance  predictions  are  presented  in  Figure  59  with  the  implemented 
direct  nitridation  reaction  and  the  rate  modihcations  proposed  in  the  previous  sec¬ 
tion.  Note  that  all  surface  model  predictions  shown  used  the  Johnston  gas  chemistry 
model.  At  surface  temperatures  of  1920  K  and  1770  K,  the  ZA  model  with  direct 
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(a)  CO2  (b)  CO 


(c)  CN 


Figure  57.  Predicted  species  surface  mass  fluxes  as  a  function  of  normalized  distance 
from  the  stagnation  point  for  the  =  2410  K  case. 
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(a)  CO2  (b)  CO 


(c)  CN 


Figure  58.  Predicted  species  mass  fractions  along  the  stagnation  streamline  for 
2410  K  case. 
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(d)  =  1770  K 


Figure  59.  Updated  comparison  of  CN  radiance  values  as  a  function  of  wall  normal 
distance  along  the  stagnation  streamline  for  lower  four  wall  temperature  cases  from 
the  X-2  experiments.  The  ZA  model  has  been  modified  to  include  direct  nitridation 
and  other  rate  coefficient  adjustments. 
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nitridation  gave  the  same  result  as  the  baseline  ZA  model,  which  is  consistent  with 
the  observations  from  Section  4.3  that  nitridation  was  only  active  above  surface  tem¬ 
peratures  of  about  2000  K.  The  ZA  model  with  all  proposed  modihcations  improved 
the  predictions  for  all  surface  temperatures  except  at  =  1770  K.  Although  not 
shown,  the  primary  oxidation  product  again  becomes  CO2  for  the  modihed  ZA  model 
at  the  lowest  surface  temperatures.  Better  agreement  is  obtained  with  the  baseline 
ZA  model  at  =  1770  K  and  is  correlated  to  the  fact  that  CO  recombines  at  the 
surface  causing  lower  CN  mass  fractions  than  those  predicted  by  the  modihed  ZA 
model. 

The  CN  radiance  error  estimates  at  a  wall  normal  distance  of  0.15  mm  are  com¬ 
pared  in  Table  18.  The  ZA  model  with  nitridation  only  provided  a  5%  improvement 
compared  to  the  baseline  model  at  =  2410  K,  whereas  with  all  proposed  modih¬ 
cations,  the  prediction  was  improved  by  39%.  It  is  also  important  to  point  out  that 
the  radiance  predictions  for  the  ZA  model  with  all  proposed  modihcations  now  fall 
within  experimental  uncertainty  estimates  for  all  surface  temperatures  considered. 

Furthermore,  the  radiance  predictions  with  all  proposed  modihcations  now  display 
a  sensitivity  to  a  varying  surface  temperature  as  was  observed  in  the  experimental 
measurements.  It  was  originally  postulated  that  a  nitridation  mechanism  would  pro¬ 
duce  this  sensitivity  to  surface  temperature,  but  it  has  now  been  shown  to  be  due  to 
oxidation  mechanisms.  Specihcally,  controlling  the  ratio  of  CO  to  CO2  produced  at 
the  surface  has  caused  the  improved  predictions  and  surface  temperature  sensitivity. 
Increasing  the  amount  of  CO  in  the  boundary  layer  for  the  higher  surface  temperature 
cases  has  increased  the  CN  concentrations  and  thereby  increased  the  CN  radiance  pre¬ 
dictions.  Note  that  the  choice  of  gas-phase  chemistry  model  is  also  going  to  heavily 
influence  the  species  produced  at  the  surface  and  within  the  gas  itself. 
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Table  18.  CN  radiance  error  estimates  at  a  wall  normal  distance  of  0.15  mm. 


Snrface  Model 

Absolnte  Error 

=  2410  K 

ZA  (.lohnston) 

64% 

ZA  (ER  Nitridation) 

59% 

ZA  Modihcations 

25% 

=  2170  K 

ZA  (Johnston) 

50% 

ZA  (ER  Nitridation) 

48% 

ZA  Modihcations 

11% 

=  1920  K 

ZA  (Johnston) 

53% 

ZA  (ER  Nitridation) 

52% 

ZA  Modihcations 

34% 

=  1770  K 

ZA  (Johnston) 

1% 

ZA  (ER  Nitridation) 

2% 

ZA  Modihcations 

37% 

5.5  Heuristic  Modeling  Approach 


Instead  of  recommending  to  remove  a  reaction  from  the  original  ZA  model,  an 
investigation  was  performed  following  a  henristic  approach  that  attempted  to  synthe¬ 
size  the  oxidation  rates  to  more  closely  match  the  X-2  experimental  resnlts.  Fignre  59 
illnstrated  that  the  baseline  ZA  model,  which  also  inclnded  a  nitridation  mechanism, 
had  excellent  agreement  with  the  experiment  at  the  lowest  wall  temperatnre  of  1770  K. 
The  three  higher  wall  temperatnres  all  needed  more  CO  to  be  prodnced  at  the  snr- 
face,  which  was  adeqnately  snpplied  with  a  modihed  ZA  model  that  removed  the 
C02  +  (s)  •«-*•  C0  +  0(s)  reaction.  The  goal  of  the  investigation  was  to  take  advantage 
of  this  snrface  reaction  by  controlling  its  efficiency  at  different  temperatnres. 

The  stndy  was  execnted  nsing  the  30  km  condition  again  from  Section  4.1  becanse 
there  was  a  portion  of  the  snrface  on  the  sphere-cone  geometry  that  passed  throngh 
the  temperatnre  range  considered  in  the  X-2  experiments.  The  snrface  pressnres 
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over  this  surface  section  were  also  around  1  atm.  The  equilibrium  constant  for  the 
CO2  +  (s)  ^  CO  +  0(s)  reaction  was  less  than  one  within  the  temperature  range  of 
1770-2410  K,  so  from  a  pure  kinetic  standpoint,  this  reaction  prefers  to  proceed  in  the 
backwards  direction  as  similarly  noted  in  Table  17.  The  forward  rate  was  expressed  as 
an  ER  reaction  type  with  parameters  7er  =  0.9,  /9  =  0,  and  E  =  Q.  It  is  unclear  why 
Zhluktov  and  Abe  [136]  had  originally  dehned  the  reaction  in  this  manner  because 
other  work  has  dehned  it  in  its  preferred  direction  of  C0+0(s)  <->■  C02+(s)  [119,  123]. 
Hence  for  the  purposes  of  this  study,  the  CO-CO2  oxidation  reaction  is  recast  as 
CO  +  0(s)  •«-*•  CO2  +  (s)  dehned  with  similar  ER  rate  parameters  but  allowing  for 
adjustment  of  the  reaction  efficiency,  7er.  The  surface  equilibrium  constant  expression 
now  changes  to: 


K.  = 


X, 


CO2 


(126) 


A'coA'oA', 

where  the  gas-phase  equilibrium  constant  for  CO  -|-  O  — CO2  was  re-computed  using 
the  CEA  database.  The  corresponding  surface  species  production  rate  and  species 
mass  balance  equations  were  updated  to  rehect  the  change  in  reaction  direction. 

The  only  other  surface  reaction  that  produces  CO  is  0(s)  +  C{b)  ^  CO  -|-  (s), 
which  is  dehned  with  an  Arrhenius  rate  expression.  The  pre-exponential  Arrhenius 
coefficient  for  this  reaction  was  also  adjusted  as  part  of  this  study.  Lastly,  the  two 
additional  rate  modihcations  proposed  in  Section  5.3,  which  were  the  increase  in  ni- 
tridation  efficiency  and  decrease  in  carbon  condensation  efficiency,  were  also  retained. 

The  results  of  the  study  are  presented  in  Eigure  60  where  the  predicted  CO  sur¬ 
face  mass  hux  is  compared  to  the  baseline  ZA  model  with  nitridation,  the  proposed 
ZA  model  from  Section  5.3,  and  the  heuristically-determined  model.  The  tempera¬ 
ture  range  considered  in  the  X-2  experiments  is  denoted  by  the  solid  vertical  bars. 
The  new  model  now  exhibits  the  desired  behavior  of  CO  production  at  high  surface 
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2410  K  1770  K 


Figure  60.  Comparison  of  the  CO  surface  mass  flux  computed  by  different  ZA  model 
modifications  for  the  30  km  re-entry  condition. 

temperatures  and  CO  destruction  at  low  surface  temperatures.  Furthermore,  the  be¬ 
havior  of  the  baseline  ZA  model  is  recovered  at  surface  temperatures  below  1770  K. 
The  new  model  with  all  proposed  modihcations  is  presented  in  Table  19. 

The  four  wall  temperature  cases  from  the  hrst  X-2  campaign  were  again  rerun  with 
this  new  model  to  determine  if  the  CN  radiance  comparisons  were  further  improved. 
A  comparison  of  the  CO2  and  CO  species  surface  mass  fluxes  is  provided  in  Figures  61- 
64  using  modihed  forms  of  the  ZA  model  and  the  new  proposed  model.  Clearly,  the 
new  model  has  signihcantly  altered  the  dominant  oxidation  processes  at  the  higher 
surface  temperatures.  Production  of  CO2  has  become  negligible  at  the  two  higher 
surface  temperatures  and  begins  to  increase  in  appreciable  quantities  at  the  two  lower 
surface  temperatures.  The  CO  oxidation  mechanism  has  completely  reversed  at  the 
higher  surface  temperatures.  The  baseline  model  has  CO  recombining  at  the  surface, 
whereas  the  new  model  has  CO  as  the  primary  ablative  product.  As  the  temperature 
drops  below  2000  K,  the  new  model  predicts  CO  production  to  reduce  and  eventually 
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Table  19.  New  proposed  gas-surface  chemistry  model. 


Reaction 

Type 

So  /  7er.  /  A 

/9 

E  (kJ  mol  ^) 

1.)  0  +  (s)  e^O(s) 

Ads 

1 

0 

0 

2.)  N  +  {s)^  N{s) 

Ads 

1 

0 

0 

3.)  20{s)  ^  O2  +  2(s) 

Arrh 

3.58  X  10^3 

1 

256.07 

4.)  02  +  (s)  ^0  +  0(s) 

ER 

1 

0 

118.06 

5.)  CO  +  0(s)  ^  CO2  +  (s) 

ER 

1 

0 

0 

6.)  0(s)  +  C'(fe)  ^CO  +  {s) 

Arrh 

4.0  X  10^2 

1 

332.56 

7.)  0  +  0{s)  +  C{b)  ^  CO2  +  (s) 

ER 

0.8 

0 

16.63 

8.)  20{s)  +  C{b)^C02  +  2{s) 

Arrh 

3.58  X  10^^ 

0 

332.56 

9.)  C  +  (s)  ^  (s)  +  C'(6) 

ER 

0.216 

0 

0 

10.)  C2  +  2{s)  ^  2{s)  +  2C{b) 

ER 

0.5 

0 

0 

11.)  C3  +  3(s)  ^  3(s)  +  3C'(&) 

ER 

0.023 

0 

0 

12.)  N2  +  (s)  ^  N  +  N{s) 

ER 

1 

0 

636.85 

13.)  TV  +  (s)  +  (:;(&)  ^  civ  +  (s) 

ER 

0.396 

0 

36.86 
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Non-Dimensional  Distance  on  Surface  (x/r„) 

(a)  CO2 

Figure  61.  Predicted  CO  and  CO 
distance  from  the  stagnation  point  for 


(b)  CO 


surface  mass  fluxes  as  a  function  of  normalized 
2410  K  using  modifled  ZA  models. 


Non-Dimensional  Distance  on  Surface  (xirj 

(a)  CO2 

Figure  62.  Predicted  CO  and  CO 
distance  from  the  stagnation  point  for 


Non-Dimensional  Distance  on  Surface  (x/rj 


(b)  CO 


surface  mass  fluxes  as  a  function  of  normalized 
2170  K  using  modifled  ZA  models. 
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Non-Dimensional  Distance  on  Surface  (x/r„) 

(a)  CO2 

Figure  63.  Predicted  CO  and  CO 
distance  from  the  stagnation  point  for 


(b)  CO 


surface  mass  fluxes  as  a  function  of  normalized 
1920  K  using  modifled  ZA  models. 


(a)  CO2 


Figure  64.  Predicted  CO  and  CO 
distance  from  the  stagnation  point  for 


Non-Dimensional  Distance  on  Surface  (x/rj 


(b)  CO 


surface  mass  fluxes  as  a  function  of  normalized 
1770  K  using  modifled  ZA  models. 
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has  CO  reacting  at  the  surface  similar  to  the  behavior  predicted  by  the  baseline  model. 
The  solution  computed  with  the  new  ablation  model  has  met  the  desired  objectives 
of  this  investigation  by  displaying  the  appropriate  oxidation  processes  that  should 
improve  the  radiation  comparisons. 

The  CN  radiances  computed  with  the  new  ablation  model  are  compared  against 
the  prior  considered  ZA  models  in  Figure  65.  The  new  ablation  model  has  the  overall 
best  agreement  with  the  experiment  at  each  considered  surface  temperature.  The 
agreement  is  not  perfect,  but  by  appropriately  controlling  the  oxidation  ratio  of  CO 
to  CO2,  the  ablation  predictions  are  trending  in  the  right  direction.  The  CN  radiance 
error  estimates  at  a  wall  normal  distance  of  0.15  mm  for  the  new  model  now  fall 
within  a  range  of  7-25%  to  the  experiments  as  shown  in  Table  20.  The  error  range 
for  the  baseline  ZA  model  was  1-64%,  so  the  new  model  has  signihcantly  reduced 
the  error  probability.  The  new  model  represents  a  substantial  improvement  over  the 
baseline  ZA  model  and  should  be  considered  one  of  the  most  advanced  hnite-rate 
models  for  further  quantihcation. 

5.6  Second  Campaign  Results 

The  measured  radiance  values  displayed  a  non-monotonic  behavior  as  the  surface 
temperatures  were  increased  for  the  second  campaign.  As  mentioned  in  Chapter  II, 
initial  UV  spectrometer  measurements  imaging  a  broader  spectral  range  failed  to  mea¬ 
sure  any  signal  from  the  C3  Swings  and  C2  Swan  bands.  Therefore,  the  experimental 
results  and  numerical  comparisons  are  focused  on  the  CN  violet  band.  The  measured 
radiance  values  from  the  previous  four  lower  surface  temperature  cases  and  the  four 
higher  surface  temperature  cases  are  shown  in  Figure  66.  The  peak  radiance  consis¬ 
tently  occurs  in  the  near-wall  boundary  layer,  but  the  peak  value  drops  signihcantly 
going  from  a  surface  temperature  of  2410  K  to  2610  K.  In  fact,  the  four  higher  surface 
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(a)  =  2410  K  (b)  =  2170  K 


(c)T^  =  1920K  (d)T^  =  1770K 

Figure  65.  Comparison  of  CN  radiance  values  as  a  function  of  wall  normal  distance 
along  the  stagnation  streamline  for  the  lower  four  wall  temperature  cases  from  the  X-2 
experiments  using  the  new  proposed  model. 
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Table  20.  CN  radiance  error  estimates  at  a  wall  normal  distance  of  0.15  mm  with  the 
new  ablation  model. 


Surface  Model 

Absolute  Error 

=  2410  K 

ZA  (.Johnston) 

64% 

ZA  (ER  Nitridation) 

59% 

ZA  Modifications 

25% 

New  Model 

22% 

=  2170  K 

ZA  (.Johnston) 

50% 

ZA  (ER  Nitridation) 

48% 

ZA  Modifications 

11% 

New  Model 

7% 

=  1920  K 

ZA  (Johnston) 

53% 

ZA  (ER  Nitridation) 

52% 

ZA  Modifications 

34% 

New  Model 

24% 

=  1770  K 

ZA  (Johnston) 

1% 

ZA  (ER  Nitridation) 

2% 

ZA  Modifications 

37% 

New  Model 

25% 
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Figure  66.  Comparison  of  experimental  radiance  values  as  a  function  of  wall  normal 
distance  along  the  stagnation  streamline  for  surface  temperatures  from  the  first  and 
second  experimental  campaigns. 


temperature  cases  all  have  similar  radiance  profiles  except  for  the  3190  K  case,  which 
roughly  has  a  20%  lower  peak  value  compared  to  the  other  three  cases. 

The  unexpected  drop  in  measured  radiances  make  it  challenging  to  compare  the 
experimental  data  to  the  numerical  predictions.  Each  ablation  model,  including  the 
proposed  model  developed  in  Section  5.5,  displays  monotonic  behavior  by  predicting 
an  increase  in  radiance  with  an  increase  in  surface  temperature  as  shown  in  Figure  67. 
The  original  Park  ablation  model  is  not  shown  because  the  predicted  radiances  were 
nearly  twice  as  high  as  the  modified  Park  ablation  model.  All  predictions  shown  used 
the  Johnston  gas-phase  chemistry  model.  It  would  be  naive  to  make  a  fair  compar¬ 
ison  between  the  measurements  and  predictions  because  there  is  some  phenomena 
occurring  in  the  experiments  that  has  yet  to  be  quantified.  Therefore,  it  would  be 
illogical  to  conclude  that  there  is  good  agreement  at  surface  temperatures  of  2760  K 
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(a)  =  3280  K  (b)  =  3190  K 


(c)T^  =  2760K  (d)T^=2610K 

Figure  67.  Comparison  of  CN  radiance  values  as  a  function  of  wall  normal  distance 
along  the  stagnation  streamline  for  the  four  high  wall  temperature  cases  from  the  X-2 
experiments. 
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and  2610  K  using  the  ZA  model  with  and  without  nitridation.  A  proper  compari¬ 
son  is  not  possible,  so  the  analysis  is  limited  to  the  differences  between  the  model 
predictions. 

When  comparing  the  predicted  radiances,  the  new  ablation  model  is  consistently 
between  the  modihed  Park  and  ZA  with  nitridation  models.  The  predicted  results 
are  as  expected  due  to  the  ratios  of  CO2,  CO,  and  CN  that  are  generated  at  the 
surface.  The  new  model  primarily  has  CN  and  CO  production  whereas  the  ZA  model 
with  nitridation  has  CO2,  CO,  and  CN  production.  The  modihed  Park  model  also 
has  CN  and  CO  production  but  does  not  allow  for  C  to  recombine  at  the  surface. 
Recalling  the  analyses  from  the  previous  chapters,  these  combinations  of  species  yield 
the  expected  results  in  Figure  67. 

An  interesting  discovery  was  found  when  analyzing  the  high  speed  video  from  each 
test  case.  For  the  cases  with  a  surface  temperature  above  3000  K,  there  appear  to 
be  carbon  particles  spallating  off  of  the  surface.  The  particles  appear  after  the  how 
has  reached  a  “steady-state”  and  gradually  become  numerous  with  time.  Figure  68 
shows  the  difference  in  the  high  speed  camera  imagery  at  a  surface  temperature  of 
2610  K  and  3280  K.  At  3280  K,  an  immense  amount  of  tiny  carbon  particles  can 
be  seen  coming  off  the  model  surface.  For  ablating  environments,  both  shock  layer 
thermodynamics  and  composition  may  be  altered  by  the  presence  of  spalled  particles. 
As  was  shown  in  the  SEM  images  in  Chapter  II,  the  initial  carbon  model  surface  is 
amorphous  with  many  defects. 

It  is  postulated  that  the  hot  gases  are  occupying  these  defects  and  mechanically 
breaking  down  the  bulk  carbon.  The  effect  of  spallation  could  be  a  reason  why  the 
measured  radiances  were  lower  at  the  higher  surface  temperatures,  which  is  a  theory 
supported  by  the  work  of  Raiche  and  Driver  [101].  Raiche  and  Driver  experimentally 
verihed  the  presence  of  solid  particles  by  measuring  the  optical  attenuation  of  a  laser 
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through  an  ablating  arcjet  flow  [101],  Upon  removal  of  the  ablating  model,  the  laser 
transmission  reached  normal  levels.  They  also  found  that  the  attenuation  increased 
quadratically  with  heating  rate.  Therefore,  it  is  proposed  that  spallation  is  likely 
causing  the  reduction  in  CN  emission  intensities  at  the  higher  surface  temperatures 
due  to  line-of-sight  blockage  effects.  However,  to-date,  no  direct  measurements  have 
been  made  to  determine  spalled  particle  size  and/or  velocity  distributions. 
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(a)  =  2610  K 


Figure  68.  Comparison  of  high  speed  camera  images  showing  effect  of  spallation.  These 
images  have  been  enhanced  to  improve  visibility. 
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VI.  Conclusions 


Despite  the  prominence  of  carbon-based  materials  for  use  in  thermal  protection 
systems,  much  uncertainty  remains  in  predicting  thermochemical  ablation  rates  at 
high  surface  temperatures.  To  address  this  uncertainty,  a  series  of  experiments  using 
pre-heated  graphite  models  with  surface  temperatures  up  to  3280  K  were  conducted 
in  the  X-2  expansion  tunnel  at  The  University  of  Queensland.  Calibrated  shock- 
layer  emission  measurements  in  the  wavelength  region  353  to  391  nm  were  taken 
to  observe  the  effect  of  surface  temperature  on  radiation  from  the  CN  violet  bands. 
Attempts  were  also  made  to  observe  emissions  from  the  C3  Swings  and  C2  Swan 
bands,  however,  no  measurable  signal  was  found.  Most  of  the  analysis  focused  on  a 
surface  temperature  range  from  1770-2410  K  because  the  higher  surface  temperature 
cases  displayed  non-monotonic  behavior  with  ablative  phenomena  that  current  CFD 
codes  cannot  model. 

Numerical  simulations  were  performed  using  US3D  with  the  Park  and  Zhluktov 
and  Abe  hnite-rate  surface  kinetic  models.  The  gas-phase  chemistry  model  used 
legacy  Park  et  al.  [92,  93,  96,  97]  rates  and  updated  rates  from  Johnston  et  al.  [46]. 
It  was  found  that  the  Johnston  rates  mostly  influenced  the  mass  fluxes  and  concen¬ 
trations  of  CO  and  CO2  but  did  not  signihcantly  affect  the  CN  radiation  spectra. 
Results  were  also  analyzed  using  a  modihed  nitridation  rate  proposed  by  Suzuki  et 
al.  [115]  in  the  Park  ablation  model,  which  signihcantly  reduced  the  amount  of  CN 
formed  at  the  surface.  The  simulation  results  were  applied  in  NEQAIR  to  reproduce 
the  experimental  radiance  prohles. 

The  ablation  products  were  quite  different  between  the  surface  kinetic  models. 
For  the  oxidation  mechanisms,  the  Park  model  only  produced  CO  whereas  the  ZA 
model  produced  CO2  and  had  CO  recombining  at  the  surface.  The  Park  model 
includes  a  nitridation  reaction  with  an  efficiency  that  does  not  replicate  experimental 


164 


results,  including  the  results  presented  in  this  work.  By  reducing  the  Park  nitridation 
efficiency,  good  agreement  was  obtained  in  comparing  measured  CN  violet  radiance 
values  at  a  surface  temperature  of  2410  K.  The  ZA  model  did  not  show  good  agreement 
at  this  surface  temperature,  and  it  was  hrst  assumed  to  be  due  to  a  lack  of  a  nitridation 
reaction. 

As  the  surface  temperature  was  decreased,  the  ZA  model  radiance  predictions 
improved  greatly.  In  fact,  the  measured  near  surface  radiance  at  =  1770  K  was 
accurately  predicted  by  the  ZA  model  with  the  .Johnston  gas  chemistry  model.  The 
two  intermediate  surface  temperatures  of  1920  and  2170  K  had  measured  radiance 
prohles  that  were  in  between  the  predictions  of  the  Park  and  ZA  models.  Logically, 
this  led  to  the  assumption  that  nitridation  became  a  relevant  surface  reaction  at  wall 
temperatures  above  1800  K. 

Modihcations  to  the  ZA  model  were  proposed  based  on  experimental  and  theoret¬ 
ical  results  with  the  goal  of  improving  the  CN  radiance  predictions.  Two  approaches 
for  modeling  carbon  nitridation  were  developed  with  rate  parameters  derived  from 
experimental  measurements  and  theory.  One  approach  assumed  that  atomic  nitrogen 
could  react  directly  with  the  bulk  carbon  to  form  gaseous  CN.  The  other  approach 
followed  a  more  traditional  gas-surface  interaction  process  and  required  atomic  ni¬ 
trogen  to  be  in  an  adsorbed  surface  state  before  reacting  with  the  bulk  carbon.  The 
experimental  test  conditions  from  which  the  rate  parameters  were  deduced  had  rel¬ 
atively  low  wall  temperatures  as  compared  to  those  that  may  be  experienced  during 
a  re-entry  flight.  Therefore,  additional  measurements  at  higher  surface  temperatures 
are  needed  to  improve  the  nitridation  rate  parameters  derived  in  this  study. 

The  nitridation  reaction  should  be  included  in  any  air-carbon  system  because 
it  will  have  an  effect  on  surface  and  radiative  heating.  Yet,  its  importance  varies 
with  surface  temperature.  It  was  concluded  that  a  direct  nitridation  mechanism 
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has  the  most  signihcant  impact  on  predicted  surface  mass  fluxes  and  species  mass 
fractions.  Direct  nitridation  was  active  between  surface  temperatures  ranging  between 
approximately  2000-4000  K,  and  the  CN  production  rate  dropped  or  approached  zero 
outside  of  this  range.  The  surface  participating  nitridation  reaction  was  shown  not 
to  influence  any  of  the  mass  fluxes  or  the  composition  of  the  boundary  layer  at  the 
stagnation  point. 

The  large  differences  between  these  nitridation  mechanisms  are  due  to  many  fac¬ 
tors  in  combination.  Under  conditions  where  desorption  rates  are  high,  such  as  at 
high  surface  temperatures,  surface  coverage  would  be  low.  Hence,  the  surface  par¬ 
ticipating  nitridation  reaction  would  be  much  less  likely  to  occur  compared  to  direct 
nitridation.  Furthermore,  direct  nitridation  would  be  uninhibited  because  the  surface 
would  be  entirely  empty.  Then  as  the  wall  temperature  cools,  there  is  a  competition 
between  the  adsorbed  atomic  oxygen  and  nitrogen  sites.  However,  if  the  wall  tem¬ 
perature  cools  too  rapidly,  it  may  result  in  the  nitridation  rate  becoming  inactive. 
There  is  a  balance  occurring  between  the  active  sites  with  adsorbed  atomic  nitrogen 
and  wall  temperature  that  will  ultimately  determine  the  nitridation  rate.  Note  that 
there  is  nothing  to  prevent  implementing  both  nitridation  mechanisms  into  a  kinet¬ 
ics  model  because  both  are  plausible  surface  reaction  processes.  Direct  nitridation 
was  chosen  for  implementation  into  the  new  ablation  model  because  it  predicted  sub¬ 
stantial  CN  production  at  high  surface  temperatures  that  agreed  with  equilibrium 
ablation  modeling. 

Additionally,  a  sensitivity  analysis  of  the  ZA  model  was  performed  to  determine 
which  surface  reactions  had  the  most  impact  on  CN  radiance  predictions.  It  was  found 
that  increasing  the  amount  of  CO  formed  at  the  surface  increased  CN  concentrations. 
Examining  the  surface  reaction  fluxes  showed  that  CO2  was  the  dominant  oxidation 
product,  contrary  to  existing  experimental  evidence  that  measured  primarily  CO  at 
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high  surface  temperatures.  Additionally,  an  equilibrium  air-carbon  mixture  indicated 
that  a  gas  would  be  composed  of  only  N2  and  CO  within  a  temperature  range  of 
1500-2500  K.  A  series  of  rate  modihcations  were  proposed,  along  with  the  initial 
recommendation  to  remove  the  CO2  +  (s)  CO  -|-  0(s)  surface  reaction,  and  then 
applied  to  the  X-2  experiments.  The  CN  radiance  predictions  improved  by  about 
19-39%  at  surface  temperatures  of  1920,  2170,  and  2410  K  compared  to  the  baseline 
ZA  model.  At  a  surface  temperature  of  1770  K,  the  modihed  ZA  model  over  predicted 
the  radiance  by  about  37%.  Overall,  the  proposed  modihcations  vastly  improved  the 
predictions  and  exhibited  a  sensitivity  to  a  varying  surface  temperature,  which  was 
lacking  with  the  baseline  model. 

The  CO2  +  (s)  <->■  CO  +  0(s)  reaction  included  in  the  baseline  ZA  model  is  a  valid 
surface  reaction,  so  instead  of  completely  removing  it,  an  attempt  was  made  to  take 
advantage  of  this  reaction  by  controlling  its  efficiency  at  different  temperatures.  It  was 
proposed  to  specify  this  reaction  in  the  reverse  order  because  it  preferred  to  proceed 
in  the  backwards  direction  under  these  simulated  conditions.  Prior  experimental 
work  had  also  dehned  this  reaction  as  CO  +  0{s)  ^  CO2  +  (s).  Therefore,  if  the 
forward  rate  for  this  reaction  is  quantihed  in  future  ablation  experiments,  it  would  be 
advantageous  to  have  the  numerical  forward  rate  specihed  in  the  same  manner  for  an 
apt  comparison.  The  Arrhenius  rate  for  the  other  oxidation  reaction  that  produced 
CO,  0(s)  -|-  C{b)  ^  CO  -|-  (s),  was  also  modihed  as  part  of  the  modeling  effort. 
The  developed  ablation  model  that  included  these  oxidation  modihcations,  and  those 
from  the  sensitivity  analysis,  further  improved  the  CN  radiation  comparisons.  The 
primary  oxidation  product  became  CO  at  the  higher  surface  temperatures  and  then 
switched  to  CO2  at  the  lower  surface  temperatures.  All  radiation  predictions  fell 
within  the  experimental  error,  and  the  error  probability  range  dropped  from  1-64% 
for  the  baseline  ZA  model  to  7-25%  for  the  new  model. 
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The  experimental  results  from  the  second  campaign  were  not  as  anticipated  be¬ 
cause  there  was  a  signihcant  drop  in  CN  radiative  intensities.  It  was  speculated  that 
the  decrease  in  measured  radiance  was  due  to  an  attenuation  of  the  optical  signal 
caused  by  spalled  carbon  particles.  Evidence  from  the  high  speed  video  supported 
this  hypothesis,  but  it  has  yet  to  be  quantihed.  Unfortunately,  a  comparison  of  the 
numerical  predictions  with  the  measured  radiances  was  practically  impossible.  How¬ 
ever,  the  predictions  behaved  as  expected  and  showed  an  increase  in  CN  radiance 
with  an  increase  in  surface  temperature. 

It  also  must  be  acknowledged  that  the  error  in  the  experimental  data  could  come 
from  multiple  sources.  An  attempt  to  quantify  this  error  was  made  by  approximating 
an  uncertainty  range  on  the  measured  radiances.  However,  these  were  truly  unprece¬ 
dented  experiments  and  some  measurement  techniques  had  to  be  adjusted  due  to 
unanticipated  difhculties.  These  techniques  will  be  improved  upon  as  testing  at  these 
extreme  conditions  becomes  more  frequent. 

Clearly,  much  work  remains  in  understanding  the  surface  kinetic  mechanisms  of  an 
air-carbon  system.  The  work  presented  in  this  dissertation  has  made  a  substantial  leap 
forward  in  determining  some  of  the  proper  mechanisms  at  relevant  hypersonic  re-entry 
conditions.  The  current  state-of-the-art  ablation  models  were  shown  to  be  inadequate 
to  accurately  predict  species  mass  fluxes  that  correlated  to  experimental  radiation 
measurements.  Furthermore,  the  models  did  not  even  match  experimental  trends 
when  surface  temperatures  were  increased.  The  Park  model  was  found  to  be  too 
simple  in  that  it  only  specihed  one-step  kinetic  processes  and  could  not  capture  non¬ 
linear  oxidation  processes.  The  ZA  model  had  competing  surface  kinetic  processes 
and  attempted  to  account  for  non-linear  Arrhenius  behavior.  However,  the  ZA  model 
had  major  dehciencies.  A  majority  of  the  novel  research  presented  in  this  dissertation 
was  driven  by  improving  the  surface  kinetic  processes,  starting  with  the  ZA  model  as 
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a  baseline,  to  better  match  experimental  resnlts.  The  new  model  that  was  developed 
has  been  shown  to  vastly  improve  the  comparison  with  experimental  resnlts.  The  new 
model  now  follows  the  trend  in  the  experimental  data  where  the  primary  oxidation 
prodnct  appears  to  switch  from  CO2  to  CO  as  the  snrface  temperatnre  rises  above 
approximately  1800  K.  Moving  forward  with  this  research,  it  is  recommended  to  nse 
the  new  model  for  fntnre  carbon  ablation  simnlations  instead  of  either  the  Park  or 
ZA  models. 

As  research  into  carbon  ablation  proceeds,  it  wonld  be  constrnctive  to  consider 
the  nse  of  higher  hdelity  mass  diffnsion  models  becanse  Tick’s  law  typically  provides 
conservative  heating  rate  predictions  at  hypersonic  re-entry  conditions.  Using  Tick’s 
law  with  a  constant  Lewis  nnmber  does  not  strictly  enforce  that  the  snm  of  the  mass 
diffnsion  flnxes  is  zero.  There  is  a  modihed  form  of  Tick’s  law  that  ensnres  total  dif¬ 
fnsion  flnx  of  zero,  and  this  form  has  been  shown  to  predict  flowheld  properties  that 
are  in  good  agreement  with  mnlti-component  diffnsion  models.  The  cnrrent  US3D 
code  does  not  have  a  version  with  modihed  Tick’s  law  implemented.  Additionally, 
implementing  the  modihed  form  of  Tick’s  law  wonld  reqnire  a  re-derivation  of  the 
snrface  mass  balance  eqnation  for  each  gas  species.  With  more  time  and  compnta- 
tional  resonrces,  the  incorporation  of  a  higher  hdelity  diffnsion  model  wonld  be  the 
next  step  in  the  research  process. 

It  was  shown  that  the  CO2  +  (s)  ^  CO  -|-  0(s)  reaction  had  a  signihcant  effect 
on  the  resnlting  radiation  predictions.  At  high  snrface  temperatnres  (greater  than 
1800  K),  CO  is  the  primary  oxidized  species.  At  relatively  low  snrface  temperatnres 
(less  than  1800  K),  CO2  may  be  the  primary  oxidized  species.  Therefore,  it  is  desired 
to  have  fntnre  carbon  ablation  experiments  analyze  the  oxidation  processes  more 
closely  to  fnrther  qnantify  the  newly  developed  ablation  model. 

Tnrther  investigation  is  warranted  to  nnderstand  why  mnch  more  CN  radiation 
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was  measured  in  the  shock  layer  than  what  was  predicted  by  the  simulations.  It  is 
recommended  that  the  X-2  nozzle  exit  conditions  be  calculated  as  an  axisymmetric 
simulation  including  nonequilibrium  gas  chemistry  effects.  The  shock  layer  radiation 
discrepancies  could  be  quantihed  with  higher  hdelity  nozzle  exit  conditions.  An  initial 
understanding  could  be  achieved  by  artihcially  adding  more  carbon  to  the  freestream 
mass  fractions  until  shock  layer  CN  concentrations  reach  those  observed  in  the  X-2 
experiments. 

Lastly,  future  experiments  need  to  quantify  the  influence  of  spallation  on  measured 
radiances.  It  is  only  speculation  that  spalled  carbon  particles  are  attenuating  the 
optical  path  to  the  spectrometer  camera,  thereby  reducing  the  CN  violet  radiative 
intensities.  If  discovered  to  be  true,  then  a  change  in  model  material  should  be 
considered  to  lessen  the  chance  for  spallation  to  occur.  For  example,  instead  of  using 
an  amorphous  graphite,  a  carbon  hber  material  may  be  more  functional.  A  uniform 
grain  structure  may  make  it  more  resistant  to  breaking  apart  under  extreme  heating 
conditions. 

The  improved  ablation  modeling  capability  presented  in  this  dissertation  may  one 
day  aid  in  the  development  of  innovative  TPS  design  to  ensure  the  survivability  and 
lethality  of  high  speed  strike  weapons  or  penetrating  cruise  missiles.  Alternatively, 
these  modeling  capabilities  may  serve  to  mature  heat  shield  materials  for  reusable 
space  access  vehicles.  At  a  minimum,  the  methods  used  to  characterize  shock  layer 
radiation  can  be  implemented  immediately  into  the  IC  for  the  purpose  of  identifying 
hostile  threats.  In  that  case,  the  hdelity  of  these  methods  could  mean  the  difference 
between  a  successful  interception  or  a  Pearl  Harbor  type  of  event. 
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